HILTON ANDRADE DE MELLO 
RONALDO SÉRGIO DE BIASI 


INTRODUÇÃO A 
FÍSICA 
DOS 
SEMICONDUTORES 


INTRODUÇÃO A FÍSICA 
DOS SEMICONDUTORES 


FICHA CATALOGRÁFICA 


(Preparada pelo Centro de Catalogação-na-fonte, 
Câmara Brasileira do Livro, SP) 


Melo, Hilton Andrade de. 


Introdução à física dos semicondutores |por | Hilton 
Andrade de Mello [e [Ronaldo Sérgio de Biasi. São Paulo, 
Edgard Blúcher; Brasília, INL, 1975. 

p. ilust. 


Bibliografia. 


1. Semicondutores |. Biasi, Ronaldo Sérgio de, 1943-— 
HI. Instituto Nacional do Livro. Ill. Título. 


CDD: 537.622 
CCF/CBL/SP—75-1026 CDU: 536.311 


Índice para catálogo sistemático (CDD): 
1. Semicondutores: Física 537.622 


Eng.º HILTON ANDRADE DE MELLO 


Engenheiro Eletrônico e Nuclear; Eng.º pela Stanford University, 
California, EUA; Chefe do Laboratório de Instrumentação e 
Controle do Instituto de Engenharia Nuclear 


Ph. D. RONALDO SÉRGIO DE BIASI 


Engenheiro Eletrônico; Ph. D. pela University of Washington, 
Washington, EU A; Professor de Fisica do Estado Sólido 
do Instituto Militar de Engenharia 


INTRODUÇÃO À FÍSICA 
DOS SEMICONDUTORES 


editora EDGARD BLUCHER Ida. 


EM CONVÊNIO COM O INSTITUTO NACIONAL DO LIVRO— 
— MINISTÉRIO DA EDUCAÇÃO E CULTURA 


01975 Editora Edgard Bliicher Ltda. 


Este livro foi coeditado com o Instituto Nacional do Livro/Mi- 
nistério da Educação e Cultura. dentro do Programa do Livro 
Didático-Ensino Superior, patrocinado pela Secretaria de Plane- 
jamento da Presidência da República. 


É proibida a reprodução total ou parcial 
por quaisquer meios 
sem autorização escrita da editora 


EDITORA EDGARD BLÚCHER LTDA. 


0 1000 Caixa PostAL 5450 — Rua Perxoro Gomipe, 1400 
En. TeLecrÁrico: BLUCHERLIVRO — FONES (011)287-2043 E 288-5285 
São PauLo — SP — BrasiL 


Impresso no Brasil Printed in Brazil 


A nossas esposas 


MARIA DE PAULA 
e 


MARÍLIA 


Conteúdo 


DALEOCLT ER A Gs NASA a ia asi e ibee paira ga sirva Bi cis ESDE IX 
Constantes ECA E od, SIL IL AO C RECESSO TE DR ITS XI 
Parâmetros de alguns semicondutores a 27ºC ......cccicciccsscicios XI 
1. Considerações clássicas sobre a estrutura atômica ................ 1 
2. Noções básicas de mecânica quântica ,,...ccciseceserarestaateaaa 8 
3. Aplicações da mecânica quântica .......ccccciiisccciiiccicrees 13 
é. Noções de Cristalogralr esa se gi E ES RÃS pp SDS q PU 27 
5. Aplicação da mecânica quântica a um cristal ......cccccilciicoioo 42 
6. Propriedades básicas dos materiais semicondutores ................ 49 
do ENSHIDUIÇÕES: CATALÍNTICAS cosmos pi pres DEE CMT AE Nai PT Pe 61 
8. Aplicação da distribuição de Fermi-Dirac a um material semicondutor 70 
9. Transporte de cargas nos semicondutores .........ccccciccicccso. 83 
TO dTunções- PN mes su virus lost nine canso apare o a aloja a prelo adoro tato jeiaa piada 1 93 
Li, Eransistores de: JONÇÃO: e pasta ira praia isso pia gre GEES ia À) SACAS trad; prsoteag 107 
12. Outros fenômenos observados nos semicondutores ............c... 114 
Apêndices 
AS "Ecóranã te; BIVOR ss ANAL LISAS E CITESA BRL A DIA 119 
B. Experiências em semicondutores ........cccccicccierescercrseero 121 


Introdução 


Desde a invenção do transistor, o estudo da Eletrônica tem-se concentrado 
cada vez mais no projeto e utilização dos dispositivos semicondutores. Em 
consegiiência, passou-se a exigir do estudante de Eletrônica um conhecimento 
das propriedades básicas da matéria que até então não era realmente necessário 
para sua profissão. 


Ao mesmo tempo, o estudo dos materiais semicondutores foi incluído no 
currículo das escolas de Física, como uma aplicação de grande interesse prático 
dos conceitos abstratos da Física do Estado Sólido. 

Neste livro, que foi escrito tendo em vista tanto os estudantes de Eletrônica 
como os de Física, procuramos relacionar as propriedades básicas dos dispo- 
sitivos semicondutores aos conceitos fundamentais da Mecânica Quântica e da 
Física do Estado Sólido. 


O Cap. | é uma introdução histórica à teoria eletrônica da matéria; as 
noções fundamentais de Mecânica Quântica são apresentadas no Cap. 2. 


No Cap. 3 são discutidas duas aplicações importantes da Mecânica Quân- 
tica: o estudo da estrutura eletrônica dos átomos e o estudo do Gás de Fermi, 


Os conceitos de cristalografia necessários aos estudo dos sólidos são expostos 
no Cap. 4; no Cap. 5, estudamos a aplicação dos princípios de Mecânica Quân- 
tica a um sólido cristalino, surgindo naturalmente o conceito de faixas de energia 
permitidas e faixas proibidas. 


O Cap. 6 é dedicado ao estudo das propriedades básicas dos materiais 
semicondutores; nesse capítulo são definidos e discutidos os conceitos de ligação 
covalente, buraco, recombinação, impureza, energia de ionização e massa efetiva. 

As distribuições estatisticas, essenciais para a compreensão dos fenômenos 
coletivos nos sólidos, são cxaminadas no Cap. 7; no Cap. 8, a distribuição de 
Fermi-Dirac é aplicada a um material semicondutor, como meio para calcular 
as concentrações de elétrons e buracos. 


O Cap. 9 trata dos fenômenos de transporte nos semicondutores, introdu- 
zindo os conceitos de mobilidade e difusão. 

Os Caps. 10 e 11 mostram que os conceitos desenvolvidos nos capítulos 
anteriores podem ser aplicados diretamente aos dois dispositivos semicondu- 
tores mais usados na prática: os diodos semicondutores (Cap. 10) e os tran- 
sistores de junção (Cap. 11). 

O Cap. 12 trata de alguns fenômenos importantes que não foram discutidos 
nos capítulos anteriores: Luminescência, Efeito Túncl c Efeito Hall. 

O Teorema de Bloch, essencial para a compreensão dos fenômenos de 
condução nos sólidos, é demonstrado no Apêndice A; o Apêndice B é uma lista 
de referências a experiências simples que ilustram os fenômenos estudados 
neste livro, 


No final de cada capítulo, o leitor encontrará uma série de exercícios. As 
respostas ou sugestões para alguns exercicios — assinalados com um asterisco — 
se encontram no final do livro. 

Queremos expressar nossos agradecimentos a todos que colaboraram 
conosco, em especial ao fisico Luiz Telmo Auler por seus comentários e suges- 
tões, ao Sr. Luis Antonio Costa Ferreira pela cuidadosa execução das ilustrações. 
e à Sr“. Eneida A. Mendonça pelo excelente trabalho de datilografia dos originais. 


H.A.M. 
R.S. B. 


CONSTANTES FÍSICAS 


Símbolo À 


Gerais 

C 

L 

No 

k 
Atômicas 
h 


h 
Ed 


Elétricas 
o) 
Hg 


Partículas 


Nome 


velocidade da luz no vácuo 
número de Loschmidt 
número de Avogadro 
constante de Boltzmann 


constante de Planck 
h/2x 
raio de Bohr 


permissividade do vácuo 
permeabilidade do vácuo 


massa de repouso do elétron 
massa de repouso do próton 
massa de repouso do néutron 
valor absoluto da carga do clétron 


constante dos gases perfeitos 


Valor e unidade 


2,998 x 10º m/s 
2.69 x 107º m ? 
6.0225 x 107º mol” ' 
1,381 x 1023 J/K 


6,626 x 1073*]-s 
1,054 x 1073 J-s 
5,292 x 107 !!'m 


8,854 x 10712 F/m 
1,257 x 107º H/m 


9,108 x 103 kg 
1,672 x 10727 kg 
1,675 x 10-27 kg 
1,601 x 101º C 


8.31 J/mol-K 


PARÂMETROS DE ALGUNS SEMICONDUTORES A 27º€ 


Densidade (átomos/cm*) 
Constante da rede (nm) 
Massa especifica (g/em?) 
Largura da faixa proibida (eV) 


Concentração intrínseca 
(portadores/cm'*) 
Mobilidade (cm?/V «s): 
Elétrons 
Buracos 
Constante dielétrica 
Densidade eletiva de estados: 
Faixa de condução: N.(cm”*) 
Faixa de valência: N, (cm?) 


Ge 


4,42 x 1072 
0,556 

5,32 

0,67 


2,40 x 10!? 
3900 

1900 

16,3 


1,04 x 101º 
6,00 x 1018 


Si 
5,00 x 1022 
0,543 
2,33 
RE 


1,45 x 10'º 
1350 

480 
11,7 


2,80 x 101º 
1,04 x 101º 


GaAs 


2,21 x 1022 
0,565 

5,32 

1,40 


9,00 x 10º 
8 600 

250 
12,0 


4,70 x 10"? 


7,00 x 1038 


1 Considerações clássicas sobre 
a estrutura atômica 


No início deste século, não havia unanimidade entre os cientistas a respeito 
da estrutura dos átomos. Thomson (1907)* propôs um modelo segundo o qual 
as partículas positivas e negativas estariam uniformemente distribuidas em 
todo o átomo. 

Experiências realizadas com particulas alfa (núcleos de hélio) emitidas 
por uma substância radioativa mostraram que essas particulas sofrem um desvio 
quando atravessam uma lâmina fina de ouro ou prata. O desvio é chamado de 
espalhamento. De acordo com o modelo de Thomson, o ângulo de espalhamento, 
(7. deveria ser sempre menor que 90º, como na Fig. 1.1 (a). 

Entretanto, as experiências de Geiger e Marsden (1909) mostraram que 
algumas das partículas alfa sofrem um desvio maior que 90º, como na Fig. 1.1 (b). 


Figura 1.1 Espalhamento de particulas alfa: 
(a) ângulo /Z < 90º: (b) ângulo (5 > 90º 


Para explicar esse fato, Rutherford (1911) sugeriu que as partículas que 
constituem o átomo não estão uniformemente distribuidas, mas que toda a 
carga positiva está concentrada em uma região central (núcleo), enquanto que 
a carga negativa está distribuida em torno da região central. 

Posteriormente, verificou-se que a região central do átomo (núcleo) é cons- 
tituída por dois tipos diferentes de partículas, o próton, que possui uma carga 
elétrica positiva, e o nêutron, que não possui carga elétrica. As cargas negativas 
distribuidas em torno do núcleo foram associadas às partículas chamadas elétrons, 
que possuem carga elétrica negativa. As cargas e massas dos três tipos de parli- 
culas estão indicadas na Tab. 1.1. 


“No final do capitulo o leitor pode encontrar a referência completa 
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Tabela 1.1 
— Partícula Massa Carga elétrica 
Elétron m, = 9,1083 x 102! kg 1,609 x 107*º coulomb 
Nêutron M, = 1838,65 m, 0 
Próton M, = 1836,12m, + 1,609 x 107!º coulomb 


O átomo em seu estado normal é eletricamente neutro; como a carga do 
próton tem o mesmo valor numérico que a carga do elétron, o número de elétrons 
em um átomo deve ser igual ao número de prótons no núcleo, Esse número é 
chamado de número atômico e representado pela letra Z. A soma do número 
de prótons e do número de nêutrons é chamada de número de massa e repre- 
sentada pela letra M, Sem nos preocuparmos com a distribuição dos elétrons 
em torno do núcleo, podemos representar esquematicamente um átomo como 
na Fig. 1.2. 


Eletrons 


Núcleo 


Z Prótons 
M-Z Nêutrons 


Figura 1.2 Representação esquemática de um 
átomo 


Nosso próximo passo será o estudo da distribuição dos elétrons em torno 
do núcleo. Antes, porém, convém mencionar um outro fenômeno, estudado por 
Planck, cuja explicação serviu de base para a formulação de duas importantes 
teorias: a teoria da emissão fotoelétrica de Einstein e a teoria do espectro do 
hidrogênio de Bohr. Trata-se do estudo da distribuição da energia irradiada 
por um corpo negro'* em função do comprimento de onda da radiação (Rayleigh, 
1900 — Jeans, 1905). Quando o problema do corpo negro é resolvido utilizando-se 
a teoria eletromagnética clássica, chega-se à conclusão de que a intensidade da 
radiação do corpo negro deve variar inversamente com a quarta potência do 
comprimento de onda. Na Fig. 1.3 vemos os resultados previstos usando a teoria 
eletromagnética clássica e os resultados experimentais. 

A fórmula de Rayleigh-Jeans concorda com os dados experimentais para 
grandes comprimentos de onda, mas conduz a resultados absurdos para pequenos 
comprimentos de onda. Planck (1901) verificou que admitindo que o corpo 


(Corpo negro é definido como qualquer corpo cuja emissividade é igual a 1, isto é, qué irradia 
energia com o máximo possivel de eficiência. O corpo negro também tem um fator de absorção uni- 
tário, isto é, absorve toda a energia que incide sobre ele 
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N Lei de 
a Rayleigh-Jeans 


T4>73>7T2>% 


ENERGIA 


COMPRIMENTO DE ONDA 


Figura 1.3 Distribuição espectral da energia irradiada pelo corpo negro. A teoria ele- 
tromagnética clássica conduz à Lei de Rayleigh-Jeans, segundo a qual R, = 27 ckTA* 


negro não irradia continuamente, como a teoria eletromagnética clássica prevê, 
mas em quanta'* de encrgia, cada quantum com uma energia hv onde v é a 
frequência da radiação e h uma constante escolhida convenientemente, era 
possível obter uma boa concordância com os dados experimentais. Planck 
atribuiu a h o valor de 6,55 x 107?” erg-s. 

É evidente que uma proposição arbitrária como a de Planck, que contra- 
riava a estabelecida teoria eletromagnética de Maxwell, não teve grande recep- 
tividade por parte dos cientistas da época. Entretanto, Einstein (1905) utilizou 
a hipótese de Planck para explicar com grande sucesso o fenômeno da emissão 
fotoelétrica, que consiste na emissão de elétrons por um metal quando este é 
exposto a uma fonte luminosa'**. 

Passemos agora ao estudo do movimento dos elétrons em torno do núcleo. 
Como o elétron possui uma carga elétrica e está em movimento, deveria, de 
acordo com a teoria eletrodinâmica clássica, irradiar energia continuamente. 
Em conseqiiência, não existiriam átomos estáveis, ao contrário do que se observa 
na prática. Essa contradição permaneceu até que Bohr (1913, 1915), contra- 
riando as leis da eletrodinâmica clássica, propôs que o elétron gira em torno 
do núcleo em órbitas estáveis, nas quais não emite energia eletromagnética. 
Bohr propôs ainda que as órbitas nas quais o elétron não irradia são órbitas 
estacionárias, tais que o momento angular *** do elétron é igual a um número 

"Quanta é o plural de quantum 

'**De acordo com Einstein, o fenômeno fotoelétrico é representado pela equação: hv = mv?/2 + 
+ à, onde hy é o quantum da luz, mv?/2 a energia cinética com que o elétron é emitido e p é a função 
trabalho do material, que corresponde à energia necessária para deslocar o elétron até a superficie 
do metal 

1***De um modo geral, chamamos de momento angular o momento da quantidade de movimento 


(m- 2) em relação ao ponto considerado (definido pelo raio vetor r): momento angular = m-v x F 
(Nessa expressão x significa produto vetorial) 
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inteiro de h/2x, onde h é a mesma constante usada por Planck: 
momento angular = mvr = nh/2n 


onde n é um número inteiro. 

Finalmente, Bohr introduziu outra hipótese, baseada nos estudos de Planck. 
Ele supôs que quando o elétron passa de uma órbita de energia E, para uma 
órbita de energia E, (E, > E,), irradia energia com uma fregiiência tal que a 
relação de Planck é satisfeita, isto é, E, —-E, = hv. 

Vamos aplicar as hipóteses de Bohr ao caso simples do átomo de hidro- 
gênio, que possui apenas um elétron. Suponhamos que esse elétron está descre- 
vendo uma trajetória circular em torno do núcleo (Fig. 1.4). 


Figura 1.4 Representação esquemática do átomo 
de hidrogênio 


Núcleo 


Para que a órbita seja estável. a força de atração que o núcleo exerce sobre 

o elétron deve ser equilibrada pela força centrifuga resultante do movimento 
de rotação; isto é!*), 

E 


meo Po Fr. EM. de] 


onde 
q é a carga do elétron, 
me a massa do elétron, 
v é a velocidade tangencial do elétron, 
r é o raio da órbita. 


A energia total do elétron é a soma da energia potencial e da energia ciné- 
tica: isto é, 
E = Eç+ E =5——— É (Eq. 1.2) 
Vamos agora introduzir a restrição de que o momento angular deve ser 
um número inteiro de h/2n: 
mvr = nh/2x (Eq. 1.3) 


“O leitor deve observar que estamos usando o sistema MKS racionalizado, no qual a força 
de atração entre duas cargas q, e q; é dada por F = q,q,/4nsyr”. Caso estivéssemos usando o sis- 
tema não-racionalizado, teriamos 
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De acordo com as Egs. 1.1, 12 e 1.3, temos 


PE pi a 14 
E =Seinho as 299 
Para duas órbitas 1 e 2, as energias serão dadas por 
4 4 
= — a ' O (ES sa air , 
Ei = Fendi É = Beni (Bg, 15) 


onde os números inteiros n, e n, caracterizam as duas órbitas. 
A diferença de energia entre as duas órbitas é dada por 


4 
E -5) E 
AE = Es E=seir (a nz = hy, (Eq. 1.6) 
Finalmente, a frequência v da radiação é dada por 
o mg 1d), - 
V= CT (az 2) (Eq. 1.7) 


O resultado anterior, obtido por Bohr, está de pleno acordo com os resul- 
tados experimentais. Na verdade, os espectroscopistas já haviam observado o 
espectro do hidrogênio, chegando a fórmulas empíricas que permitiam prever 
as radiações monocromáticas emitidas. Essas fórmulas eram usadas para cal- 
cular, não o comprimento de onda ou a fregiência, mas o chamado nimero 
de onda (frequência da radiação dividida pela velocidade da luz). Rydberg, por 
exemplo, mostrou que as séries observadas por Lyman, Balmer e Paschen podiam 
ser enquadradas na fórmula geral seguinte: 

; | l 
número de onda = R (nm N ) (Eq. 1.8) 
onde M e N são números inteiros e R é a chamada constante de Rydberg (R = 
= [09,737 cm” !). 

Em termos de números de onda, a expressão obtida por Bohr (Eq. 1.6) 

pode ser escrita: 


à ade fl, 
número de onda = ETR (ma ) (Eq. 1.9) 
ou 
E dr 
número de onda = R (1-5) (Eq. 1.10) 
1 2 
onde: 
4 
ro mg 
R Seghºc 


A constante R' que aparece na Eq. 1.10 concorda excepcionalmente bem 
com o valor experimental (R) obtido por Rydbereg. 

Na Fig. 1.5 ilustramos algumas transições entre órbitas permissiveis para 
o átomo de hidrogênio. 

Neste ponto pode parecer que o problema está completamente resolvido. 
Entretanto quando o espectro do hidrogênio é examinado com mais cuidado, 
verifica-se que as linhas previstas pela teoria de Bohr são na realidade compostas 


[7] Introdução à física dos semicondutores 


Figura 1.5 Algumas transições correspondentes às linhas do espectro do hidrogênio 


de várias linhas muito próximas, que constituem à chamada estrutura fina do 
espectro. Esse desdobramento das linhas do espectro, que não é previsto na 
teoria original de Bohr, pode ser explicado a partir dos princípios da Mecânica 
quântica. 


EXERCÍCIOS 


(*1.1. Calcule, usando a Eq. 1.4, o raio da órbita de menor energia do átomo de hidro- 
gênio, o chamado raio de Bohr. 

1,2. Calcule as fregiiências das radiações correspondentes à série de Lyman c à série 
de Balmer. 
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2 Noções básicas de mecânica 
quântica 


Como vimos no capitulo anterior, provara-se que a teoria de Bohr continha 
os fundamentos necessários à solução do problema do espectro dos átomos. 
Conceitualmente, a maior dificuldade dessa teoria estava associada à própria 
idéia de quantização, usada por Planck, Einstein e Bohr, esse último quanti- 
zando inclusive o momento angular. Qual a justificativa para a hipótese de que 
a energia é o momento angular tinham que ser múltiplos de uma certa constante? 
Por que a mecânica clássica resolvia alguns problemas com precisão suficiente, 
enquanto em outros problemas levava a resultados absurdos? Existiria uma 
teoria mais geral, da qual a mecânica clássica seria apenas um caso particular? 

Um grande passo para resolver essas dificuldades foi dado em 1924, quando 
de Broglie (1925, 1926) apresentou sua teoria ondulatória da matéria. De acordo 
com de Broglie, a matéria não é apenas corpuscular, mas qualquer partícula 
também possui propriedades ondulatórias, sendo caracterizada por um compri- 
mento de onda dado por 


dee: (Eq. 2.1) 


onde h é a constante de Planck, m a massa e v a velocidade da partícula. 

Essa hipótese de de Broglie não pode ser demonstrada teoricamente, mas 
Davidson e Germer (1927) mostraram que os elétrons são difratados por uma 
substância cristalina e os comprimentos de onda observados são exatamente 
os previstos pela teoria de de Broglie. 

É interessante notar que a teoria proposta por de Broglie dá um caráter 
todo especial às orbitas estacionárias imaginadas por Bohr. De acordo com as 
Egs. 1.3 e 2.1, o perimetro da circunferência de uma órbita de Bohr é dado por 
nh/2n a ho 

mt me 

Assim a concepção de de Broglie dá às órbitas previstas por Bohr o signi- 
ficado de que a circunferência da órbita é igual a um número inteiro de compri- 
mentos de onda do elétron (onda estacionária). A Fig. 2.1 ilustra esse fato. 

Em 1925, Schroedinger criou a Mecânica ondulatória (Schroedinger, 1926), 
partindo da equação clássica que descreve as ondas estacionárias e usando a 
concepção de de Broglie para definir o comprimento de onda do elétron. Com 
isso chegou a uma equação, atualmente conhecida como equação de Schroe- 
dinger, que, convenientemente interpretada, conduz à solução geral do problema 
das órbitas dos elétrons. 

Na mesma época em que Schroedinger publicou seus trabalhos sobre 
Mecânica ondulatória, Max Born e Werner Heisenberg criaram a chamada 
Mecânica matricial”, que à primeira vista era independente da formulação 


c=2n=27 = nà (Eg. 2.2) 


HA formulação matricial envolve o uso de operadores que geralmente são estudados em 
cursos mais avançados 
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Figura 2.1 Onda estacionária 


Vir 2 


de Schroedinger. Estudos subseqgiientes demonstraram a total equivalência das 
duas teorias, que atualmente são englobadas sob a designação geral de Mecânica 
quântica. Neste capítulo adotaremos a formulação de Schroedinger, já que o 
uso da Mecânica matricial exige uma matemática mais sofisticada. 

Do ponto de vista didático, o maior problema do estudo da Mecânica 
quântica é o caráter postulatório da apresentação. Em outras palavras, a equação 
de Schroedinger tem como única justificativa o fato de que sua aplicação conduz 
a resultados teóricos que concordam com os dados experimentais. 

Outro ponto que contribui para a dificuldade da Mecânica quântica é o 
seu caráter probabilístico. Vamos estudar, por exemplo, uma função que se 
relaciona à probabilidade de encontrarmos um elétron em um determinado 
estudo. Na Mecânica quântica nunca sabemos com certeza se o elétron está 
no estado À ou no estado B; tudo que conhecemos é a probabilidade de que 
esteja em cada um desses estados. 

Schroedinger propôs que o movimento de um elétron submetido a um 
potencial V obedece à seguinte equação: 

h” E Cid ôº h dl), 


dm lgx ? gy + o | Wir D=-Vo, Wir) =— (Eg. 2.3) 


Oy 
onde h = h/2n; x,y e z são as coordenadas do elétron; r é o raio vetor que define 
a posição do elétron; t é o tempo; m é a massa do elétron e V(r, t) é o potencial a 
que está submetido o elétron. 

A função W(r, t) que aparece na Eq. 2.3 é chamada de função de onda do 
elétron. Como a quantidade i = ./—l aparece explicitamente no segundo mem- 
bro da Eq. 2.3, a função de onda é em geral uma função complexa. Mas isso 
equivale a dizer que não podemos atribuir-lhe significado físico algum. O que 
possuí significado físico é o módulo da função W(r.t). Postula-se!” que o qua- 
drado do módulo da função de onda representa a densidade de probabilidade — 


teNão existe qualquer argumento teórico que justifique a interpretação atribuída a W(r,1). 
Por isso, aconselhamos o leitor a aceitar a Equação de Schroedinger da mesina forma como aceitou 
a Lei da Gravitação Universal de Newton 
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probabilidade por unidade de volume e por unidade de tempo — de encontrarmos 
o elétron entre rer+dreentretet+ dt Em outras palavras, se chamarmos 
de p(r.t) a probabilidade de encontrarmos o elétron entre Fr er + dr e entre 
te t+ dt, teremos 


p(r, 1) = [W(r, 1) |  dxdydzdt. (Eg. 2.4) 


Se chamarmos de P(r,t) a densidade de probabilidade, teremos 
P(r. e) = [yr = (rt) prt), (Eq. 2.5) 


onde W*(r.1) é o complexo conjugado de W(r, 1). 

A Eq. 2.3 é a Equação de Schroedinger para o caso geral em que o poten- 
cial V é função das coordenadas espaciais e do tempo. Quando o potencial é 
função apenas das coordenadas espaciais, a equação fica bem mais simples, 

Suponhamos que o elétron está submetido a um potencial V(r) independente 
do tempo. Nesse caso, podemos usar o método de separação de variáveis para 
resolver a equação diferencial 2.3, A idéia é tentar” escrever à função de onda 
w(r, 1) como “produto de duas funções, uma que só depende das coordenadas 
espaciais, W(r), e outra que só depende do tempo, T(1). 


Wlrs 1) = W(P) = TU). (Eq. 2.6) 

Substituindo a Eg. 2.6 na Eq. 2.3, obtemos 

2 2 2 2 
cd a + Es + a ur) TO) = VET) T() = Ba (CD TN, (Eq. 2.7) 

m|ôx” Oy" dz ot 
ou 

hº | 0º a -, a = 

To fã apo 3 + E + =| v|- VP yr) TIO = = ui - y (Eq. 2.8) 


e. finalmente. 


: (Eq. 2.9) 


O primeiro membro da Eq. 2.9 é apenas função das coordenadas espaciais, 
enquanto o segundo membro é' apenas função do tempo. Como não fizemos 
qualquer restrição à Eq. 2.9. ela deve ser válida para todos os pontos do espaço, 
em qualquer instante. Para que isso aconteça, é preciso que os dois membros 
sejam iguais a uma mesma constante. Chamando essa constante de E, devemos 


ter 
hn [ 0º õº õº 
RAE Me de RE TRÊ O Ah Wir) ) 
É dx] o da sige-t AM tia 


Wlr) i T(t) Ot 


'*'Muitas vezes o processo de separação de variáveis não permite a solução da equação. O 
leitor pode estudar o assunto em um livro de equações diferenciais 
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Obtemos, portanto, duas equações independentes: 


h? 02 Fela 02 ; ” e Fé 

RSS lã Ee VE sa Wlr) + Vir) ylr) = Etr, (Eq. 2.11) 
dx” Oy Oz 

oT( 

ot 

A equação original foi desmembrada em duas equações, uma que só de- 
pende das coordenadas espaciais do elétron e outra que só depende do tempo. 
A Eq. 2.11 é chamada de Equação de Schroedinger independente do tempo. 

A interpretação da constante E introduzida na Eg. 2.10 é importante, No 
segundo termo da Eq. 2.11, a função de onda Wy aparece multiplicada pela energia 
potencial V(r). Assim para que a equação seja dimensionamente correta E deve 
ser uma encrgia. Além disso, como E é uma constante, essa energia deve ser 
constante. Ora, de acordo tanto com a Mecânica clássica como com a Mecânica 
quântica, a energia total de um sistema é constante (Princípio da Conservação 
da Energia). Assim, a constante que aparece na solução da Equação de Schroe- 
dinger é interpretada como sendo a energia total do eletron (energia cinética 
mais energia potencial). 

A solução da Eg. 2.12 é 


TO, (Eq. 2.12) 


T(t) = et, (Eq. 2.13) 
Assim, de acordo com a Eg. 2.6, temos 
Wlr. 1) = yr) err. (Eq. 2.14) 
E a densidade de probabilidade é dada por'* 
P(r, 0) = lr, 0) We, 0) = WE) lr). (Eq. 2.15) 


A função W(r) possui uma propriedade importante. Suponhamos, por exem- 
plo, que estamos estudando o movimento de um elétron em torno de um núcleo. 
É evidente que em certo instante o elétron estará em algum ponto do espaço. 
Assim a probabilidade de encontrarmos o elétron em algum ponto do espaço deve 
ser unitária, de acordo com o conceito clássico de probabilidade. Nesse caso 
devemos impor à W(r) a seguinte condição, conhecida como condição de nor- 


malização: K ã 
W(r*-y(r) dxdydz = 1 
todo o espaço 
Resumindo, podemos dizer que o movimento de um elétron submetido a 


um potencial V que só depende das coordenadas espaciais é descrito pelas se- 
guintes equações: 


(Eq. 2.16) 


a2 


2 2 2 — 
' [oa + E + e] Wlr) + [E- vô] Wm) =0; (Eq. 2.17) 


2m 
T() = e”, (Eq. 2.18) 
fee “ylr) dxdyde = 1: (Eq. 219) 
WE. 1) = TE) = Wir) et (Eq. 2.20) 


(WObserve que os fatores e-1E!h, que aparece em W(r, t), e e!!! que aparece em W*(r. 1), se can- 
celam na Eq. 2.15 
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EXERCÍCIOS 


2.1. Calcule o comprimento de onda de um elétron que foi submetido a uma diferença 
de potencial V 

2.2. Por que o modelo de Bohr não explica completamente os resultados experimentais? 

2.3. rr que não podemos atribuir sentido físico algum a y? Que significado é dado 
a hy]*? 

24. Se dia de um átomo de hidrogênio é substituído por um méson pi negativo, 
temos o chamado atomo piônico. Supondo que a teoria de Bohr se aplica ao átomo 
piônico, calcule a frequência da radiação emitida pelo méson pi ao passar do 
primeiro estado excitado para o estado fundamental. (A massa do méson pi é 
igual a 273m,, onde m, é a massa do elétron). 

(82.5, Considere uma partícula de massa m submetida a uma força IF que depende apenas 
do deslocamento da partícula em relação a uma certa posição de equilíbrio. Um 
sistema desse tipo é chamado de oscilador harmônico. Definindo a força aplicada 
à partícula atravês da equação F = —-kx, onde x é o deslocamento da partícula 
em relação à posição de equilíbrio e k é uma constante de proporcionalidade, é 
possivel demonstrar que as funções de onda do sistema são dadas por 


Vl) = Ae RH 0) n=1,200 
onde 


mk 


ex FE * H, é o chamado polinômio de Hermite de ordem n. 


x 


Determine a constante A,. 
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3 Aplicações da mecânica quântica 


Agora que dispomos das ferramentas básicas da Mecânica quântica, po- 
demos aplicá-la a dois exemplos importantes. Em primeiro lugar, estudaremos 
a estrutura eletrônica dos átomos, mostrando quais são os estados possíveis para 
os elétrons de um átomo. No segundo exemplo, vamos examinar as proprie- 
dades do chamado Gás de Fermi, de importantes aplicações para o estudo dos 
materiais semicondutores. 


3.1. ESTRUTURA ELETRÔNICA DOS ÁTOMOS 


Consideremos o caso do átomo de hidrogênio, representado esquemati- 
camente na Fig. 3.1. 


E x | + 
Figura 3.1 Representação esquemática do C es , 
átomo de hidrogênio Eletron 


Próton 


Nosso problema é resolver a Eg. 2,3 para o sistema da Fig. 3.1. Sabemos 
que, para o sistema representado, o potencial é apenas função de r e é dado por 
Tg” 


VOl=am6 


(Eq. 3.1) 


O problema se resume portanto à solução de uma equação diferencial. 
Infelizmente a solução completa é bastante trabalhosa, de modo que nos li- 
mitamos a reproduzir aqui os resultados, deixando os pormenores para um 
livro que trate especificamente de Mecânica quântica (Sherwin, 1959). 

Durante a solução do problema, aparecem três parâmetros: para que haja 
solução, esses parâmetros, chamados de n,l e m, devem guardar uma certa re- 
lação entre si: 

n= 1,23, 
| = iZan=] 
m=0,+1,+2,...,+l 


O número n é o mesmo parâmetro usado empiricamente por Bohr para 
quantizar o momento angular, Os parâmetros | é m só podem ser calculados 
pela Mecânica quântica. 

Aparentemente poderíamos caracterizar o estado de um elétron pelos três 
números n,1 e m. Entretanto um outro fenômeno, que somente a teoria da rela- 
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tividade conseguiu explicar de modo satisfatório, exige um quarto número para 
completa caracterização do estado do elétron. Trata-se do spin do elétron, que 
pode ser imaginado como o movimento de rotação do elétron em torno de si 
mesmo; em outras palavras, um novo tipo de momento angular entra em jogo!*. 

A Mecânica quântica relativística demonstra que o parâmetro associado 
ao spin do elétron (m,) só pode assumir os valores +1/2 e —1/2. 

Em resumo podemos dizer que quatro parâmetros, os chamados numeros 
quânticos, são necessários para descrever completamente o estado de um elétron: 
nm em, . É possível mostrar também que a energia total E do elétron'** é 
função apenas do número quântico n (também chamado de número quântico 
principal) e do número quântico 11***); 


E = função (n, 1). (Eq. 3.2) 


Isso equivale dizer que quando dois elétrons têm o mesmo número quântico 
principal e o mesmo número quântico |, eles têm a mesma energia. Os elétrons 
podem entretanto estar em estados diferentes, bastando para isso que um dos 
outros números quânticos seja diferente. Na Tab. 3.1 vemos todos os estados 
possiveis de um elétron cujo número quântico principal é 2. 


Tabela 3.1 Estados possiveis do elétron de um átomo para n = 2 
n ) m m, 

2 0 0 +1/2 Estados com a mesma energia E 
2 (O) O) —1/2 (n=21=0) 

Ps | — +1/2 

2 | | -1/2 

2 asi 0 +1/2 Estados com a mesma energia E 
2 I O) -1/2 (n=21=1) 

2 1 +1 + 1/2 

2 l +1 —1/2 


De acordo com a Tab. 3.1, existem 8 diferentes estados com n = 2. De 
modo geral, é possivel demonstrar que o número de estados diferentes com o 
mesmo número quântico principal n é igual a 2n?. Um grupo de estados que 
possui o mesmo número quântico principal é chamado de camada. Assim, os 
estados disponiveis do elétron do átomo de hidrogênio podem ser dispostos em 
camadas. Quanto menor o número quântico associado a uma camada, menor 
a energia dos estados correspondentes. O estado em que a energia é mínima 
(no caso do hidrogênio, n = 1) é chamado de estado fundamental. 

Nos átomos dos elementos mais complexos do que o hidrogênio, que pos- 
suem mais de um elétron, os estados disponiveis também podem ser dispostos 
em camadas, cada uma associada a um valor do número quântico principal n. 


*“ Anteriormente tinhamos apenas considerado o momento angular do elétron no seu mo- 
vimento em torno do núcleo 
**'Na ausência de um campo magnético externo 
te**Exceto no caso do hidrogênio, em que a energia só depende de n 
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Cada um desses estados so pode ser ocupado por um elétron, devido ao impor- 
tante principio de exclusão de Pauli. Esse princípio estabelece que num mesmo 
sistema físico dois eletrons não podem ocupar o mesmo estado. Em outras pa- 
lavras, os quatro números quânticos não podem ser iguais para dois elétrons. 
Assim, por exemplo, se dois elétrons no mesmo átomo possuem os mesmos 
valores de n,l e m, o valor de m, pata esses elétrons não pode ser o mesmo; um 
tem que term, = 1/2 e o outro m, = 1/29. 

Em qualquer átomo a ordem de ocupação dos estados é tal que a energia 
total do átomo, isto é, a soma das energias dos estados ocupados, é a menor 
possível. De modo geral, isto significa que as diferentes camadas são ocupadas 
sucessivamente, à medida que aumenta o número de elétrons. Em certos ele- 
mentos, entretanto, como a prata, o índio, o estanho, etc., a interação entre os 
elétrons favorece energeticamente a ocupação de uma camada antes que a ca- 
mada anterior esteja completa. 

Os diferentes estados dos elétrons de um átomo são designados de vários 
modos. A maneira mais conveniente, naturalmente, seria especificar os quatro 
números quânticos. Existem, entretanto, dois sistemas de nomenclatura que 
datam dos primeiros dias da espectroscopia. O primeiro é usado quando estamos 
interessados apenas no número quântico principal. Nesse caso as camadas são 
designadas por letras sucessivas do alfabeto, começando pela letra K, Assim, 
a camada K corresponde a n = 1, a camada La n = 2, etc. O segundo é usado 
para especificar os números quânticos n e |. Nesse sistema de nomenclatura, 
cada grupo de estados é indicado por um número e uma letra. O número é o 
número quântico principal. À letra indica o valor de |. Por motivos históricos, 
os estados correspondentes a | = O são indicados pela letra s (do inglês sharp); 
os estados em que | = 1 são indicados pela letra p (principal): os estados em 
que |! = 2 são indicados pela letra d (diffuse). Os estados em que | = 3 são indi- 
cados pela letra f/ ( fundamental). De | = 4 em diante, os estados são indicados 
em ordem alfabética, começando pela letra g. Assim o grupo de estados 2p inclui 
todos os estados em que n=2 e |= 1; grupo 5d compreende os estados em 
que n=5el=2,cete, Na Tab. 3.2 vemos os estados normalmente ocupados 
dos diferentes elementos. 


3.1.1. Espectros dos elementos 


A distribuição da Tab. 3.2 se aplica aos átomos no estado normal, não- 
-excitado, em que a energia é minima (estado fundamental). Quando um elétron 
em um certo estado recebe energia suficiente, pode passar para um estado de 
maior energia que se encontre desocupado. A energia necessária para causar 
essa transição pode ser fornecida por uma onda eletromagnética ou por um 


(*'As partículas que têm spin + 1/2, +3/2,... são chamadas particulas anti-simetricas (o elétron, 
o próton e o nêutron estão nessa categoria). As particulas que têm spin O, +1, +2,... são chamadas 
particulas simétricas (os fótons estão incluidos nesse grupo). O Princípio de Exclusão de Pauli não 
se aplica apenas aos elétrons, mas a todas as particulas anti-simétricas. Assim, O principio de exclusão 
de Pauli pode ser enunciado do seguinte modo: “Duas particulas idênticas anti-simétricas não podem 
ocupar o mesmo estado”, É interessante ainda observar que as partículas anti-simétricas obedecem 
à chamada estatistica de Fermi-Dirac, (Cap. 7), sendo portanto chamadas de fermions. As particulas 
simétricas obedecem à estatistica de Bose-Einstein, sendo portanto chamadas de hosons. O elétron 
é, consequentemente, um férmion. Essa nomenclatura é muito usada nos livros de Mecânica Quântica 
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Tabela 3.2 Estrutura eletrônica dos elementos 


Número quântico principal n Í 4 
Nomenclatura das camadas K L M N 
0 


B- d, ido sd E o NM, 
3p 3d 4 dp 


Número quântico | 
Notação espectroscópica 1 
Z Elemento 


a 
p 
bo 

a] 
Ls 
a 


1 Hidrogênio (H) 1 

2 Hélio (He) 2 

3 Lítio (Li) 2 1 e 

4 Berilio (Be) 2 2 — 

5 Boro (B) 2 + 

6 Carbono (C) 2 2 2 

7 Nitrogênio (N) 2 2 3 

8 Oxigênio (0) 2 2 

9 Flúor (F) 2 2 5 

10 Neônio (Ne) 2 2. 8 

11 Sódio (Na) 1 o — 

12 Magnésio (Mg) 2. — 

13 Alumínio (Al) 2 q 

14 Silício (Si) Configuração 2 2 

15 Fósforo (P) do neônio Z à 

16 Enxofre (8) 2 4 

17 Cloro (Cl) 2 5 

18 Argônio (A) 2 6 

19 Potássio (K) — 1 
20 Cálcio (Ca) — 2 
21 Escândio (Sc) LN - 
22 Titânio (Ti) E ARIR- 
23 Vanádio (Vv) 3 2 
24 Crômio (Cr) 5 mil 
25 Manganês (Mn) o po 
26 Ferro (Fe) q A 
27 Cobalto (Co) Configuração E RR 
28 Niquel (Ni) do argônio da 
29 Cobre (Cu) 10 1 — 
3% Zinco (Zn) 10 2 — 
31 Gálio (Ga) 38: + 2 1 
32 Germânio (Ge) Ly Du SB 
33 Arsênio (As) JK 4 “a 
34 Selênio (Se) 10 2d 
315 Bromo (Br) lg 3 5 
36 Criptônio (Kr) 10. 2 6 
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Tabela 3.2 (continuação) 


Número quântico principal n 4 5 
Nomenclatura das camadas N [0] 
Número quântico | 2» 3 0 l 
Notação espectroscópica 4d 4f 5s 5p 
Configuração 
Z Elemento das camadas 
internas 
37 Rubídio (Rb) — — I 
38 Estrôncio (Sr) — — 2 
39 Ítrio (Y) Ll— 2 
40 Zircônio (Zr) Configuração 2 — 2 
41 Nióbio (Nb) do criptônio 4 — | 
42 Molibdênio (Mo) 5 — 1 
43 Tecnécio (Tc) 6 — 1 
44 Rutênio (Ru) 7 — 1 
45 Ródio (Rh) 8 — ] 
46 Paládio (Pd) 10 — — 
47 Prata (Ag) a) | ma 
48 Cádmio (Cd) — 2 — 
49 Índio (In) Configuração — 2 | 
50 Estanho (Sn) do paládio — 2 2 
51 Antimônio (Sb) — 2 3 
52 Telúrio (Te) — 2 4 
53 Todo (1) — 2 5 
54 Xenônio (Xe) — 2 6 
55 Césio (Cs) — 
56 Bário (Ba) 
57 Lantânio (La) — 
58 Cério (Ce) 2 
59 Praseodimio (Pr) 3 
60 Neodímio (Nd) - 
61 Promécio (Pm) As camadas Is 5 As camadas 
62 Samário (Sm) a 6 Ssa Sp 
ai a 4d contêm F- 
63 Európio (Eu) 46 elétrons 7 contêm 8 
64 Gadolinio (Gd) 7 elétrons 
65 Térbio (Tb) 8 
66 Disprósio (Dy) 9 
67 Hólmio (Ho) 10 
68 Érbio (Er) o! 
69 Túlio (Tm) 12 
70 Itérbio (Yb) 13 


71 Lutécio (Lu) 14 


oliva 


a 
[A 


SS EI TI DS TS DS dO 3 | DI SI DS SS DI dO ma 
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Tabela 3.2 (continuação) 


Número quântico principal n 5 6 7 
Nomenclatura das camadas (8) P Q 
Número quântico | 2 3 fis «ll 2 0 
Notação espectroscópica 5d 5Sf 6 6 6d 7s 
Configuração 
Z Elemento das camadas 
internas 

72 Háfnio (Hf) 2 — 2 

73 Tântalo (Ta) 3 — 2 

74 Tunestênio (W) As camadas 4 — 2 

75 Rênio (Re) Isaspcontim 5 — 2 

76 Ósmio (Os) 68 elétrons G& — 2 

77 Irídio (1) 9 — 0 

78 Platina (Pt) q es Jd 

79 Ouro (Au) — Le Fe 

80 Mercúrio (Hg) — o 

81 Tálio (TI) — 2 i 

82 Chumbo (Pb) — 2 2 

83 Bismuto (Bi) — 2 3 

84 Polônio (Po) — 2 4 

85 Astatínio (At) As camadas — 2 5 

86 Radônio (Rn) Isa Sd contêm = of 

87 Frâncio (Fr) 78 elétrons — 2 6 — 1 
88 Rádio (Ra) — 2 6 — 2 
89 Actínio (Ac) — 2 6 1 2 
90 Tório (Th) WU) 2 6 í1(2) 2 
91 Protactinio (Pa) M1) 2 6 102) 2 
92 Urânio (U) 3 2 6 i 2 
93 Netúnio (Np) s(494 2 6 +1) 2 
94 Plutônio (Pu) 6(5) 2 6 —L1) 2 
95 Americio (Am) 7 2 [5] — 2 
96 Cúrio (Cm) 7 2 6 1 2 
97 Berquélio (Bk) Bo Donf — 2 
98 Califórnio (Cf) 10 2 (o) E 2 
99 Einsteinio (Es) Io do fi — 2 
100 Férmio (Fm) 2 2 6 — 2 
101 Mendelévio (Md) 13 2 o fiy = 2 
102 Nobélio (No) 14 2 6 — 2 
103 Laurêncio (Lw) 14 2 6 ] 2 


choque com outra partícula. Depois de um certo tempo, o átomo volta espon- 
taneamente ao estado fundamental, emitindo o excesso de energia sob a forma 
de radiação eletromagnética. 

Assim os átomos de um elemento, quando excitados, emitem ondas eletro- 
magnéticas com comprimentos de onda bem característicos, que variam de 
elemento para elemento. O conjunto de radiações emitidas pelos átomos de um 
elemento é chamado de espectro do elemento. 
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3.1.2. Representação esquemática de um átomo 


Vamos concluir nosso estudo da estrutura eletrônica dos átomos apresen- 
tando duas representações gráficas muito usadas na prática. À primeira dessas 
representações está ilustrada na Fig. 3.2 para o caso do silício, cuja estrutura 
eletrônica, de acordo com a Tab. 3.2, é K?LSM*, onde os números representam 
o número de elétrons em cada camada. Nessa representação gráfica, as camadas 
são representadas por círculos e os elétrons em cada camada por pontos. O 
circulo central representa o núcleo. O número indicado no interior do núcleo 
representa a carga elétrica total do núcleo, que em um átomo no estado normal'* 
é igual à soma das cargas de todos os elétrons, 


Figura 3.2 Representação esquemática de um átomo 
de silicio 


O átomo de germânio, cuja estrutura eletrônica é K?LºM'SN*, é repre- 
sentado como na Fig. 3.3. 


Figura 3.3 Representação esquemática de um (32) 
átomo de germânio tu 


Como os fenômenos que vamos estudar neste livro envolvem apenas os 
elétrons da camada externa (os chamados elétrons de valência), podemos usar 
uma representação mais simples, em que apenas os elétrons de valência são 
indicados. Nesse caso, associamos à parte central do átomo uma carga igual 
à diferença entre a carga positiva do núcleo e a carga negativa dos elétrons res- 
tantes. Como tanto o silício como o germânio possuem quatro elétrons de va- 
lência (na camada M e na camada N, respectivamente), essa representação 
simplificada é a mesma para os dois elementos (Fig. 3.4). 


32. GÁS DE FERMI 


Duas razões importantes justificam o estudo do Gás dc Fermi neste livro, 
Em primeiro lugar, o Gás de Fermi é um exemplo simples e ilustrativo da apli- 


(+Não-ionizado 
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Eletron de 


PAi Valência 


Figura 3.4 Representação esquemática simplificada dos 
átomos de silício e germânio 


(O) 


cação da Mecânica quântica a um problema prático. Em segundo lugar, usa- 
remos uma das equações deduzidas nesta seção para o estudo do comporta- 
mento dos elétrons em um material semicondutor. 

O Gás de Fermi é definido como um conjunto de partículas totalmente 
livres, mas sujeitas ao Princípio de Exclusão de Pauli. Em muitos metais, os 
elétrons de condução se comportam como se constituissem um Gás de Fermi. 
Consideremos por exemplo o caso do cobre, cuja estrutura eletrônica é 
1522522pº3s?3pº3d!º4s! (os sobrescritos indicam o número de elétrons em cada 
estado). O elétron 4s está fracamente ligado à rede cristalina, podendo se mo- 
vimentar com liberdade quase completa no interior do material. Quando usamos 
o modelo do Gás de Fermi para estudar as propriedades elétricas do cobre, 
estamos supondo que os elétrons 4s estão inteiramente livres, isto é, que sua 
energia potencial é nula. É claro que na prática os elétrons “livres” do cobre 
interagem com a rede cristalina e com outros elétrons “livres”, mas essa inte- 
ração é tão fraca que pode ser desprezada em primeira aproximação. 


3.2.1. Gás de Fermi em uma dimensão 


Vamos estudar inicialmente o comportamento do Gás de Fermi em uma 
dimensão. Suponhamos um elétron de massa m confinado a uma linha de com- 
primento L por barreiras infinitas de potencial nas extremidades da linha (Fig. 
3.5). 


sds BARREIRAS o 
INFINITAS DE POTENCIAL 


Figura 3.5 Gás de Fermi em uma dimensão 


V=0 
X=0 X=L 
Nosso problema é resolver a Equação de Schroedinger (Eq. 2.17) para o 
caso em que V= 0 no interior da linha”, Nesse caso, temos 


hº jlx) o E 
mm aa * EV) =0. (Eq. 3.3) 


Como a função de onda J(x) é uma função continua e como W(x) deve ser 
zero para x = 0 c para x 2 L. já que o elétron não pode atravessar uma bar- 
reira infinita de potencial, a solução da Eq. 3.3 deve satisfazer às condições 


*A condição V = O significa que o elétron está livre 
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de contorno: 
WO) = WL) = 0. (Eq. 3.4) 


Como estamos interessados em obter todos os estados disponiveis para os 
elétrons, vamos caracterizar cada estado por um índice n colocado em iWy(x) e E, 
isto é, vamos escrever 


h2 d2y lx) 
— — + + E = (); Eq. 3.5 
do dd nn) (Eq. 3.5) 
WO) = ylL) = 0. (Eq. 3.6) 
O leitor pode facilmente resolver essa equação diferencial, obtendo 
b=AsenDs n=123.:,5 (Eq. 3.7) 
h2 (nn? 
E, = Ze) (Eq. 3.8) 
Yamos proceder inversamente e simplesmente mostrar que 
nTx 
y, = Asen EE : 
é uma solução da equação. 
De fato, 
Mp RA (o 
dx L L 
dy, nm nax 
dx? E? L 
Portanto 


2 2 
E E (n i) | A sen E =-E, W x). 


Assim para um elétron confinado a uma linha de comprimento L, os estados 
disponíveis são caracterizados por uma função de onda da forma 


wy, = A sen aaa (Eq. 3.9) 
e por uma energia da forma 
hº 2 
E, = mt) (Eq. 3.10) 


Para calcular a constante A, usamos a condição de normalização: 


L 
) [Wnloo)]? dx = 1. (Eq. 3.11) 
O) 


22 Introdução à física dos semicondutores 


No nosso caso, temos 


A? ds =, (Eq. 3.12) 
[1H] 
2 
PP dio le fe pec (Eq. 3.13) 
2 hTx L L 
en  * E) 


2 
A= JE: (Eg. 3.14) 


Em conseqjiiência, as expressões finais para as funções de onda e para as energias 
permitidas são as seguintes: 


p= |osenTO n=1,2,3...: (Eq. 3.15) 


2 2 
E = 2a) ER UÉ, (Eq. 3.16) 


De acordo com a Eq. 3.16, a energia de um elétron confinado à linha não 
pode variar continuamente; os únicos valores permitidos para a energia do 
elétron são da forma 


E, = na, (Eg. 3.17) 


A ht (5): 
2miL 

Suponhamos agora que estamos interessados em acomodar N elétrons na 
linha. De acordo com o Princípio de Exclusão de Pauli, o mesmo estado não 
pode ser ocupado por mais de um elétron. Para cada valor de n na Eg. 3.15, 
existem dois valores possíveis para O spin do elétron, +1/2 e -1/2, e portanto 
dois estados disponíveis. À temperatura do zero absoluto (0 K), o sistema se 
encontra no estado fundamental (estado de mínima energia) e portanto os estados 
disponíveis são preenchidos em ordem crescente de energia. Assim, por exemplo, 
se existem 13 elétrons na linha, a distribuição desses elétrons entre os diferentes 
estados a OK será a seguinte: 


onde: 


1 2 elétrons, 

2 2 elétrons, 

3 2 elétrons, 
=4 2 elétrons, 

S 2 elétrons, 

6 2 elétrons, 
n= 1 elétron, (nível incompleto). 
n>7 O clétrons, (níveis vazios). 


A energia mais elevada para a qual ainda existem estados ocupados a OK 
(E,, no exemplo anterior) é chamada de energia de Fermi, 
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Repare que como a energia é uma função quadrática do número quântico n 
(Eg. 3.17), o número de estados disponiveis dentro de uma dada faixa de ener- 
gias AE diminui à medida que n aumenta (e portanto à medida que E aumenta). 
Para analisarmos quantitativamente a variação do número de estados com a 
energia, definimos a chamada densidade de estados D(E) como o número de 
estados disponiveis por unidade de energia. 

Vamos agora calcular a densidade de estados do Gás de Fermi em uma 
dimensão. Nosso ponto de partida pode ser a Eq. 3.16. Considerando a energia 
E, como uma função contínua de nº”, podemos diferenciar a Eg. 3.16, obtendo 


h? (ny 
dE, = Sal) 2n dn, (Eg. 3.18) 
dn Lya 


Como existem dois estados disponiveis para cada valor de n, a densidade 
de estados é dada por 


dn E? 1 4L (my? 
= 7 aum e — == tm , 9 
D(E) =2 JE, E 5) Eae" Do 2) (Eq. 3.20) 
ou, finalmente, 
la Es 
= —— 2 p= 1/2 
DE) ==, V'2m ET u2 (Eq. 3.21) 


A Fig. 3.6 é um gráfico da variação da densidade de estados D(E) com a 
energia para o Gás de Fermi em uma dimensão. À temperatura de OK, todos 
os estados com energias entre O e E, (energia de Fermi) estão ocupados e todos os 
estados com energias maiores que E, estão vazios. 


D(E) 


Figura 3.6 Variação da densidade de estados 
com a energia para o Gás de Fermi em uma 
dimensão 


1 
Er E 


tsEmbora não seja estritamente válida, já que, como vimos, o número n só pode assumir va- 
lores inteiros, essa aproximação é usada em muitos problemas de Mecânica quântica e em geral 
conduz a resultados corretos 
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3.2.2. Gás de Fermi em três dimensões 


Embora matematicamente mais complexo, o problema do Gás de Fermi 
em três dimensões não passa de uma extensão do problema unidimensional. 
Suponhamos, por exemplo, que os elétrons estão confinados a um cubo de 
aresta a (Fig. 3.7). 


Figura 3.7 Gás de Fermi em três dimensões 


Nosso problema é resolver a equação de Schroedinger em três dimensões 
h?[ 0? õ? 
>m| 9x2 | avi 
2m| Ox Oy 


2 
+ | Wlx, y, 2) + Ex, y,7) = 0, (Eg. 3.22) 


sujeita às seguintes condições de contorno: 


WO, y, 2) = Wa, y, 2) = 0; 
W(x, 0, 2) = W(x,a, 2) = 0; (Eq. 3.23) 
WO, y, 0) = lx, y, a) = 0. 
As Egs. 3.23 são uma consegiiência do fato de que a função de onda deve 
se anular nas faces do cubo. 


Por analogia com a solução unidimensional, Eqgs. 3.7 e 3.8, a solução da 
Eg. 3.22 é dada por 


naAx o nn) naz 


Wlx, y.Z) = A sen sen = Sen E é (Eq. 3.24) 
E, = Keri + nº + n2] (Eq. 3.25) 
n ?ma? x y E) da 


onde 


No Non, =1,2,3,.0 


Fica a cargo do leitor verificar que as funções de onda (Eq. 3.24) satisfazem 
à equação de Schroedinger (Eq. 3.22). 
Para calcularmos a constante 4 usamos a condição de normalização seguinte: 


| | | [Wlx, y, 2)[P dxdydz = 1. (Eg. 3.26) 
o do o 
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Substituindo a Eq. 3.24 na Eq. 3.26, temos 


A? | sen2lz2X gr | sent gy | sento =1. (Eq. 327) 
0 ia 0 q 0 q 


Como cada uma das integrais é igual a a/2, temos 


| 8 
A = = (Eq. 3.28) 
e consequentemente 


8 nAx n,Zy n.Az 
o x pH s02, 
Wlx, YZ) = I qi en “= (Eq. 3.29) 


Assim o estado de um elétron confinado a uma região tridimensional é carac- 
terizado por quatro números quânticos: três números quânticos espaciais, n,. 
N,s N.. que podem assumir os valores 1,2,3,... e um número quântico de spin, 
m,. que pode assumir os valores +1/2 e —1/2. 

Vamos agora calcular a densidade de estados do Gás de Fermi em três 
dimensões. 

Nosso ponto de partida será a equação que nos dá as energias permitidas 
para o elétron: 


hº nº 
E=— (n2 +nº + n2), (Eq. 3.30) 
2ma 


Se representarmos os valores de n,, n, e n, em um sistema de eixos car- 
tesianos, o lugar geométrico dos pontos de mesma energia E será dado por 


5 ma?  2ma e 
no +n;+ni= PER E [e pr] : (Eq. 3.31) 
o que corresponde à superficie de uma esfera de raio 
RE 
PAR denad EM2, (Eq. 3.32) 


hr 


conforme vemos na Fig. 3.8. 


Figura 3.8 Superficie de energia constante de 
um Gás de Fermi em três dimensões 
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Como os números quânticos n,, n, e n, são positivos por definição, todos 
os estados de energia menor que E estão compreendidos em um volume igual 
a 1/8 do volume de uma esfera de centro na origem e raio R. Nesse caso, o nú- 
mero de estados com energia menor que E é dado por 


1 4 3 [2m |? 
ta= 25 aR= [| ER, (Eg. 3.33) 


onde o fator 2 se deve à existência de duas orientações possíveis para o spin 
do elétron. 
Procedendo como para o caso unidimensional, diferenciamos a Eqg. 3.33, 
obtendo : EA 
ce Di 1/2 : 
D(E) = JE — SÊ [| Elté, (Eq. 3.34) 
A Eg. 3.34 nos dá a densidade dc estados de um Gás de Fermi em três di- 
mensões. Podemos definir também a densidade de estados por unidade de volume 
como sendo a relação entre a densidade de estados e o volume no qual os elé- 
trons estão confinados. Como supusemos que os elétrons estavam confinados 
a um cubo de lado a (volume a”), a densidade de estados por unidade de volume 
do Gás de Fermi em três dimensões é dada por 


h2 


Como veremos no Cap. 8, a Eq. 3.35, embora deduzida para o Gás de Fermi, 
pode ser reinterpretada de modo a representar a densidade de estados por uni- 
dade de volume dos elétrons de condução de um semicondutor. 


3/2 
DE) = 5 DIE) = = EM2, (Eq. 3.35) 


EXERCÍCIOS 


3.1, Faça uma lista dos estados possíveis de um elétron de um átomo cujo número quân- 
tico principal é 3. 
3.2. Mostre que a energia de Fermi do Gás de Fermi em uma dimensão é dada por 


RR ÓLA 
“rm QL)” 


onde N é o número total de elerrons. 

3.3. Neste capítulo, para obtermos as densidades de estados nos casos unidimensional 
e tridimensional, consideramos a energia como uma função de n. Analise cuida- 
samente a validade dessa aproximação. 

3.4. Estude o comportamento de vm Gás de Fermi confinado a uma semi-esfera de centro 
na origem dos eixos cartesianos, raio R ce face plana apoiada no plano xy(z > 0). 
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4 Noções de cristalografia 


4,1. SUBSTÂNCIAS CRISTALINAS 


Chamamos de substância cristalina qualquer substância cujos átomos estão 
dispostos regularmente no espaço (Fig. 4.1). 


Figura 4.] Estrutura de uma substância cristalina 


Um modo interessante de caracterizar um cristal como o da Fig. 4.1 é de- 
finir uma estrutura puramente geométrica, a cujos pontos são associados um 
ou mais átomos para reproduzir o cristal. Esse método de representar as 
substâncias cristalinas está ilustrado na Fig. 4.2. 


Pontos 


(a) tc) 
Figura 4.2 A estrutura cristalina que aparece em (c) é formada associando-se um átomo 
(b) a cada ponto de uma estrutura puramente geométrica (a) 


A estrutura puramente geométrica é chamada de rede cristalina da subs- 
tância considerada”. O conjunto de átomos que é associado a cada ponto da 
rede cristalina (no caso da Fig. 4.2, apenas um átomo) é chamado de base. 
Qualquer cristal pode ser definido através de uma rede cristalina associada a 
uma base apropriada. 

A Fig. 4.3 mostra a rede cristalina c a basc correspondente para o caso 
do cloreto de sódio (NaCl). Repare que neste exemplo a base é constituída por 
dois átomos diferentes, Na e CI. 


“Crystal lartice, em inglês 
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tc) 
Figura 4.3 Cristal de cloreto de sódio. (a) rede cristalina; (b) base; (c) cristal completo 


É importante observar que o mesmo cristal pode ser representado por mais 
de uma rede cristalina. Nesse caso, naturalmente, as bases associadas a cada 
rede serão diferentes. 


42. CÉLULA UNITÁRIA E SISTEMAS CRISTALINOS 


Um fato importante a respeito das redes cristalinas é que qualquer rede 
cristalina tridimensional pode ser gerada pela translação de um paralelepípedo 
ao longo de três eixos mutuamente perpendiculares. Assim, por exemplo, a rede 
cristalina da Fig. 4.2 pode ser gerada pela translação de um cubo (Fig. 4.4). 


Figura 4.4 Geração de uma rede cris- 
talina pela translação de um cubo 


O paralelepípedo que gera a rede cristalina é chamado de célula unitária 
da rede. A mesma rede cristalina pode ser gerada por paralelepípedos de di- 
ferentes tamanhos. O menor paralelepípedo capaz de gerar a rede cristalina é 
chamado de celula primitiva da rede. 

De acordo com as relações entre os lados e os ângulos da célula unitária, 
as redes cristalinas podem ser divididas em sete sistemas cristalinos. Cada um 
desses sistemas por sua vez compreende um ou mais tipos de rede cristalina, 
de acordo com as simetrias apresentadas pela rede. Existem ao todo 14 tipos 
possíveis de rede cristalina tridimensional. 

As Fipgs. 4.5 a 4.11 mostram as células unitárias (que não são necessariamente 
células primitivas) dos 14 tipos de redes cristalinas. 

Convém lembrar que, como a rede cristalina, a célula unitária é uma estru- 
tura puramente geométrica, a cujos pontos deve ser associada uma base ade- 
quada para gerar o cristal. É evidente que podemos definir qualquer cristal por 
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al 
ls ET ge = 
[b) (c) 


a=b=c 

X=p=Y=90" 
Figura 4.5 Sistema cúbico; (a) cúbico simples; (b) cúbico de corpo centrado; (c) cúbico 
de faces centradas 


Figura 4.6 Sistema tetragonal; (a) tetragonal sim- 
ples; (b) tetragonal de corpo centrado 


(a) (b) 
a=bzct 
A =B= $ =90º 
7 dE DS 
- AM TE 
(a) Cb) (c) (d) 
azbacr 
d=p=4=90º 


Figura 4.7 Sistema ortorrômbico; (a) ortorrômbico simples; (b) ortorrômbico de bases 
centradas; (c) ortorrômbico de corpo centrado; (d) ortorrômbico de faces centradas 


SH os 


agzbzc azbzáco 
a=9=90"%p dzpey 
Figura 4.8 Sistema monoclinico; (a) monoclinico Figura 4.9 Sistema triclinico 


simples; (b) monoclinico de bases centradas 


uma célula unitária com os átomos devidamente posicionados, já que uma 
estrutura desse tipo é capaz de gerar todo o cristal por translação. Por que, 
então, nos darmos ao trabalho de definir separadamente uma rede cristalina e 
uma base? A razão é que em um tratamento matemático mais profundo dos 
cristais, a definição do cristal pela rede e pela base torna as equações mais simples. 
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À a=bzc 
Q=p=e «= PR = 90º 
A=P=ó <120º, 490º g=120º 
Figura 4.10 Sistema trigonal Figura 4.11 Sistema hexagonal 


Por outro lado, como estamos interessados apenas nas aplicações da cris- 
talografia ao estudo das propriedades dos materiais semicondutores, vamos nos 
ater à representação das células unitárias já com os átomos posicionados, o que 
permite uma visualização mais simples da estrutura cristalina. A Fig. 4.12 ilustra 
o que acabamos de explicar. 


tb) 


(a) 
Figura 4.12 (a) Célula unitária cúbica de corpo centrado; (b) base a ser associada à 
rede cristalina; no caso, átomo de cromo; (c) célula unitária do cromo com os átomos 
posicionados 


A partir deste momento vamos apresentar as estruturas dos cristais na 
forma da Fig. 4.12(c), deixando a cargo do leitor a identificação da rede cristalina 
e da base correspondente. Outros elementos que possuem a mesma estrutura 
cristalina que o cromo (cúbica de corpo centrado) são o lítio, o bário, o nióbio, 
o césio e o tungstênio. 


Cs Cs 


cs 


Cs Cs 
Figura 4.14 Estrutura cristalina do 
cloreto de césio 
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A Fig. 4.13 mostra a estrutura cristalina do aluminio (cúbica de faces cen- 
tradas). Outros elementos com a mesma estrutura são o cobre, o ouro, o chumbo, 
o niquel c a prata. 

A Fig. 4.14 mostra um exemplo de cristal composto, o cloreto de césio. 
A rede cristalina do cloreto de césio é cúbica simples; a base consiste em dois 
átomos, um de cloro e um de césio. 

Finalmente apresentamos na Fig. 4.15 a célula unitária do diamante, que 
nada mais é que o carbono na forma cristalina”, 

Aproveitamos esse exemplo para introduzir o parâmetro a, a chamada 
constante da rede; como vemos na Fig. 4.15, a constante de rede é igual ao lado 
da célula unitária. 

É interessante observar que a estrutura da Fig. 4.15 é exatamente a mesma 
do germânio e do silício, duas substâncias de vital importância para a tecnologia 
dos dispositivos semicondutores; apenas a constante da rede é diferente. 


Figura 4.15 Estrutura cristalina do diamante Figura 4.16 Ligações tetraédricas 


A Fig. 4.16 mostra a parte inferior esquerda da Fig. 4.15, onde estão loca- 
lizados os átomos 4,B,C,D e E. 


Os átomos B, C,D c E formam um tetraedro regular, com o átomo A exa- 
tamente no centro. Assim, o átomo A está à mesma distância dos quatro átomos 
mais próximos, B,C,D e E, Como veremos mais tarde, entre o átomo A e os 
átomos B,C,D e E existe uma ligação especial, chamada ligação covalente, que 
mantém os átomos firmemente presos à rede cristalina. Essa ligação, que no 
caso de cristais como o diamante, o germânio e o silício é também chamada 
de ligação tetraedrica (porque os átomos vizinhos de qualquer átomo consi- 
derado formam um tetraedro) tem uma profunda influência sobre as proprie- 
dades dos materiais semicondutores. 


(*Chamamos a atenção do leitor para a diferença extraordinária entre o carbono na sua forma 
cristalina (diamante) e amorfa (grafita) 
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43. CÁLCULO DE NÚMERO DE ÁTOMOS DE UM CRISTAL 
POR UNIDADE DE VOLUME 


Conhecida a célula unitária de um cristal e a constante de rede, o número 
de átomos do cristal por unidade de volume pode ser calculado facilmente. 
Consideremos, por exemplo, uma substância cristalina cúbica de corpo cen- 
trado (Fig, 4.17). 


As 


- ÀS 


Figura 4.17 Exemplo do cálculo do número de áto- 
A mos por unidade de volume em um cristal 


As 
E 

Como vimos, a célula unitária deve ser transladada ao longo de três di- 
mensões mutuamente perpendiculares para gerar o cristal no espaço. Assim um 
átomo como A,, localizado no vértice do cubo, pertence simultaneamente a 
8 células adjacentes. Para efeito de cálculo, apenas 1/8 desse átomo pertence 
à célula unitária da Fig. 4.17. Como existem 8 vértices, a contribuição dos átomos 
dos vértices para uma célula unitária é igual a 8 x 1/8 = 1 átomo. 

Por outro lado, o átomo As, localizado no centro do cubo, pertence exclu- 
sivamente à célula da ilustração. Assim, o número total de átomos pertencentes 


à célula é dado por: N = 1 (contribuição dos átomos dos vértices) +1 (con- 
tribuição do átomo central) = 2 átomos. 


Como o volume da célula é a?, onde a é a constante de rede, o número de 
átomos por unidade de volume de um cristal como o da Fig. 4.17 é igual a 2/a?. 


Vamos agora calcular o número de átomos por unidade de volume para 
os materiais que nos interessam mais de perto, o germânio e o silício, cuja estru- 
tura cristalina é a representada na Fig. 4.15. 


No caso da estrutura da Fig. 4.15, cada átomo situado no vértice do cubo 
pertence simultaneamente a 8 células, exatamente como o exemplo anterior. 
Assim a contribuição total dos átomos dos vértices é igual a 1 átomo. 


Os átomos situados nos centros das faces do cubo pertencem a duas células 
vizinhas: como existem 6 faces, a contribuição total para uma célula é igual a 
6 x 1/2 = 3 átomos. 


Os 4 átomos situados no interior da célula pertencem exclusivamente à 
célula considerada. 
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Podemos sumarizar o nosso cálculo como segue. 


Número de átomos nos Contribuição de cada átomo Contribuição total para 
vértices da célula = 8 dos vértices para a célula a célula: 
=1/8 8x 1/8=1 
Número de átomos nos Contribuição de cada átomo 6x 1/2=3 
centros das faces da dos centros das faces para a 
célula = 6 célula = 1/2 
Número de átomos no Contribuição de cada átomo 4x1/1=4 
interior da célula = 4 interior para a célula = 1/1 


Número total de átomos 
da célula 8 


Como o volume da célula é a?, o número por unidade de volume será 
número de átomos/unidade de volume = 8/uº 
Os valores numéricos correspondentes são os seguintes: 


Constante 


e rede (nm)* Número de átomos/cm* 
t 


Substâncias 


diamante 0,356 1,77 x 102? 
germânio 0,566 4,42 x 102? 
silício 0,543 5,00 x 107? 


Esse processo para o cálculo do número de átomos por unidade de volume 
de um cristal utilizando-se a célula unitária e calculando-se a contribuição de 
cada átomo para a célula pode ser aplicado a qualquer cristal. 


4.4. POSIÇÃO E ORIENTAÇÃO DOS PLANOS DE UM CRISTAL 


Os átomos de um cristal definem uma série de planos. A Fig. 4.18, por 
exemplo, mostra três dos planos definidos pelos átomos de um cristal cúbico 
de faces centradas. 

O conhecimento da posição e orientação desses planos em uma dada amostra 
pode ser extremamente importante. Muitos cristais, por exemplo, possuem os 
chamados “planos de clivagem”, que apresentam a seguinte propriedade: quando 
golpeamos a amostra na direção de um plano de clivagem, a fratura ocorre 
exatamente nesse plano. Quando a direção do golpe não coincide com nenhum 
dos planos de clivagem, a fratura é irregular e a amostra pode ser estilhaçada. 
A lapidação das pedras preciosas é feita utilizando-se esse princípio. 

A maneira mais imediata de definir um plano de um cristal seria indicar 
as coordenadas dos traços desse plano nos três eixos de referência. Assim, por 
exemplo, o plano P da Fig. 4.19 seria definido pelos pontos: A(a, 0, 0). B(O, a, 0) 
e €(0,0,a), onde a é a constante de rede. 

Em cristalografia, entretanto, prefere-sc definir os planos através dos cha- 
mados índices de Miller, obtidos pelo seguinte processo (Fig. 4.20): 


1 nm (1 nanometro) = 107º m 
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Figura 4.18 Planos de um cristal 


kigura 4.19 Definição de um plano 
cristalino por seus traços 


a) Calculam-se as distâncias da origem aos pontos 4,B e € de interseção 
do plano com os cixos de referência, tomando como unidade a constante de 
rede. No caso da Fig. 4.20, 


OA =1, OB = 2, OC=2. 
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Figura 4.20 Exemplo de cálculo dos indices de 
Miller 


3a 
X 


b) Tomam-se os reciprocos dos valores encontrados: 
1/04 = 1, W/0B= 12, 1/00 = 1/2. 
c) Reduz-se os resultados a números inteiros que guardem entre si a mesma 
relação que os valores do item (b). 
No caso da Fig. 4.20, multiplicamos todos os números por 2, obtendo: 
2x1=2, 2x 12=1, 2x1/2=1. 
Os números obtidos são chamados de indices de Miller e representados 


entre parênteses: (2 1 1). É evidente que os índices de Miller não definem apenas 
um plano, mas toda uma família de planos paralelos. 


Figura 4.21 Plano com traço no semi- 
-eixo negativo 
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Quando o plano intercepta um dos eixos no semi-eixo negativo, atribuimos 
um sinal negativo à coordenada correspondente. A Fig. 4.21 ilustra esse caso. 
O cálculo seria o seguinte: 


ÕÃ = 2, 0B =, OC =3; 
1/04 = —1/2 1/0B= 1/2, 1/0C=1/3; 
6x(-1/))=-3, 6x12=3, 6x13=2. 
Índices de Miller: (332). 
Repare que o sinal negativo é colocado sobre o número. 
A Fig. 4.22 mostra alguns planos de um cristal cúbico simples e os índices 


de Miller associados a esses planos. Repare que quando um dos indices de Miller 
é zero, isto significa que o plano é paralelo ao eixo correspondente. 


ÍNDICES DE MILLER (j00) ÍNDICES DE MILLER (020) 


ÍNDICES DE MILLER (200) 


Figura 4.22 Alguns planos de um cristal cúbico: (a) plano (1 00), (b) plano (0 2 0), 
(c) plano (1 1 1), (d) plano (2 0 0) 


É fácil demonstrar que os índices de Miller de um plano são proporcionais 
aos co-senos diretores da perpendicular traçada da origem ao plano considerado; 
na Fig. 4.23 estão indicados os parâmetros envolvidos. 
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OP-= perpendicular ao plano P 


Figura 4.23 Relação entre os indices de Miller e os 
co-senos diretores da perpendicular baixada da ori- 
gem ao plano considerado 


Os co-senos diretores da reta definida por OP são dados por 
cosa =d/a, cosB=d/b, cosy = de, 
onde d é a distância da origem ao plano considerado. 


Assim os co-senos diretores são proporcionais a 1/a, 1/b, 1/c, do mesmo 
modo que os índices de Miller. 


4.5. DIREÇÕES EM UM CRISTAL 


Convencionalmente, os índices que caracterizam uma direção em um cristal 
são definidos pelos três menores números inteiros que guardam entre si a mesma 
relação que as componentes de um vetor de módulo arbitrário, paralelo à di- 
reção considerada. 

A Fig. 4.24 mostra um exemplo. A direção dada é a direção «d. Um vetor 
v, traçado paralelamente a d, possui as seguintes componentes: 


m=1, m=1, vw=7/2 
Multiplicando por 2 as componentes de v, obtemos 
2%X1=2 2Z2x1l=2, 2x7)W)2=7. 
Direção d: [227]. 
Repare que os índices que definem a direção são colocados entre colchetes. 


Pode-se demonstrar (Prob. 4.5) que em um cristal cúbico o plano (hk) 
é perpendicular à direção [hk/]. Os exemplos da Fig. 4.25 ilustram esse fato. 


4.6. IMPORTÂNCIA DOS ÍNDICES DE MILLER 


Provavelmente o leitor ainda não está convencido da utilidade dos indices 
de Miller. Por que indicar um plano cristalino pelos índices de Miller, quando 
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Figura 4.24 Direção em um cristal 


PLANO P: (100) 


= PLAN 4 
DIREÇÃO DO EIXO DOS X: [100] ro a DO 


DIREÇÃO d: [111411 


Figura 4.25 Exemplos de planos é 
direções com os mesmos indices 


seria muito mais simples defini-lo pelas coordenadas dos seus traços? Na ver- 
dade, a importância dos indices de Miller está intimamente relacionada ao 
conceito de rede reciproca, que infelizmente foge ao escopo deste livro. Mesmo 
assim, fizemos questão de definir os indices de Miller para que o leitor possa 
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compreender como são indicados os planos de um cristal nos artigos científicos. 
Por exemplo, se o leitor encontrar a expressão “um cristal de germânio, com 
o plano (100) paralelo à direção do campo magnético aplicado. ..”, saberá 
imediatamente que o campo magnético é paralelo a uma das faces do cubo da 
célula unitária do cristal de germânio. 


4.7. DETERMINAÇÃO DA ESTRUTURA DE UM CRISTAL 


A tecnologia atual ainda não permite, exceto em raros casos, a observação 
direta da estrutura de um cristal. Como é possível, então. determinar a posição 
dos átomos no interior do cristal, isto é, as dimensões e a composição da célula 
unitária? 

A solução para esse problema consiste em fazer incidir um feixe de radiação 
(fótons, nêutrons, elétrons) sobre o cristal. Quando o comprimento de onda da 
radiação é da mesma ordem de grandeza que a constante de rede do cristal, o 
feixe é difratado pelos átomos do cristal. A análise do feixe difratado permite 
calcular matematicamente as posições dos átomos do cristal. 

Vejamos agora como a análise do feixe difratado nos permite determinar 
a estrutura de um cristal. 

Vamos supor que façamos incidir um feixe de raios X em uma amostra 
cristalina (Fig. 4.26). 


FEIXE INCIDENTE FEIXE/DIFRATADO 


Figura 4.26 Feixe de raios X 
incidindo em um cristal 


A Fig. 4.27 é uma vista de perfil do cristal da Fig. 4.26, que mostra como 
os raios X são difratados pelo conjunto de planos paralelos P,, P,,...do cristal. 

Como vemos na Fig. 4.27, os raios são refletidos pelos átomos localizados 
nos planos P,, P,... e em seguida são registrados em uma placa fotográfica. 
Chamando de d a distância entre os planos e de 0 o ângulo entre o feixe incidente 
e o conjunto de planos paralelos, a distância entre os pontos de incidência no 
filme de dois raios refletidos por planos vizinhos. como P, e P, , é igual a (d cos 8). 
Ora, como q é extremamente pequeno, o ponto de incidência no filme dos raios 
refletidos R,, R,, R,.... é praticamente o mesmo. Em consegiiência a imagem 
produzida no filme pelos raios R,. R,, R,, etc. será um ponto. 
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FILME 


Figura 4.27 Vista de perfil do 
cristal da figura anterior, mos- 
trando o conjunto de planos pa- 
ralelos P,, P;,... 


Para que os raios refletidos produzam uma imagem, entretanto, é preciso 
que estejam em fase, isto é que a diferença de percurso entre eles seja igual a 
um múltiplo inteiro do comprimento de onda da radiação utilizada. De acordo 
com a Fig. 4.27, a diferença de percurso entre os raios refletidos por dois planos 
vizinhos é dada por 


e= 40 + OB = 2dsen 0. (Eq. 4.1) 

Assim à condição para que os planos paralelos P,. P,.... produzam uma 
imagem na placa fotográfica é expressa pela seguinte equação: 

2d senO = nÃ, (Eg. 4.2) 


onde d é a distância entre os planos, O é o ângulo de incidência da radiação, n 
é qualquer número inteiro e À é o comprimento de onda da radiação. A Eg. 4.2 é 
chamada de Lei de Bragg em homenagem ao físico inglês W.L. Bragg, que a 
deduziu em 1913 (Bragg, 1913). 

De acordo com a Lei de Bragg, um conjunto de planos paralelos do cristal 
só produz uma imagem na placa fotográfica para certos ângulos O da radiação 
incidente, que satisfaçam à Eq. 4.2. Por outro lado, como vimos nos parágrafos 
anteriores, um cristal possui muitas famílias de planos paralelos. Assim, variando 
o ângulo de incidência, podemos obter as imagens correspondentes a várias 
famílias de planos. A análise dessas imagens permite deduzir a estrutura do 
cristal e sua orientação em relação ao feixe incidente. 

Na verdade a Lei de Bragg exprime apenas uma condição necessária para 
que o feixe difratado produza uma imagem. Isso porque sua dedução se baseia 
apenas na periodicidade da rede cristalina, sem levar em conta os átomos que 
estão associados a cada ponto da rede. Assim, a Lei de Bragg pode ser satisfeita 
sem que apareça ponto algum na placa fotográfica. 

A influência da base do cristal pode ser incorporada ao estudo do fenômeno 
de difração através da introdução de dois fatores. O primeiro, chamado fator de 
estrutura, representa a influência puramente geométrica dos átomos da base e 
portanto não depende do tipo do átomo considerado. O segundo, o chamado 
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fator de forma, é uma medida do poder de espalhamento dos átomos que cons- 
tituem a base. Para uma discussão elementar desses dois fatores, consulte por 
exemplo Kittel (1967). Um estudo mais completo pode ser encontrado em 
Azaroff (1968). 


EXERCÍCIOS 


(814.1. O ouro tem uma massa específica de 19,32 g/cm*, um peso atômico de 197,0 e 
uma estrutura cúbica de faces centradas. (a) Calcule o número de átomos de ouro 
por centimetro cúbico; (b) Calcule a constante de rede. 

4.2. O sistema tetragonal de faces centradas é equivalente a um dos 14 tipos de rede 
cristalina ilustrados nas Figs. 4.5 a 4.11. Qual é esse tipo? 

4.3. Calcule a distância entre dois átomos vizinhos em um cristal de silício. 

(84.4. Calcule o valor do ângulo existente entre as ligações tetraédricas de um cristal 

de Si. Esse valor difere para o caso do Ge e do diamante? 

4.5. Prove que em um cristal cúbico o plano definido pelos índices (hkl) é sempre per- 
pendicular à direção [hkl]. 

4.6. Mostre as posições atômicas do Ga e do As nos planos (1 10)e (100) de um cristal 
de GaAs. 

tx14.7, Quais os planos que contêm os três tipos de direções de baixo índice, [100], 
[110]e [111]? 

4.8. Em uma experiência de difração, quais as condições que devem ser satisfeitas 
para que se tenha uma imagem na chapa fotográfica, correspondente a um con- 
junto de planos paralelos? 

4.9. Lembrando que a luz visivel está na faixa de comprimentos de onda de 400 a 
800 nm, explique por que não se pode usar luz visivel para estudar a rede cris- 
talina de um cristal por difração. 
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5 Aplicação da mecânica quântica 
a um cristal 


Neste capítulo vamos aplicar as idéias desenvolvidas nos capitulos ante- 
riores ao estudo dos estados de energia dos elétrons em um cristal. O problema, 
em sua forma geral, é extremamente complexo, pois os elétrons, além de serem 
atraídos pelos núcleos dos átomos que formam a rede cristalina, são repelidos 
pelos outros elétrons. No caso de um cristal, a equação de Schroedinger assume 


a forma 
Xi EE ++ +aajv+[E- vê |v| = 0 (Eg. 5.1) 


2m ôy; O 


que é uma generalização da Eq. 2.17. 

Na Eg. 5.1, os índices i e j são usados para distinguir as coordenadas dos 
diferentes elétrons e o somatório se estende a todos os elétrons do cristal. 

A energia potencial v(r, :)) pode ser decomposta em dois termos, um que só 
depende da posição do elétron i, Vir). e que corresponde à encrgia potencial 
do elétron “i no campo elétrico dos núcleos dos átomos do cristal, é um segundo 
termo, Vir, :) que corresponde à energia potencial do elétron i no campo clé- 
trico de todos os outros elétrons e que portanto depende da posição de todos 
os elétrons do cristal. 

Se o termo Vil) da energia potencial fosse desprezível, a Eq. 5.1 assumiria 


a forma 
h2[ 0? 0? 0? —+ 
, ns ge E- Vtr. et) q. 5. 
Et E + a + 5=|* +[ | ul 0 (Eq. 5.2) 


Vamos agora demonstrar que à solução da Eq. 5.2 pode ser escrita na forma 


Wa Px PD) VA) Ur) lr), (Eg. 5.3) 
onde ip ,(r;), Wolro).-..s Wir) ---« Wlr,) satisfazem a equações do tipo 


h a2 a2 az as e 
olaa + | Wu) + E vi =0 k=12..0,nº (Eq.54) 
a 


2m| xt Oy 


Para isso vamos substituir a Eg. 5.3 na Eq. 5.2. Fazendo isso, temos 


h2[ 0º 02 02 = ». Es 
Stoa + Fm + | Ve) par) br)... 


lr) + [E “6 | [7 CA A(O PRACA AA) =0. (Eg.5.5) 


Aplicação da mecânica quântica a um cristal 
Ora, é evidente que 


2 Fa 22 a ã LA n Ev 
RE Ee tr” Ay? EA | Walr,) Wolrs) aces dr; ... vAr,) = 
1 


2m[ôx7 Oy 
= - io A mibr es 
E — blro).. vdr) ... ui)l5a E dy + | Wlr,) 
ola” ê dE + am Ma ANA AA RAÇA 


= e pr) 7 pa ... AT) ... Ur) x 
02 02 0? EA 
x [5 = x + | Walr's). 
Ou, de modo geral, 


2 2 
E ló + é Ma na =| vo (rj): Walro) NAÇÃE EE TAÇA) — 


dx Oy? 


= Er ATUA ER : Via) Wir iCaa) o ATA) x 


o? o? Fi - 
x [5 + Era + =| Wáro. 


Ora, de acordo com a Eqg. 5.4. 
hº[ 0? 0 p (dé - - 
Em o” | ATA - [nã - E] wlr,). 
Substituindo a Eq. 5.7 na dg 5.6, temos 


h? 0? Fal q? = a p 3 
2m [53 à ôyê Ci =] bylry) Walro) ses bAr,) = ww Wdr,) E 


= [Mi ES dO bee DUTO Dis Pe) bro) = 


= [Vir = EJA Paseoosra) 
Substituindo a Eq. 5.8 na Eq. 5.2. temos 


É [ros W + [E- no v=0, 


equação que é satisfeita para 
E=53 E, 
Assim a solução da Eg. 5.2 é 
Wo Pcs rd = EDP). r), 
ca energia total é dada por 
E=E,+E, ++ E, 
onde E,, E,..... E, são as energias individuais dos elétrons. 


a3 


(Eq. 5.6) 


(Eq. 5.7) 


(Eq. 5.8) 


(Eq. 5.9) 


(Eq. 5.10) 


(Eq. 5.11) 
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O uso da Eq. 5.10 implica que a probabilidade de encontrar o elétron i 
no ponto r,. |y,|?, depende apenas de F,, isto é, da coordenada do elétron i em 
relação a uma origem arbitrária. Naturalmente, se existe uma interação entre 
o elétron i e o elétron j, a densidade de probabilidade do elétron i envolve tam- 
bém as coordenadas do elétron j. Mesmo assim, na aproximação de “um elétron”, 
muito usada na prática, supomos que uma função produto como a da Eq. 5.10 
é uma boa aproximação para a função de onda. 

Nessa aproximação de “um elétron”, a equação de Schroedinger se reduz a 
hº [ 0º E 
2m|ôx? dy? 
onde V(7) é o potencial no qual o elétron está se movendo (potencial dos núcleos 
dos átomos do cristal) e o indice k serve para distinguir as diferentes soluções. 
A Eg. 5.12 se aplica a qualquer elétron do cristal, isto é, o somatório da Eq. 5.1 
foi transformado em um sistema de n equações idênticas, quando desprezamos 
o termo Vir). Assim podemos estudar o movimento de cada elétron como se 
ele estivesse sozinho na rede cristalina (núcleos) e dai o nome de aproximação 
de “um elétron”. 

'* Devido à regularidade da disposição dos núcleos em um cristal, o potencial 
V(r) é periódico no espaço. Essa periodicidade, por sua vez, impõe certas res- 
trições à forma das soluções da Eq. 5.12. Essas restrições são expressas pelo 
Teorema de Bloch, que é demonstrado no Apêndice A. 

O Teorema de Bloch diz que a solução mais geral da equação de Schroedinger 
para um elétron em um cristal é da forma 


pur) = Eu, (0), (Eq. 5.13) 


onde u, (7) é uma função com a mesma periodicidade espacial que a rede cris- 
talína. 

A única restrição imposta à função (”) pelo Teorema de Bloch é que 
ela deve ter a periodicidade da rede cristalina. Sua forma detalhada depende 
da forma do potencial V(r) da Eg. 5.12. Vamos agora examinar o caso simples de 
um cristal de uma dimensão apenas e verificar qual a forma assumida pela função 


tg (1). 


2 
+ [ve 0) + LE VW (DO, (Ea. 517) 


5.1, CRISTAL DE UMA DIMENSÃO 


Consideremos um átomo ionizado e um elétron. O potencial do elétron 
no campo do átomo ionizado é dado por 


q 
V = Eh (Eq. 5.14) 
onde r é a distância entre o elétron e o núcleo do átomo ionizado. Esse po- 
tencial está representado na Fig. 51. 

Consideremos agora um “cristal unidimensional”, isto é o conjunto for- 
mado por uma linha de átomos ionizados e uma série de elétrons. O potencial 
de qualquer um desses elétrons no campo dos átomos ionizados está repre- 
sentado na Fig. 5.2. 
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Vtn 


Figura S.1 Energia potencial de um elétron 
nas vizinhanças de um átomo ionizado 


Figura 5.2 Energia potencial dos elétrons no 
interior de um cristal unidimensional 


Como para um potencial desse tipo não existe solução analitica da equação 
de Schroedinger, vamos substituí-lo por um potencial que, em certas condições, 
se parece com o potencial do cristal, e para o qual a equação de Schroedinger 
apresenta uma solução simples. Nesse modelo, mostrado na Fig. 5.3 e sugerido 
por Krônig e Penney (1931), o potencial periódico consiste em um número 
infinito de barreiras de potencial de largura b, com regiões de potencial zero 
de largura a entre as barreiras. A altura das barreiras é V, (Fig. 5.3). 

A equação de Schroedinger do sistema será 

242 
m a +[E-Vo)]Jy =0, (Eg. 5.15) 


onde E é a energia do elétron, m a massa do elétron e V(x) a energia potencial. 
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Vix) 


Figura 5.3 Energia potencial dos elétrons no 
modelo de Kroônig-Penney 


——>» 
x 


Na região O<x <a, V=0,a solução da Eq. 5.15 assume a forma 
Wylx) = Ae” + Be” Ps, (Eg. 5.16) 
onde fp = /2 Em/h. Na região 0>x>-b, V=V, e E-V,<O; assim, a solução 
assume a forma 
Wolx) = Ce + De", (Eg. 5.17) 


onde « = ./ 2m(K, — E)/h. Os coeficientes 4, B,C e D são relacionados entre si 
pelas condições de fronteira. Como o valor de E- V(x) é finito na fronteira, a 
derivada dy/dx e a própria função y devem ser continuas na fronteira x = 0. 
Em outras palavras, 


A+B=C+D, 
iB(A-B) = «C-D). 
No caso de um cristal de uma dimensão, a Eq. 5.13 se torna 
Wlx) = e ulx). (Eg. 5.19) 
Como (x) tem a mesma periodicidade que a rede cristalina. 
ul) = ulx + a + b). 


(Eq. 5.18) 


Então 

lx t+a+ b) q eita + bpikx (x) E, eita tbhy(x). 
Assim, à solução na região a < x <a+b é da forma 

bslx + a + b) = ey (x) 
Val) = €HPy dx—a—b). 

As condições de fronteira em x = a são portanto 

Aeita + Be” ts ao (Ce + De""e'tta tb), 

iB(Ae' -Be” iliay a d(Ce * Tá De" eitta+ b), 


Para uma solução não-trivial, o determinante dos coeficientes deve ser 
nulo e portanto (Prob. 5.1) 
2 p2 


F senh ab sen fa + coshabcos Ba = cosk(a + b). (Eq. 5.21) 


ou: 


(Eq. 5.20) 
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Em principio, a Eq. 5.21 nos dá a energia para qualquer valor de k, ou o 
valor de k associado a uma energia particular E. Embora a Eg. 5.21 seja com- 
plicada demais para permitir uma solução analítica, podemos observar que a 
um valor de k correspondem vários valores diferentes de energia. 

Para simplificar o problema, vamos fazer b tender a zero e W, tender a 
infinito, de modo que o produto bV, se mantenha finito. Nesse caso, a Eq. 5.21 
se reduz a 


p a + cos Ba = cos ka (Eq. 5.22 
d 
onde abas 
= 3 q 


Na Fig. 5.4, desenhamos o lado esquerdo da Eq. 5.22 em função de fa. 
Como [cos ka| < 1, só são permitidos os valores de fa para os quais o lado 
esquerdo da Eq. 5.22 esteja entre —] e 1. Dentro dessa faixa, para cada valor de 
ba existe um valor correspondente de ka (Fig. 5.4). 


Figura 5.4 Representação gráfica da 
função: 


— Ssenfa 
F(fa) = Po + cos fia 


Como f está diretamente relacionado à energia E do elétron, é possivel 
obter, a partir da Fig. 5.4, um gráfico de E em função de ka. Esse gráfico apa- 
rece na Fig. 5.5. 

Os trechos 4B, CD, EF, etc. das Figs. 5.4 c 5.5 são chamados de fuixas 
permitidas de energia porque os elétrons que possuem essas energias podem pro- 
pagar-se livremente no cristal. Por outro lado, os elétrons com energias entre 
Be C,Ce D, etc., devem ter um k imaginário para satisfazerem à Eq. 5.22. Assim 


Figura 5.5 Energia E dos elé- 
trons em função de ka 
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o movimento dos elétrons é muito atenuado e eles encontram grande dificuldade 
para se propagarem no cristal. Essas faixas de energia são chamadas de faixas 
proibidas de energia. 


5.2. CONDUTORES, ISOLANTES, SEMICONDUTORES 


O que dissemos no parágrafo anterior pode ser generalizado para um cristal 
tridimensional. Assim ficamos sabendo que um elétron em um cristal só pode 
se propagar se sua energia estiver dentro de uma das faixas permitidas. 

Em um cristal, os elétrons tendem a ocupar sempre os estados de menor 
energia. Assim as faixas permitidas mais baixas estão sempre cheias. Suponhamos 
que o número de elétrons no cristal seja suficiente para encher completamente 
um certo número de faixas permitidas, deixando as outras vazias. Nesse caso 
os elétrons do cristal, quando submetidos a um campo elétrico, não podem se 
propagar, pois qualquer movimento aumentaria sua energia, colocando-os den- 
tro da faixa proibida. Nesse caso dizemos que o cristal é um isolante. Se, por 
outro lado, os elétrons enchem apenas parcialmente uma das faixas permitidas, 
os elétrons podem mover-se com facilidade sob a ação de campos elétricos, 
passando para estados de energia ligeiramente maior. Nesse caso, o cristal pode 
ser um condutor, um semicondutor ou um semimetal. Em um condutor, as faixas 
permitidas parcialmente cheias possuem um grande número de elétrons. Em 
um semicondutor ou semimetal, as faixas parcialmente cheias possuem um pe- 
queno número de elétrons ou estão quase completamente cheias; em ambos os 
casos, a condução de corrente é menor que nos metais. A diferença entre um 
semicondutor e um semimetal é que na temperatura do zero absoluto um semi- 
condutor se torna isolante, enquanto um semimetal continua a conduzir corrente. 


EXERCÍCIOS 


5.1. Forme um determinante com os coeficientes 4,B,C e D das Egs. 5.18 e 5.20, Em 
seguida, forme um novo determinante somando € subtraindo as colunas da seguinte 
forma: 1! =(1-2/2i;2'=(1+2)/2;3=(3+ 4)/2;4' =(4-3)/2, onde os números 
1,2,3 e 4 representam as colunas do antigo determinante e os números 1',2,3' e 4 
representam as colunas do novo determinante. Prove a Eq. 5.21 expandindo o novo 
determinante. 

5.2. Faça um gráfico de F(fa) = (p/Ba) sen Ba + cos fa, Eg. 5.22, em função de fa, para 
p = 37/2, como na Fig. 5.4. Em seguida, use o gráfico para obter os valores de ka 
para vários valores de fa. Use os valores de ka e fia para traçar o gráfico de 
E = P“a?/n? em função de ka (Fig. 5.5). 

5.3. Explique a diferença entre os níveis permitidos de energia em um átomo e os níveis 
permitidos de energia em um cristal. 
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6 Propriedades básicas dos materiais 
semicondutores 


Como vimos no capitulo anterior, os semicondutores se caracterizam por 
apresentarem, à temperatura ambiente, uma resistividade intermediária entre a 
dos metais e a dos isolantes. Na prática, os semicondutores mais usados são 
o germânio e o silício, dois elementos cuja estrutura cristalina já foi descrita 
no Cap. 4. 

Neste capítulo vamos examinar mais de perto as propriedades que fazem 
do germânio e do silício materiais semicondutores. 


6.1. ESTABILIDADE DAS SUBSTÂNCIAS CRISTALINAS 


A Mecânica quântica permitiu explicar a estabilidade das substâncias cris- 
talinas à luz da estrutura eletrônica dos átomos. Para compreendermos por 
que os átomos se unem para formar um cristal, consideremos primeiramente 
o caso dos gases raros. 

Os gases raros são os elementos mais estáveis da natureza; praticamente 
não participam de reação química alguma. Ora, examinando a estrutura ele- 
trônica dos gases raros (Tab. 6.1), notamos uma peculiaridade; com exceção 
do hélio, a última camada de todos os gases raros contém exatamente 8 elétrons. 


Tabela 6.1 Estrutura eletrônica 
dos gases raros 


Hélio : Kº 

Neônio : K?Lº 

Argônio : K?L8M$ 
Criptônio: K?LEM!SNS 
Xenônio : K2LSM'SNISOS 
Radônio : KZL“MISNSIO!S p* 


Na verdade, é possível demonstrar que não se trata de uma coincidência; 
para ser quimicamente estável, um átomo de número atômico maior que 2 deve 
ter oito elétrons na última camada (camada de valência). É por isso que todos 
os elementos, exceto os gases raros, são quimicamente ativos: unindo-se em 
moléculas ou cristais, os átomos dos elementos ativos redistribuem seus elé- 
trons de tal forma que cada átomo fica efetivamente com & elétrons na camada 
de valência. 

Assim. por exemplo, o sódio e o cloro possuem as seguintes estruturas 
atômicas: 

Sódio: K?L*M', 
Cloro: KºLºM", 

Quando esses elementos se unem para formar o cloreto de sódio, cada 

átomo de sódio cede um elétron a um átomo de cloro. Assim a rede cristalina 
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do cloreto de sódio é constituida na realidade por ions de cloro e de sódio, com 
a seguinice estrutura eletrônica: 


fon de sódio: K2LS, 
fon de cloro: K?LSMº, 


Repare que tanto os íons de sódio como os ions de cloro possuem 8 elétrons 
na camada de valência. 

Esse tipo de ligação, em que um átomo cede um ou mais elétrons a outro 
átomo, é chamado de ligação iônica. 


6.2. LIGAÇÕES COVALENTES 


Nos cristais de germânio e silício, os átomos são unidos por um tipo mais 
complexo de ligação, a chamada ligação covalente. Para compreendermos me- 
lhor a natureza da ligação covalente, consideremos o tetraedro formado pelos 
átomos mais próximos de um átomo qualquer de um cristal de germânio ou 
silício (Fig. 6.1). 

Como vimos no Cap. 3, tanto o silício como o germânio possuem 4 elétrons 
na camada de valência. De acordo com a discussão precedente, é de se esperar 
que quando os átomos de um desses elementos se unem para formar um cristal, 
os elétrons sejam redistribuídos de modo que cada átomo fique com 8 elétrons 
de valência. 

No caso de uma estrutura como a da Fig. 6.1, uma redistribuição do tipo 
ligação iônica seria obviamente impossível, já que todos os átomos são do mesmo 
elemento. O que acontece é que cada elétron da camada de valência passa a 
pertencer simultaneamente a dois átomos. Esse tipo de ligação é chamado de 
ligação covulente. 

Vamos esclarecer melhor a situação. O átomo A da Fig. 6.1 possui 4 elétrons 
na camada de valência. Cada um desses elétrons pertence também a um dos 
átomos vizinhos. Do mesmo modo, um elétron de cada átomo vizinho pertence 
também ao átomo A. Assim o átomo A passa a dispor de oito elétrons em sua 
camada de valência. O mesmo ocorre para todos os outros átomos do cristal. 


Si al 
Si 
E 
865; B 
Figura 6.1 Átomo de silício localizado Figura 6.2 Representação esquemática 
no centro do tetracdro formado pelos das ligações covalentes de um cristal de 
4 átomos vizinhos germânio ou silício 


(SA estrutura cristalina do germânio e do silício já foi discutida no Cap. 4 
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Na Fig. 6.2 as ligações covalentes estão representadas por segmentos de 
reta que unem os átomos envolvidos. Cada segmento de reta representa um 
dos elétrons que participam da ligação. 

Embora permita uma visualização mais precisa da estrutura do cristal, a 
representação tridimensional da Fig. 6.2 é desnecessariamente complexa. Para 
estudarmos as propriedades básicas dos semicondutores, podemos usar uma 
representação bidimensional, como a da Fig. 6.3. 


| | poa | létrons da Ligação 
Covalente 


Figura 6.3 Representação bidimen- = 0===O 
sional de um cristal de Si ou Ge | | 


Repare que a representação bidimensional reproduz corretamente as po- 
sições relativas dos átomos e ligações covalentes, mas que os ângulos entre as 
ligações são alterados. Entretanto esses ângulos não são importantes para as pro- 
priedades que vamos estudar. 


6.3. ELÉTRONS E BURACOS 


Como os elétrons de valência dos semicondutores formam pares covalentes, 
eles não são normalmente livres como no caso dos metais. Na verdade em tem- 
peraturas muito baixas, próximas do zero absoluto, em que as ligações cova- 
lentes se mantêm praticamente intactas, o cristal se comporta como um isolante. 

Entretanto à medida que a temperatura cresce, algumas ligações covalentes 
se rompem, devido à agitação térmica”), 

Na Fig. 6.4 vemos um cristal em que uma das ligações covalentes deixou 


de existir. 
O eletron livre 
Ligação Covalente 
Incompleta 


Figura 6.4 Quebra de uma ligação 
covalente 


LO E ——e] (8 E 


'siFenômeno conhecido como geração térmica e estudado no parágrafo 6.3.1 
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Em consegiência, um dos elétrons ficou livre para vagar por todo o cristal. 
e ficou um “vazio” em uma das ligações covalentes do cristal. A ausência do 
elétron na ligação covalente está representada por uma linha tracejada. Esse 
“vazio” é chamado de buraco. Os buracos também podem servir para o trans- 
porte de cargas, pelo mecanismo que explicaremos em seguida. 

Quando uma ligação covalente está incompleta, isto é, quando existe um 
buraco, o elétron de um átomo vizinho pode deixar sua ligação covalente e vir 
preencher esse “vazio”. Com isso, entretanto, o elétron deixa um buraco em 
sua posição original, que por sua vez pode ser preenchido por outro elétron. 
Assim um buraco pode “caminhar” por um cristal, na direção contrária à do 
movimento dos elétrons que procuram “enchê-lo”. 

No que se refere ao transporte de cargas, o buraco se comporta como uma 
carga positiva, com o mesmo valor absoluto da carga do elétron. Na verdade, 
para simplificar o estudo, podemos considerar o buraco como uma entidade 
real; é preciso ter sempre em mente, entretanto, que se trata de um artifício e 
que os buracos não passam de “vazios” nas ligações covalentes do cristal. Outro 
ponto que vale a pena destacar é o seguinte: o movimento dos elétrons e o mo- 
vimento dos buracos se realiza de modo diferente — enquanto os elétrons ca- 
minham no cristal de forma mais livre, os buracos só podem passar de ligação 
covalente para ligação covalente. É por isso que a velocidade com que os elé- 
trons se movem no interior do cristal é em geral maior do que a dos buracos'*). 


6.3. GERAÇÃO E RECOMBINAÇÃO 


6.3.1. Geração térmica 

Definimos como geração térmica o processo de criação de pares elétron- 
“buraco em um cristal pela rutura de ligações covalentes devido à agitação 
térmica da rede cristalina. À medida que a temperatura de um cristal aumenta, 
a partir do zero absoluto (0 K), os átomos da rede começam a vibrar em torno 
de suas posições ideais, fazendo variar as distâncias entre os átomos e os ângulos 
entre as ligações covalentes. 

Quanto maior a temperatura, maior a amplitude das oscilações. Parte da 
energia das oscilações é transferida para os elétrons de valência dos átomos 
da rede. Alguns dos elétrons recebem energia suficiente para abandonarem a 
ligação covalente de que participam. Naturalmente a rutura de uma ligação 
covalente resulta na criação de um elétron livre e um buraco. 

Para descrevermos a criação de pares elétron-buraco pela agitação térmica 
da rede, definimos a chamada taxa de geração térmica, G(T), que representa o 
número de pares elétron-buraco criados no cristal, por unidade de volume e 
por unidade de tempo. Para um determinado material semicondutor, a taxa de 
geração térmica é função apenas da temperatura. 


6.3.2. Recombinação 


Definimos como recombinação a associação de um elétron a um buraco, 
com a reconstituição de uma ligação covalente e a liberação de uma certa quan- 


HO parâmetro conhecido como mobilidade, que será definido no Cap. 7, relaciona a velo- 
cidade dos elétrons (ou buracos) ao campo elétrico aplicado. Em geral, a mobilidade dos elétrons 
é maior que a dos buracos 
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tidade de energia. Essa energia pode ser irradiada (emissão de fótons) ou absor- 
vida pela rede cristalina (emissão de fônons). 

Embora os mecanismos da recombinação sejam bastante complexos, o fe- 
nômeno pode ser descrito por uma equação fenomenológica bastante simples. 
Como o processo de recombinação depende da presença simultânea de um elétron 
e um buraco em um mesmo ponto do cristal, a taxa de recombinação, R, é pro- 
porcional às concentrações de elétrons c buracos no cristal, isto é,* 

R = np (Eq. 6.1) 
onde R é a taxa de recombinação, n e p são as concentrações de elétrons e bu- 
racos, respectivamente, e x um parâmetro que depende da temperatura e do 
mecanismo de recombinação considerado. 

A Eq. 6.1 é chamada de Lei da Ação das Massas, por analogia com a lei 
que descreve a evolução das reações químicas. 


6.3.3, Equilibrio térmico 


Dizemos que um sistema está em equilibrio térmico quando sua tempe- 
ratura é constante e uniforme. 

Como vimos nos parágrafos anteriores, as concentrações de elétrons e bu- 
racos em um semicondutor dependem de dois processos competitivos, a geração 
térmica e a recombinação, Para uma dada temperatura, a taxa de geração tér- 
mica G(T) tem um valor constante. Assim se o semicondutor estiver em equi- 
librio térmico, o número de pares elétron-buraco gerados por unidade de tempo 
será constante. Ora, o número de elétrons e buracos no cristal não pode aumentar 
indefinidamente, já que, como vimos, a taxa de recombinação é diretamente 
proporcional às concentrações de elétrons e buracos. Em conseqiiência, quando 
mantemos um semicondutor a uma temperatura constante durante um tempo 
suficientemente longo, as concentrações de elétrons e buracos assumem valores 
tais que a taxa de recombinação se torna exatamente igual à taxa de geração 
térmica: 

GT) = R = np. (Eg. 6.2) 

Naturalmente a Eq. 6.2 só é válida para um semicondutor em equilíbrio 
térmico. 


6.4. SEMICONDUTORES INTRÍNSECOS 


Até agora, supusemos implicitamente que o material semicondutor a que 
nos referiamos era um semicondutor puro, isto é, que não continha átomos de 
outros elementos (impurezas). Quando um semicondutor não contém impurezas 
em quantidade suficiente para alterar-lhe as propriedades, dizemos que se trata 
de um semicondutor intrínseco. 

Em um semicondutor intrinseco, o número de elétrons é igual ao número 
de buracos, pois cada vez que uma ligação covalente é rompida ou formada, 
um par elétron-buraco é gerado ou eliminado no processo. Usando o índice i 
para indicar que se trata de um material intrínseco, podemos escrever 

n=P=h. (Eg. 6.3) 

SA taxa de recombinação é definida como o número de pares que se recombinam por uni- 

dade de volume e por unidade de tempo 
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O parâmetro n, é chamado de concentração intrinseca dos portadores e é 
função da temperatura e do material considerado. A concentração intrínseca 
dos portadores de um semicondutor será calculada no Cap. 8. 


6.5. SEMICONDUTORES EXTRÍNSECOS OU IMPUROS 


A adição de pequenas quantidades de outras substâncias a um semicondutor 
pode modificar consideravelmente as propriedades do material. Essas impurezas 
podem substituir os átomos da rede cristalina (impurezas substitucionais) ou 
ocupar posições vagas entre os átomos da rede (impurezas intersticiais). 

Na prática, as propriedades dos semicondutores são controladas através 
da adição de certas impurezas em concentrações cuidadosamente medidas. As 
impurezas substitucionais mais usadas são certos elementos pentavalentes (arsê- 
nio, fósforo, antimônio) e trivalentes (boro, alumínio, gálio, índio). As Figs. 6.5 
e 66 mostram as representações esquemáticas dos átomos dessas impurezas. 


- A — 
Je N 
f N ; - ; ? 
Figura 6.5 Representação esquemática do átomo de 
/ 4 uma impureza pentavalente. Estão indicados os cinco 
| |] elétrons de valência e a carga positiva equivalente asso- 
N / ciada à parte central do átomo. Nessa categoria estão 
N / o fósforo, o arsênio, o antimônio, etc. 
bes ” 
a 
E ta, 
/ x . 
/ N Figura 6.6 Representação esquemática do átomo de 
uma impureza trivalente. Estão indicados os três elétrons 
| (+3) ] de valência c a carga equivalente associada à parte 
1 central do átomo. Nessa categoria estão o boro, o alu- 
N / mínio, o gálio, o índio, etc. 
E Pá 
sua Eae 


Vamos analisar primeiro o caso em que acrescentamos a um semicondutor 
puro (por exemplo, o silício) uma pequena quantidade de uma impureza pen- 
tavalente (por exemplo, o arsênio). Nesse caso, alguns átomos de arsênio ocupam 
o lugar dos átomos de silício no cristal. Entretanto um átomo de arsênio só 
tem quatro átomos vizinhos de silício para estabelecer ligações covalentes. 
Assim apenas quatro átomos de valência do arsênio ficam retidos nas ligações 
covalentes; o elétron que sobra fica praticamente livre para vagar pelo cristal, 
indo fazer companhia aos elétrons livres já existentes, provenientes da quebra 
das ligações covalentes dos átomos de silício. Na Fig. 6.7 vemos o que acontece 
no interior do cristal. 

Como os átomos das impurezas pentavalentes cedem elétrons ao cristal, 
essas impurezas são chamadas de impurezas doadoras ou impurezas tipo N (a 
letra N é a inicial de negativa, referindo-se à carga do elétron). 

Quando acrescentamos ao silício uma impureza trivalente — índio, por 
exemplo —, ocorre exatamente o oposto, pois como cada átomo de índio só 


Propriedades básicas dos materiais semicondutores 55 


Figura 6.7 Impureza pentavelmente adi- 
cionada à um cristal de Si ou Ge. O elétron 
extra fica praticamente livre para se des- 
locar no cristal 


possui três elétrons de valência e pode estabelecer ligações covalentes com 
quatro átomos vizinhos de silício, uma dessas ligações covalentes fica incom- 
pleta, isto é, fica “faltando” um elétron. Ora, como já vimos, a “falta” de um 
elétron em uma ligação covalente pode ser considerada como a “presença” de 
um buraco, À Fig. 6.8 mostra o efeito da adição de uma impureza trivalente a 
um cristal de silício ou germânio, 


Figura 6.8 Impureza trivalente adicionada 
a um cristal de Si ou Ge. A ligação cova- 
lente incompleta corresponde a um buraco 
introduzido pela impureza 


Como os átomos das impurezas trivalentes cedem buracos ao cristal (re- 
tiram elétrons do cristal), essas impurezas são chamadas de impurezas aceitadoras 
ou impurezas tipo P (a letra P é a inicial de positiva, referindo-se à carga do 
buraco). 

Os elétrons livres e buracos de um semicondutor são chamados de por- 
tadores, pois é através deles que se realiza qualquer transporte de cargas de um 
ponto para outro do material. Como vimos, em um semicondutor intrínseco 
(puro) o número de elétrons é igual ao número de buracos. 

Nos semicondutores tipo N, isto é, que contêm impurezas doadoras. pen- 
tavalentes, o número de elétrons é maior, naturalmente, que nos semicondutores 
intrínsecos. Além disso, o número de buracos é menor, pois havendo excesso 
de elétrons, os buracos têm maior oportunidade de se recombinarem. Nos se- 
micondutores tipo P, contendo impurezas aceitadoras, trivalentes, ocorre o 
oposto; como a concentração de buracos é maior que a concentração de elétrons, 
o número de elétrons é menor que no material intrínseco. 
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O que dissemos anteriormente é demonstrado matematicamente no Cap. 8, 
onde provamos que 


n'p= ni. (Eq. 6.4) 


A Eq. 6.4 é válida para qualquer semicondutor, seja intrinseco, tipo N ou 
tipo P. Essa equação nos diz simplesmente que o produto das concentrações 
de elétrons e buracos em um semicondutor é igual a uma constante, n?. Como 
vimos, a grandeza n, pode ser definida como a concentração intrínseca de por- 
tadores para o material semicondutor em questão. 

Assim o fato de acrescentarmos uma impureza a um semicondutor intrin- 
seco não só aumenta o número total de portadores, mas também altera a pro- 
porção entre elétrons e buracos. Nos materiais tipo N, os elétrons são chamados 
de portadores em maioria e os buracos de portadores em minoria. Nos materiais 
tipo P, os papéis se invertem. 

Como já vimos, para transformar um elétron firmemente preso a uma 
ligação covalente em elétron livre é necessário fornecer-lhe energia. Essa energia 
pode ser fornecida, por exemplo, sob a forma de calor. Assim quando dizemos 
que a agitação térmica rompe ligações covalentes, na verdade estamos afirmando 
que a energia térmica eleva alguns elétrons de valência a um nível mais alto 
de energia, tornando-os livres. Por que não acontece o mesmo nos isolantes? 
Por que nos metais, mesmo a temperaturas baixissimas, existem muitos elétrons 
livres? Essas perguntas foram respondidas pela Mecânica quântica, como ve- 
remos no próximo parágrafo. 


6.5. NÍVEIS DE ENERGIA NOS SEMICONDUTORES 


No capitulo anterior, verificamos que em um sólido cristalino existem 
certas faixas de energias permitidas, separadas por faixas de energias proibidas, 
isto é, energias que os elétrons não podem assumir. No nosso caso, o que nos 
interessa são apenas duas faixas seguintes: a faixa dos elétrons firmemente 
presos às ligações covalentes, chamada de faixa de valência, e a faixa dos elétrons 
de valência que se tornaram clétrons livres, e que é chamada de faixa de condução. 

Em um metal, não existem faixas proibidas entre a faixa de valência e a 
faixa de condução; as duas se superpõem parcialmente. Isso significa que em 
um metal, mesmo a temperatura muito baixa, há elétrons com suficiente energia 
para serem elétrons livres. Isso explica por que os metais são condutores. Nos 
isolantes, por outro lado, a faixa de valência está completa e é separada da faixa 
de condução por uma faixa proibida. Isso significa que para termos um elétron 
livre em um isolante, é preciso fornecermos a um elétron da faixa de valência 
suficiente energia para que ele passe para a faixa de condução. Entretanto a 
largura da faixa proibida nos isolantes é muito grande, isto é, a energia neces- 
sária para colocar um elétron na faixa de condução é elevada. Assim os isolantes 
praticamente não apresentam elétrons livres. 

As Figs. 6.9(a) e 6.9(b) mostram como as bandas de energia estão dispostas 
nos metais e nos isolantes. ; 

Nos semicondutores também existe uma faixa proibida, mas cla é muito 
mais estreita que no caso dos isolantes. A temperaturas muito baixas, os semi- 
condutores se comportam como isolantes, pois seus elétrons não têm energia 
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Figura 6.9 Faixas de energia (a) nos metais; (b) nos isolantes 
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suficiente para passarem para a faixa de condução. Nesse caso as ligações co- 
valentes se mantêm intactas. 

Quando aumentamos a temperatura, entretanto, alguns elétrons adquirem 
energia suficiente para transporem a faixa proibida e passarem para a faixa de 
condução, deixando buracos na faixa de valência. 

As Figs. 6.10a), (b) e (c) mostram o que acontece nos três tipos de semi- 
condutores. 


Condução | 


E to 


Figura 6.10 Faixas de energia nos semicondutores 


A Fig. 6.10(a) corresponde ao caso de um semicondutor intrinseco 
(tn =p =n). Como vimos, tanto os elétrons da faixa de condução como os 
buracos da faixa de valência contribuem para o transporte de carga. 

Na Fig. 6.10(b) vemos a situação em um semicondutor do tipo N. Repare 
que a inclusão da impureza doadora deu origem a um nivel permitido de energia 
(E,) dentro da faixa proibida, próximo da faixa de condução. É nesse nível que 
estão os elétrons em excesso dos átomos doadores. Assim, à temperatura ambi- 
ente, quase todos esses elétrons adquirem energia suficiente para passarem para 
a faixa de condução e se tornarem “livres”. Com isso os átomos dos doadores 
ficam positivamente carregados, como está indicado na Fig. 6.10(b). A situação 
da Fig. 6.1Xb) é a mesma que a da Fig. 6.7. 

Nos cristais do tipo P, em que as impurezas são aceitadoras, o nível per- 
mitido (E,) está próximo da faixa de valência. Assim, os elétrons da faixa de 
valência podem passar com grande facilidade para o nível dos aceitadores, 
deixando buracos na faixa de valência. Com isso, os accitadorcs ficam carregados 
negativamente, como está indicado na Fig. 6.10(c), que equivale à Fig. 6.8. 

A Tab. 6.2 mostra as energias de ionização das impurezas mais comuns, 
expressas em elétrons-volt (Conwell, 1958). 
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Tabela 6.2 Energias de ionização em elétrons-volt (eV) 


Semicondutor 
Impurezas — > — — 
Silício Germânio 
D P 0,044 0,0120 
0) 
A As 0,049 0,0130 Condução 
D E, 
0) Sb 0,039 0,0096 
R 
A Li 0,033 0,0010 
s 
A B 0,045 0,010 
c 
E 
I Al 0,057 0,010 
T 
A 
D Ga 0,065 0,011 
(6) 
R 
A In 0,160 0.011 
S 


6.6. MASSAS EFETIVAS DOS ELÉTRONS E BURACOS 
EM UM SEMICONDUTOR 


Como vimos nos parágrafos anteriores, existem dois tipos de portadores 
de carga em um material semicondutor, isto é, os elétrons e os buracos. Evi- 
dentemente o estudo do movimento desses portadores sob a ação de forças 
externas (produzidas, por exemplo, pela aplicação de campos elétricos) é extre- 
mamente complexo, pois os portadores estão também submetidos à influência 
dos átomos da rede cristalina, dos átomos das impurezas e dos outros portadores. 

Uma maneira simples de resolver o problema é definir massas fictícias 
para os elétrons e buracos tais que a segunda lei de Newton (F = ma) seja sa- 
tisfeita. Em outras palavras, escolhemos valores tais para as “massas” do elétron 
e do buraco que as relações 

F,=m,a (Eq. 6.5) 
F,=m,a, (Eg. 6.6) 
sejam satisfeitas. 

Nas equações anteriores, F é a força externa aplicada, a, e a, são as ace- 
lerações dos elétrons e buracos c m, e m, são as chamadas massas efetivas dos 
elétrons e buracos, respectivamente!'*! (Shockey, 1950). 


(HA massa efetiva do buraco é representada pelo simbolo m, porque o buraco em inglês é 
chamado de hole 
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Na verdade o uso do conceito de massa efetiva só se justifica porque em 
um semicondutor as forças internas do cristal são várias ordens de grandeza 
maiores que as forças externas normalmente aplicadas. Em consequência, o 
uso de uma massa efetiva reflete adequadamente o efeito do cristal sobre o mo- 
vimento dos portadores. 

Repare que ao introduzirmos o conceito de massa efetiva estamos supondo 
que o cristal foi substituído por um meio homogêneo, no qual os elétrons e 
buracos se deslocam com as respectivas massas efetivas, satisfazendo às leis da 
Mecânica clássica e sujeitos apenas a forças externas. 

As Figs. 6.1] e 6.12 ilustram melhor essa idéia. 


Figura 6.11 Elétron de massa real m, deslo- 
cando-se no interior de um cristal semicon- 
dutor sob o efeito de um campo elétrico 
externo. Devemos aplicar as leis da Mecânica 
quântica 


q 
Figura 6.12 O cristal foi substituido por um Me 
meio homogêneo que não oferece resistência 
alguma ao movimento. A massa “real” foi subs- 
tituida pela massa efetiva. Podemos aplicar as 
leis da Mecânica clássica 


o — 


E 
Na prática, as massas efetivas são determinadas com o auxilio da resso- 
nâneia de ciclotron dos elétrons e buracos!” (Lax e Mavroids, 1960). 
Voltaremos ao conceito de massa efetiva no Cap. 8, quando calcularmos 
as concentrações de elétrons e buracos em um semicondutor. 


EXERCÍCIOS 


6.1. Verificamos a possibilidade da introdução de impurezas em um material semicon- 
dutor, transformando-o em um semicondutor tipo P ou tipo N, dependendo do 


HO nome “ressonância de ciclotron” provêm do fato de que os elétrons e buracos descrevem 
órbitas circulares no interior do cristal, sob o efeito de campos elétricos e magnéticos. O mesmo 
ocorre com as partículas carregadas no interior de um “ciclotron” (acelerador de partículas) 
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6.2. 


6.3. 
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tipo de impureza adicionada. Entende-se por solubilidade sólida de uma impureza a 
máxima concentração dessa impureza que pode ser “dissolvida” no silício (ou ger- 
mânio) sólido, a uma dada temperatura. 


a) Procure explicar fisicamente por que colocamos entre aspas o termo “dissolvida”; 

b) A máxima solubilidade do arsênio no silício é da ordem de 2 x 10?! átomos/cm* 
(Trumbore, 1960). Qual a máxima concentração (4) de As que pode ser utilizada 
em um dispositivo semicondutor de silício? 


Considere um material semicondutor em equilíbrio térmico, iluminado por um feixe 
de luz cuja frequência pode ser variada. Um fotodetetor é usado para medir a luz que 
atravessa o cristal. Pense cuidadosamente e procure explicar o que deverá indicar 
o detetor, à medida que a frequência da luz for variada. 

Considere um material semicondutor a uma dada temperatura T: Suponha que um 
feixe de luz constante é aplicado a esse semicondutor. Chamando de G(T) e G(1) 
as taxas de geração devido à temperatura e à iluminação, respectivamente, qual a 
equação de equilíbrio do sistema? 
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7 Distribuições estatísticas 


Já sabemos quais são os estados permitidos para os elétrons em um semi- 
condutor. Neste capítulo, vamos investigar qual é a probabilidade de que esses 
estados estejam ocupados. Aproveitaremos a oportunidade para examinar as 
distribuições estatísticas de um modo geral. 


7.1. NECESSIDADE DAS DISTRIBUIÇÕES ESTATÍSTICAS 


O melhor meio de ilustrar a necessidade das distribuições estatísticas para 
o estudo das propriedades da matéria é apresentar algumas perguntas-chave, 
conforme segue. 

a) Suponhamos que estamos estudando um sistema constituido por mo- 
léculas que não interagem entre si, a não ser através de choques (esse sistema 
é chamado de gás perfeito). A pergunta que temos a fazer é a seguinte: qual é a 
distribuição de energia das moléculas desse gás? Em outras palavras, quantas 
moléculas estarão em uma determinada faixa de energias? Qual a energia média 
das moléculas? 

b) Nos sólidos cristalinos, os elétrons de condução só podem ocupar de- 
terminados niveis de energia. Como os elétrons se distribuem entre esses níveis? 


Basta uma ligeira reflexão do leitor para verificar que perguntas como 
essas só podem ser respondidas estatisticamente. Existe um ramo da Física, 
chamado Mecânica estatística, que trata exclusivamente desses problemas. No 
restante deste capítulo vamos apresentar alguns resultados obtidos pela Me- 
cânica estatística. 

A primeira observação a fazer é que os resultados obtidos dependem das 
restrições impostas ao sistema em estudo. Por exemplo, para o gás perfeito, 
definido anteriormente, o número de estados permitidos é muito maior que o 
número de particulas presentes; assim, não há sentido em aplicarmos o Princípio 
de Exclusão de Pauli. Por outro lado, devemos aplicá-lo para calcularmos a 
distribuição de energia dos elétrons de condução em um sólido cristalino, já 
que nesse caso o número de elétrons é da mesma ordem de grandeza que o nú- 
mero de estados permitidos. Com essa idéia em mente, vamos estudar algumas 
distribuições importantes. 


7.2. SISTEMAS DILUÍDOS — DISTRIBUIÇÃO DE BOLTZMANN 


Sistemas diluídos podem ser definidos como aqueles em que o número de 
partículas é muito menor que o número de estados permissíveis. Um exemplo 
de sistema diluído é o gás perfeito. Vamos ver agora como as N partículas de 
um sistema diluído se distribuem entre os n estados permissíveis de energia. 

Nosso cálculo será bascado nas seguintes hipóteses: 

a) As partículas são distintas. O conceito de partículas distintas e idênticas 
é muito importante para a formulação estatística do problema, Quando dizemos 
que as partículas são distintas, isso significa que podemos localizar uma de- 


62 Introdução à física dos semicondutores 


terminada partícula e ter certeza de que se trata da partícula 4 e não de B,C ou 
outra partícula qualquer. Suponhamos, por exemplo, que o sistema consiste de 
apenas três partículas, 4, Be C. Se essas partículas são distintas, então existem 
três grupamentos distintos de duas partículas, ou seja AB, AC e BC. Grupamentos 
como AB e BA são considerados idênticos, porque contêm as mesmas par- 
tículas. Por outro lado, se as particulas fossem idênticas, não poderíamos dis- 
tinguir o grupamento AB do grupamento AC, por exemplo. A Fig. 7.1 ilustra 
essa diferença!” 


lasc] | o ] lararr| | o ] 
[o] ac] Lo] lara 
las] [c | |2 Part) |4 PART] 
Lc | [as] la part] [2 ART 


b 
| B | [ac | to) 
I BE | | A | Figura 7.1 Diferentes modos de distribuir três par- 
tículas entre dois estados. (a) particulas distintas. (b) 
A | BC | partículas idênticas. No caso (b), como não podemos 
distinguir entre as particulas, indicamos apenas o 
(a) número de particulas em cada estado 


b) O número total de partículas é conservado. Essa hipótese não se apli- 
caria, por exemplo, a um sistema contendo particulas instáveis. 

c) À energia total do sistema é conservada. Isso equivale a supormos que 
o sistema está livre de influências externas. 

Vamos agora distribuir as N partículas do sistema entre os n estados per- 
mitidos, de tal forma que N, partículas sejam alojadas no estado 1, N, no estado 
2 e assim por diante (Fig. 7.2). 


estado 1 estado 2 estado i estado n 
El ana Es 
N, part. N, part. N, part. N, part. 


N+N; +... +N,+... + N,=N = Número total de partículas. 
Figura 7.2 Distribuição de N partículas entre os estados 1,2.....n 


Vejamos inicialmente de quantos modos é possível alojar as particulas da 
forma indicada na Fig. 7.2. O número de modos de colocar N, partículas no 
primeiro estado é C4''**!. Após colocarmos as N, partículas no estado 1, restam 
N-N, partículas; existem CJ? ,, modos possíveis de colocar essas N - N, par- 


LO leitor deve observar que estamos permitindo que mais de uma partícula ocupe o mesmo 
estado, isto é, estamos ignorando o Princípio de Exclusão de Pauli 
*MICY! representa combinação de N elementos N, a N, 
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tículas no estado 2. Do mesmo modo, existem C4º ,, -w, modos possíveis de 
colocar as partículas restantes no estado 3, e assim por diante. 

O número total de modos possíveis de se colocar N, partículas no estado 1, 
N, no estado 2, etc., é dado por 


ua N “Na N ... Nu = 
W=CyxCy,xCy Sis X Dad Md v= 5a = 


C=Na 


N! (N=NO! (NON,=Nj)! 
2“ NAN=NO! O ASIN-N, No! NAMN-N,-No NJ! * 
(N=Ni=NoeN od)! o! andio 
“NAN, N$ =. NO! MN, Ngl ooo N ea 


Vamos introduzir agora duas hipóteses adicionais que podem ser deduzidas 
rigorosamente no contexto da termodinâmica, mas que são suficientemente 
intuitivas para serem aceitas sem demonstração: 


d) Quando o sistema está em equilíbrio, a distribuição de partículas entre 
os diferentes estados é a mais provável. 

e) A distribuição mais provável de partículas é a que pode ser obtida maior 
número de vezes, isto é, é a que maximiza o valor de W na Eq. 7.1. 


Estamos agora de posse de todas as informações necessárias para che- 
garmos à chamada distribuição de Boltzmann. Essa distribuição também é 
chamada de distribuição clássica, pois ignora o Princípio de Exclusão de Pauli. 

Vamos recordar as hipóteses que conduzem à distribuição de Boltzmann: 


a) as partículas são distintas; 
b) o número total de particulas é conservado: 
c) a energia total é conservada; 
d) em equilíbrio, a distribuição de particulas é a mais provável; 
e) a distribuição mais provável é a que maximiza W. 
Nosso problema é portanto maximizar a função W (Eq. 7.1) usando como 
condições adicionais 
NAN AN... +... + N,=EN,= Número total de particulas = 


= Constante; (Eq. 7.2) 
NE, + N,E, + N,E, +... = LEN,E, = Energia total = 
= Constante, (Eq. 7.3) 


onde E, é a energia do estado i. 

Um artifício bastante conveniente para facilitar o desenvolvimento mate- 
mático é maximizar InW* em lugar de W, o que evidentemente não altera a 
solução do problema. 

A condição necessária para que InW seja máximo é 


dnW = DOM an, =0. (Eq. 7.4) 
i i 
Além disso de acordo com as Egs. 7.2 e 7.3, devemos ter 
LN; = constante .. EdN, — 0; (Eq. 7.5) 
1 | 
3 N,E, = constante v. 5 EdN,=0. (Eq. 7.6) 
i f 


1 
*InW = Logaritmo neperiano de W 
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Neste ponto temos que fazer alguns comentários sobre o problema ma- 
temático que enfrentamos. Suponhamos que queremos maximizar uma função 
Fx y; 2), 

Ora, sabemos que 

df = dx + dy + a 

Se os incrementos dx, dy e dz são dias a condição necessária 

para que f seja máxima é 
q o af cai 
dx O Jd 

No nosso caso, queremos maximizar lnW, mas os incrementos dN, não são 
independentes, pois devem satisfazer às Eqs. 76 e 7.7. É óbvio portanto que 
não podemos anular diretamente as derivadas de InW em relação às variáveis 
N;,, como fizemos para encontrar o máximo da função f(x, y, z). 

Torna-se então necessário utilizar um artifício, os chamados coeficientes 
indeterminados de Lagrange. A idéia é transformar o problema para o caso mais 
simples em que os incrementos são independentes entre si. 

Vamos multiplicar as Egs. 7.5 e 7.6 por « e f, respectivamente, e somar 
o resultado à Eq. 7.5. Fazendo isso, obtemos 


X[m +u+ p£.| dN, = 0. (Eg. 7.7) 
i 


0N, 

Ora, podemos sempre escolher dois dos incrementos, dN, e dN,, por 
exemplo, de tal forma que as Egs. 7.5 e 7.6 sejam satisfeitas, qualquer que seja 
o valor dos outros incrementos. Em seguida, escolhemos os coeficientes q e B 
de forma que os coeficientes de dN, e dN, na Eq. 7.7 se anulem, isto é, fazemos 

OInW 
0N, 
InW 
oN, 

Ora, nesse caso, a Eq. 7.7 se reduz a 


z Es + 0 + pE, am, =, (Eq. 7.8) 
iél,2 


+u+ BE, =0, 


+24 + BE, =0. 


em que os incrementos dN,, dN,,... são todos independentes. Em conse- 
qiiência, podemos igualar separadamente a zero cada uma das parcelas do 
somatório, obtendo'*! 


OlnW 
Tp (Eq. 79) 
Mas de acordo com a Eq. 7.1, 
InW [à ê 
om = o MN! = ln] = = say SON = og daN (50.710) 


iWEstritamente falando, a Eg. 7.9 só é válida para i = 3,4,5,...; entretanto assim como esco- 
lhemos os incrementos dN, e dN,, poderíamos ter escolhido qualquer outro par de incrementos; 
assim a equação pode ser aplicada, qualquer que seja o valor de i 
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Para qualquer sistema de dimensões práticas, N, é um número extrema- 
mente grande; em consequência, podemos usar a aproximação de Stirling para 
a função fatorial: 


- InN,! = NAnN,-1). (Eq. 7.11) 
Então 


= nN,. (Eq. 7.12) 


De acordo com as Egs. 7.9 e 7.12, temos 


InN,=0 + BE,; 
N = et BE; Es ee Es, 


Fazendo & = A, temos finalmente 
N, = AH, (Eq. 7.13) 


que é a chamada distribuição de Boltzmann. Repare que na distribuição de 
Boltzmann o número de partículas N, com energia E, depende exponencialmente 
da energia E, do estado considerado. 

Em geral, a Eq. 7.13 é escrita sem o indice i, isto é, na forma 


N = Ag, (Eq. 7.14) 


Uma dedução puramente matemática, como a apresentada anteriormente, 
não permite calcular o valor de $, que deve ser obtido através de considerações 
fisicas. É possível provar, entretanto, utilizando conceitos da termodinâmica, 
que 4 =—1/kT, onde k é a constante de Boltzmann e T a temperatura absoluta. 

Vamos agora aplicar a Eg. 7.14 a um problema específico; o cálculo da 
distribuição de energia potencial das moléculas de um gás perfeito. 

Chamando de z a cota do ponto onde está localizada a molécula, sua energia 
potencial é dada por 

E = mgz, (Eg. 7.15) 


onde m é a massa da molécula e g a aceleração da gravidade. 
De acordo com a Eg. 7.14, o número de moléculas com cota z (e conse- 
quentemente com energia mgz) é dada por 


N, = der, (Eq. 7.16) 


A Eq. 7.16 está representada graficamente na Fig. 7.3. 
N 


Figura 7.3 Distribuição de Boltzmann 
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Vamos calcular o valor da constante 4. Chamando de N,. o número total 
de molécula do gás, temos 


=x su 


Nr = |n.de - | e menr — qkT 
“rt ez mg 
e portanto 
N 
dl ie = (Eq. 7.17) 


De acordo com as Egs. 7.16 e 7.17, a distribuição de energia potencial das 
moléculas de um gás perfeito é dada por 
— Namg 


N.o= 


=p e“, (Eq. 7.18) 


7.2. SISTEMAS CONCENTRADOS — DISTRIBUIÇÃO 
DE FERMI-DIRAC 


No parágrafo anterior, estudamos o caso de sistemas diluídos de partículas 
distintas; consideremos agora o caso de um sistema constituído por partículas 
idênticas; além disso, suponhamos que o sistema é concentrado, isto é, que há 
necessidade de aplicar o Princípio de Exclusão de Pauli. A distribuição que 
resulta é a chamada distribuição de Fermi-Dirac'*. 

As hipóteses que conduzem à distribuição de Fermi-Dirac são as seguintes: 
a) as partículas são idênticas; 

b) o número total de particulas é conservado; 

c) a energia total é conservada; 

d) o Princípio de Exclusão de Pauli é respeitado; 

e) a distribuição mais provável é a que maximiza W 

Nosso primeiro problema é calcular o valor de Wí Seja N, o número de 
estados permitidos com energia E,. e n, o número de partículas distribuídas 
entre estes N, estados (Fig. 7.4). 

E IR , estados 
* |n, partículas 


a » estados 


E ? 
2 ln, particulas 


E N ;, estados 
n; partículas 


Figura 74 Distribuição das n, partículas entre os estados permissíveis 


(Como a distribuição de Fermi-Dirac leva em conta o Princípio de Exclusão dc Pauli, que 
não tem análogo algum na Mecânica Clássica, também é chamada de distribuição quântica 
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Consideremos o grupo E, ; como em cada um dos estados N, só pode haver 
uma das particulas n/*!, o número total de modos de colocar as n, partículas 
nos N, estados é dado por 


N;! 


W=Cij= n; (N,-nj! 


(Eq. 7.19) 


O problema agora é maximizar a função W usando como condições adi- 
cionais: 

En, = constante; (Eg. 7.20) 

En,E, = constante. (Eg. 7.21) 


Deixamos a cargo do leitor o desenvolvimento completo usando os coe- 
ficicntes de Lagrange (Problema 7.1). O resultado é o seguinte: 


N; =) 2 
np = [+ dh (Eg. 7.22) 


Fazendo f = No vemos que f corresponde à probabilidade de encon- 


1 


trarmos n, partículas com energia E, (para a qual existem N, estados). 
De acordo com a Eq. 7.22, temos 


(= w Ire 
Novamente é possivel mostrar que fi = — EP obtendo-se então 
7a l : 
J a ] dh Et EukT 
Definimos uma nova variável E, tal que 
E, 
PER 


Nesse caso, temos 


l 


f o 1 + gfE ErdlkT ” 


ou eliminando o indice i: 


l 


fa T+ E EnAT (E g. 7.25) 


A Fig. 7.5 mostra o gráfico de / em função da energia para várias tem- 
peraturas. 

Repare que para T =0 todos os estados com E < E, estão ocupados 
(J = 1) e todos os estados com E > E, estão vazios (/ = 0). Repare ainda que 


*ilsso é uma consegiência do Principio de Exclusão de Pauli 
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Figura 7.5 Distribuição de Fermi- 
-Dirac 


em qualquer temperatura, / = 1/2 para E = E,, isto é, que a probabilidade 
de ocupação do nível E, é sempre igual a 1/2. 

A energia E, é chamada de nível de Fermi. 

É fácil demonstrar (Prob. 7.2) que para sistemas diluídos, isto é, quando 
n, É N,, a distribuição de Fermi-Dirac se reduz à distribuição de Boltzmann. 


7.3. OUTRAS DISTRIBUIÇÕES 


A distribuição de Boltzmann e a distribuição de Fermi-Dirac não são as 
únicas distribuições possíveis. Dependendo das hipóteses iniciais, outras dis- 
tribuições podem ser obtidas. Por exemplo: se supusermos que as particulas 
são idênticas, mas não obedecem ao Princípio de Exclusão de Pauli (como é 
o caso dos fótons) chegaremos à chamada distribuição de Bose-Einstein. O 
leitor interessado em maiores pormenores deve consultar um livro de Mecânica 
Estatistica (Kittel, 1958). 


EXERCÍCIOS 


7.1. Demonstre a Eg. 7.22. 

7.2. Mostre que para sistemas diluídos a distribuição de Fermi-Dirac se reduz à dis- 
tribuição de Boltzmann. 

7,3. Determine a distribuição de energia cinética das moléculas de um gás perfeito 

(distribuição de Maxwell). Para calcular a constante multiplicativa, suponha que 
o gás contém N moléculas. 

7.4. Neste exercicio vamos mostrar com um exemplo particular que a constante 5 da 
Eq. 7.14 é igual a —1/kT. Considere um gás perfeito em um recipiente a uma tem- 
peratura T sob a ação da gravidade. Definindo: 


n(z) — número de moléculas com cota z por unidade de volume, 


m — massa de cada molécula; 
g — aceleração da gravidade; 
p — massa especifica do gás; 
k — constante de Boltzmann; 


T — temperatura absoluta. 
a) Prove, sem usar a distribuição de Boltzmann, que 
n(z) = Ae” msitT 


onde 4 é uma constante. 
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b) Compare essa expressão com a Eq. 7.14 e verifique que 
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8 Aplicação da distribuição de 
Fermi-Dirac a um material 
semicondutor 


8.1. APROXIMAÇÕES PARA A FUNÇÃO DE DISTRIBUIÇÃO 
DE FERMI-DIRAC 


A distribuição de energia dos elétrons em um sólido é governada pelas 
leis de estatística de Fermi-Dirac, já examinadas no capitulo anterior. O prin- 
cipal resultado dessa estatística é a função de distribuição de Fermi-Dirac, que 
nos dá a probabilidade de que um estado eletrônico com energia E esteja ocupado 
por um elétron: 


: l 
J(E) = Tx E EMT (Eg. 8.1) 
Essa função contém um parâmetro, E,, que é chamado de nível de Fermi. 
Uma definição rigorosa do nível de Fermi o descreve como o potencial químico 
dos elétrons em um sólido. Para nossos propósitos, entretanto, é suficiente 
observar que o nível de Fermi é a energia para a qual a probabilidade de ocupação 
desse nível de energia é exatamente 1/2(/ = 1/2 se E = E,). 

A função de distribuição de Fermi-Dirac está representada na Fig. 8.1(a) 
para o caso de um semicondutor intrínseco. Na faixa de condução existe um 
grande número de estados disponíveis. Entretanto a probabilidade de ocupação 
desses estados é pequena; assim, haverá um número relativamente pequeno de 
elétrons na faixa de condução. A faixa de valência também tem um grande nú- 
mero de estados disponíveis. Quase todos esses estados estão ocupados por 
elétrons, já que a probabilidade de ocupação desses estados é grande; assim, 
haverá um número relativamente pequeno de estados desocupados, isto é, de 
buracos, na faixa de valência. 

A função de distribuição de Fermi-Dirac é simétrica com relação ao nível 
de Fermi E,.. Assim se o número de estados disponíveis nas faixas de condução 
e de valência é o mesmo e se o número de elétrons na faixa de condução é igual 
ao número de buracos na faixa de valência (o que sabemos ser verdade no caso 


E E E 
sa Er -ses0ee Ec E 
Er 
Er 

Er 

ip p=y —— Ey —svooe- Er 
FE) fg) FE) 
(a) (b) (c) 


Figura 8.1 Aplicação da distribuição de Fermi-Dirac a semicondutores: intrínseco (a), 
tipo N (b) c tipo P (c) 
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de um semicondutor intrinseco), o nível de Fermi estará no centro da faixa 
proibida. O nível de Fermi de um semicondutor intrinseco é chamado de nível 
de Fermi intrínseco e designado pelo simbolo E,. 

Em um semicondutor tipo N, o número de elétrons na faixa de condução 
é maior que no caso intrínseco. Como, entretanto, a densidade de estados dis- 
poniveis de energia é o mesmo que no caso intrinseco, conclui-se que, em um 
semicondutor tipo N, o nivel de Fermi, e com ele toda a função de distribuição 
de Fermi-Dirac, será deslocado para cima, para mais perto da faixa de con- 
dução. Ao contrário, em um semicondutor tipo P o nível de Fermi e a função de 
distribuição de Fermi-Dirac se deslocam para baixo. Os dois casos estão ilus- 
trados nas Figs. 8.1(b) e 8.1(c). 

Para energias muito acima ou muito abaixo do nivel de Fermi, a distri- 

buição de Fermi-Dirac pode ser aproximada pelas seguintes expressões: 


f(E) = et ErIRT (para |E-E,|>kT E> E,), (Eq. 8.2) 
S(E)= 1-e tt EMT (para |[E-E|>kT E<E,). (Eq. 8.3) 


O segundo termo do segundo membro da Eg. 8.3 pode ser considerado 
como a probabilidade de ocupação de um nível de energia E por um buraco. 


8.2. CONCENTRAÇÃO DE PORTADORES EM UM CRISTAL 


Normalmente, a concentração de impurezas em um semicondutor tipo N 
ou tipo P não é suficiente para que o nível de Fermi se aproxime muito da faixa 
de condução ou de valência. Nesse caso as aproximações anteriores são válidas 
e é fácil calcular a concentração dos portadores em função do nível de Fermi. 


8.2.1. Determinação das concentrações de elétrons e buracos 
em função do nível de Fermi 


A probabilidade de que um estado da faixa de condução esteja ocupado, 
no caso em que E-E, > kT. é dada por 


JE) = e = Eh, (Eq. 8.4) 


Por outro lado, a densidade de estados por unidade de volume na faixa de 
condução é dada pela Eg. 3.35 


1 /2m*? 1 
DJE) = à (5) [E- Ec] 2, (Eq. 8.5) 
onde tomamos como referência de energia o topo da faixa de valência (Fig. 8.2). 
Repare que a Eq. 3.35 foi deduzida para o Gás de Fermi. O que fizemos 
foi supor que o efeito do cristal sobre os portadores pode ser descrito por uma 
massa efetiva, m,, usada em lugar da massa “real”, m (veja o parágrafo 6.6). 
Ora, conhecendo a probabilidade de ocupação dos estados da faixa de 
condução, /(E), e a densidade de estados disponiveis nessa mesma [aixa, D(E). 
podemos calcular a concentração de portadores, n, através da equação: 


n= | D(E) f(E)dE. (Eq. 8.6) 


Es 
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Faixa de 


Condução 
E=EG+4E 


Refêrencia 


Faixa de 


Valência 


Figura 82 Diagrama ilustrando a referência de ener- 
gias (topo da faixa de valência) 


Substituindo as Eqgs. 84 e 8.5 em 8.6, obtemos 


1 (2m'? ki 
n=5 ( a) pree [ [E-Eç]'2e FNT dE, (Eq. 8.7) 
Es 
Fazendo E — E; = x, temos 
| (2m O er-eonkT - 1/2, 
=— E id x Rd. Eq. 8.8 
"E 2a ( hê ) E p a E ARS) 


A integral anterior pode ser encontrada em qualquer tabela de integrais. 
O resultado é 


| Hg dy = 3 mkT. (Eg. 8.9) 
Substituindo a Eq. 89 em 8.8 e lembrando que h = h/2r, obtemos 
3/2 
n= (Ee) Pa (Eq. 8.10) 


A Eg. 8.10 pode ser usada para calcular a concentração de elétrons a partir 
dos parâmetros do material considerado. A Eg. 8.10 em geral é escrita na forma 


p= NB SE, (Eq. 8.11) 


9 3/2 
N,= (E) (Eq. 8.12) 


onde 


é uma constante chamada de densidade efetiva de estados na faixa de condução. 

Para calcular a concentração de buracos, basta observar que a probabi- 
lidade de ocupação dos níveis da faixa de valência por buracos é dada pelo 
segundo termo da Eq. 8.3, que pode ser escrito na forma: 


SAE) = ef FAT (para [E-Ep| > kT) (Eg. 8.13) 
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e que a densidade de estados por unidade de volume na faixa de valência é dada 
por 
1 Im N32 
DXE)=—=(=*| (CE), Eq. 8.14 
“= a (68) CRP (Eq. 8.14) 
onde m, é a massa efetiva dos buracos. 
Deixamos a cargo do leitor o restante da manipulação matemática para 
provar que 


de p=N,e "Mr, (Eg. 8.15) 
on "3/2 
N,=2 (a) (Eq. 8.16) 


é chamada de densidade efetiva de estados na faixa de valência. 

Por outro lado, se tomarmos como referência um nível arbitrário de energia, 
localizado abaixo do topo da faixa de valência (Fig. 8.3), as Egs. 8.11 e 8.15 po- 
derão ser escritas na forma mais usual 


n= N er EMT, (Eq. 8.17) 
p=N,e EM, (Eq. 8.18) 
onde E. e E, são as energias correspondentes à base da faixa de condução e ao 
topo da faixa de valência, respectivamente. Repare que a escolha do nível de 


referência é totalmente arbitrária, já que as Egs. 8.17 e 8.18 envolvem apenas 
diferenças entre os diversos níveis. 


Faixa de 
Condução 


Nivel de | 


) iZ 
Figura 8.3 Referência arbitrária para a energia Earl 


Faixa de 
Valência 


8.2.2. Produto np 
Multiplicando a Eg. 8.11 pela Eg. 8.15, obtemos 
np=NN,e Fe, (Eg. 8.19) 
A Eq. 8.19 mostra que em qualquer material semicondutor, o produto da 
concentração de eletrons pela concentração de buracos não depende das impurezas 
presentes. 
8.2.3. Concentrações de portadores em um semicondutor intrínseco 


Em um semicondutor intrinseco. a concentração de elétrons é igual à con- 
centração de buracos, isto é 
n= Pp = hs (Eq. 8.20) 


onde n, é chamada de concentração intrínseca de portadores. 
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Combinando a Eq. 8.20 com a Eq. 8.19, obtemos 

ni =N,N é fot”, (Eg. 8.21) 

ho=/ N;N e Patr, (Eg. 8.22) 


A Eg. 8.22 permite calcular a concentração intrínseca de portadores, co- 
nhecidos os parâmetros do material e a temperatura em que este se encontra. 


ou 


8.2.4. Taxa de geração térmica 


Podemos também usar as Egs. 8.10, 8.15 e 8.16 para calcular a taxa de 
geração térmica de um semicondutor em equilíbrio térmico. Com efeito, de acordo 
com a Eq. 6.2, temos 


GT) = qnp. (Eq. 8.23) 
Substituindo n e p por seus valores, temos 
GT) = AT?" FaltT, (Eg. 8.24) 
onde 
mk 
4=3R2 (5) (mm)? (Eq. 8.25) 


é uma constante independente da temperatura. 

A Eq. 8.24 mostra que a variação da taxa de geração térmica com a tem- 
peratura é da forma T*e”“7, Essa dependência pode ser verificada experimen- 
talmente. 


8.2.5. Nivel de Fermi de um semicondutor intrínseco 


Como vimos, as seguintes expressões são válidas para qualquer semicon- 
dutor: 
n= N  dtBr EoRT 
p= N,e7t8r— EVT. 
Para um semicondutor intrínseco, n =p =n,. Assim chamando o nivel 
de Fermi do semicondutor intrínseco de E,, temos 


(E EST — —(Ei— Ey)/kT 
Ne dE =N,e , 


u k7 N, 


Ora, de acordo com as Egs. 8.12 e 8.16, a única diferença entre N e N, 
está no fato de que as massas efetivas dos elétrons e dos buracos não são iguais. 
Como em geral a diferença entre as duas massas não é grande, o nível de Fermi 
do material intrínseco está localizado aproximadamente no centro da faixa 
proibida. Com efeito, se tomarmos N, = N,., In(N,/N)) = 0 e portanto 
Ee à Eos (Eq. 8.27) 


E, = 
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8.2.6. Concentrações de portadores em função de n, e E, 


As concentrações de elétrons e buracos em um semicondutor tipo P ou 
tipo N também podem ser expressas em função da concentração intrínseca, n,, 
e do nível de Fermi intrínseco, E,. Com efeito, combinando as Egs. 8.17, 8.18, 
8.22 e 8.26, é fácil demonstrar que 


n= nelbr- EMT, (Eq. 8.28) 
p= en TRE, (Eg. 8.29) 


Naturalmente, quando E, = E, nas equações anteriores, obtemos n =p =n,, 
a concentração de portadores em um material intrinseco. 


8.3. Neutralidade de carga elétrica nos semicondutores 


Consideremos um bloco isolado de um semicondutor. É evidente que o 
bloco como um todo deve ser neutro, isto é o número de cargas negativas no 
interior do bloco deve ser igual ao número de cargas positivas. Essa neutralidade 
global é chamada de neutralidade macroscópica. 

A neutralidade macroscópica não impede que existam pequenas regiões no 
cristal em que a neutralidade local (a chamada neutralidade microscópica) não 
é respeitada. A Fig. 8.4 mostra um exemplo. Se considerarmos a pequena região 
em tomo do átomo de arsênio, encontraremos um saldo de carga positiva. 
Entretanto, o elétron cedido à faixa de condução pelo átomo de arsênio não 
pode sair do cristal, de modo que o cristal como um todo continua neutro. 


Cristal de Si | || ue Elétron livre 
24 


06 
si 


Vejamos agora as implicações da neutralidade macroscópica para o equi- 
librio dos portadores em um semicondutor. 


Figura 8.4 Neutralidade macroscó- 
pica. O cristal como um todo con- 
tinua neutro 


8.3.1. Semicondutores intrinsecos 


No caso dos semicondutores intrínsecos, a condição de neutralidade pode 
ser expressa pela relação n = p, que é satisfeita automaticamente devido à na- 
tureza dos processos de geração e recombinação de portadores!” 


8.3.2. Semicondutores extrinsecos 


No caso dos semicondutores extrinsecos (tipo N ou tipo P), a condição 
de neutralidade se torna mais complexa, pois além dos portadores, os átomos 
ionizados das impurezas também participam do equilibrio de cargas. 

Suponhamos um semicondutor que possui uma concentração N, de impu- 
rezas aceitadoras e uma concentração N , de impurezas doadoras. Essas impurezas 


*Mais rigorosamente, deveriamos escrever qn = qp para a condição de neutralidade, mas 
como as cargas dos portadores são iguais em valor absoluto, clas podem ser ignoradas 
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Figura 8.5 Diagrama ilustrando os niveis de energia E,., 
Ep, Es. Ep 
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introduzem, na faixa proibida, os niveis de energia E, e E,, respectivamente 
(Fig. 8.5). 

No caso mais geral, nem todos os átomos doadores e aceitadores estarão 
ionizados. Vamos chamar de N a concentração de átomos doadores ionizados!* 
e de N$ a concentração de átomos doadores neutros. Podemos imediatamente 
escrever 

N$ = No-N$:. (Eg. 8.30) 
Do mesmo modo, 


N7=N,-Nº, (Eq. 8.31) 


onde Nj é a concentração de átomos aceitadores ionizados** e N4 a con- 
centração de átomos aceitadores neutros. 
Nesse caso a condição de neutralidade macroscópica assume a seguinte 


forma: 
n+ Na =p + No. (Eg. 8.32) 


84. CONCENTRAÇÕES DE PORTADORES EM UM 
SEMICONDUTOR EXTRINSECO 


As Egs. 8.17 e 8.18 só permitem calcular as concentrações de portadores 
se conhecermos a posição do nível de Fermi, E, , em relação ao topo da faixa de 
valência ou em relação à base da faixa de condução. No caso de um semicondutor 
intrinseco, como vimos, a posição do nível de Fermi pode ser calculada facil- 
mente (Eg. 8.26), a partir da condição n = p. 

No caso de um semicondutor extrinseco, em que n & p, a situação é mais 
complexa. Para determinarmos o nível de Fermi, temos que usar a equação 
de neutralidade (Eq. 8.32). Para isso, entretanto, é preciso exprimir N4 e N5 
em função do nível de Fermi, É o que veremos a seguir. 


8.5. CÁLCULO DAS CONCENTRAÇÕES DE IMPUREZAS 
IONIZADAS EM FUNÇÃO DO NÍVEL DE FERMI 


Vamos começar estudando o caso das impurezas aceitadoras. Podemos 
escrever 


Ná = SuN (Eq. 8.33) 


“O sinal positivo indica que os doadores ionizados são ions positivos 
!**0 sinal negativo indica que os aceitadores ionizados são ions negativos 
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onde f, é a probabilidade de que o estado de energia E, introduzido pela impu- 
reza aceitadora na faixa proibida esteja ocupado (ou seja, a probabilidade de 
que o átomo da impureza esteja ionizado). 

Ora, de acordo com a Eq. 8.1, a probabilidade de que um estado de energia 
E, esteja ocupado por um elétron é dada por 


| 

fi= 1x EA EE (Eq. 8.34) 
De acordo com as Egs. 8.33 e 8.34, 
Ns 


NZ (Eq. 8.35) 


-. 1 + FA ENAT” 


equação que permite calcular a concentração de átomos aceitadores ionizados 
em função do nível de Fermi. 

No caso das impurezas doadoras, temos 

NS = SoNo: (Eq. 8.36) 

onde f, é a probabilidade de que o estado de energia E, introduzido pela impu- 
reza doadora na faixa proibida esteja ocupado (ou seja, a probabilidade de que 
o átomo da impureza não esteja ionizado). 

Nesse caso, portanto, 


Ng =(1-fo)No. (Eq. 8.37) 
O valor de f, é dado por uma equação análoga à Eq. 8.34: 
1 
Jo (Eq. 8.38) 


— L+ g!Ep= Er)lkT 


De acordo com as Egs. 8.37 e 8.38, 
pia — Ma 
1+ etEr- Ep)ikT 
equação que permite calcular a concentração de átomos doadores ionizados 
em função do nível de Fermi. 
Substituindo n.p. N4 e Nj por seus valores na Eq. 8.32, obtemos 


(EB N 
Ne (Ee ErykT 4 Tr Er Ear” (Eg. 8.40) 


(Eg. 8.39) 


4 Es — (Er EJkT 
1 + eta” EpdkT N,e ç + 


Teoricamente a Eq. 8.40 permite calcular o nivel de Fermi a uma dada 


Lemnperatuta, st CONDCUCILVS US paraiicuvo Lrjy Irys 27 4s 273 Ds DA» Dry Ep 
Na prática, entretanto. a Eq. 8.40 não possui solução analítica, de modo que é 
preciso recorrer a métodos numéricos ou gráficos. Um método gráfico que 
mostra com clareza a influência dos diferentes parâmetros sobre a posição do 
nível de Fermi e a concentração dos portadores será examinado a seguir. 


8.6. MÉTODO DE SCHOCKLEY PARA A DETERMINAÇÃO DAS 
CONCENTRAÇÕES DE PORTADORES EM UM 
SEMICONDUTOR INTRÍNSECO 
Voltemos à Eq. 8.17, que exprime a concentração de elétrons em função 
do nível de Fermi: 


p= NA ET. (Eg. 8.41) 
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Tomando os logaritmos de ambos os membros da Eq. 8.41, obtemos 


— Ep 


E 
ga qr + a 
SEM (Eg. 8.42) 


log,,n = log, N, 


Assim, se traçarmos um gráfico de n em função de E,, usando uma escala 
logaritmica para n e uma escala linear para E, , obteremos uma linha reta (Fig. 
8.6). 


Figura 8.6 Variação do logaritmo 
da concentração de elétrons com o 
nivel de Fermi 


A Eg. 8.42 mostra que uma variação de uma década na concentração de 
elétrons, n, corresponde a uma variação de 2,3kT no nível de Fermi, Ep. 
Analogamente, a variação da concentração de buracos com o nivel de 
Fermi pode ser expressa pela equação 
E,- E, k 
kT 


A Eq. 8.43 está representada graficamente na Fig. 8.7, 


log,oP = l0g,9N, + (Eq. 8.43) 


Figura 8.7 Variação do logaritmo 
da concentração de buracos com o 
nível de Fermi 


Por outro lado a concentração de doadores ionizados em função do nivel 
de Fermi é dada pela Eq. 8.39: 


ps 1 
No = Noyy q (Ep ERRT (Eq 844) 


Vamos analisar por partes a variação de Np com Ep. Quando E, = E,, 
N; = Np/2. Quando E, < E,. (E,-Ep)kT >0 e portanto N$ = N,. Por 
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outro lado, quando E, > E,,e Fo EMkT > 1 e portanto 


N 
Fm D = (Er EplkT 
No * o 1Ep- EPT O Nope 


log o Nj = log,oNp + Fome. (Eq. 8.45) 
A Eq. 8.45 tem a mesma forma que a Eq. 8.43, o que mostra que quando 
E, > Ep, N5 varia com E, da mesma forma que p. 
Os casos extremos E, < E, e E, > Ep nos dão as assintotas da curva 
que representa a variação de log,oNp com E,. Para obtermos alguns pontos 
intermediários, podemos usar a seguinte tabela, baseada na Eg. 8.44: 


E, = E, fa tn e (En EMT Nó N 

11 kT Li 0,33 Np/1,33 = 0,75 Np 

0,1 kT -0,1 0,90 Np/1,90 = 0,53 Np 

O kT 0 1 Np/2 =0,50N, 
-0,1 kT +0,1 1,10 No/2,1 =048Np 
—1,1 kT +1I0 3,00 No/4 =025N, 
-1,39kT + 1,39 4,00 No/5 =02N, 
1,95 kT +1,95 7,00 No/8 =0,125Np 
-23 kT +2,30 -9,00 Ny/10 =01N, 


Podemos agora traçar o gráfico da variação de log Nj com Ep. (Fig. 8.8). 


log ND 


| Assintotas 


== Nó | 


+ 


23kT/ década 


| Er=ED 
+ 
ND=ND/2 


ED 
Eb>> Er 
ND=ND Her >>Ep 


Figura 8.8 Variação de log N$ com Ep. Para E, > Ep, essa variação é linear, com coe- 
ficiente angular de uma década de variação de Ny para cada 2,3kT de variação de Ep 
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Na prática, três pontos na região intermediária são suficientes para se 
conseguir uma precisão razoável no traçado da curva. 

A análise do comportamento de N, em função de E, é inteiramente aná- 
loga e fica a cargo do leitor (Prob. 8.9). 

Agora que descrevemos a variação de n,p, N$ e N4 com o nível de Fermi, 
vejamos em que consiste o método de Schockley. A idéia é traçar em um mesmo 
gráfico (Moll, 1964) as curvas de n,p, N5 e N4 em função de E, e resolver a 
equação n + N4 =p + N5, procedendo da seguinte forma (Fig. 8.9): 


1. Traçam-se as curvas de n,p,N; e Np; 

2. Somam-se ponto a ponto os valores de ne N4; 

3. Somam-se ponto a ponto os valores de p e Np; 

4. Determina-se o ponto de interseção das curvas log(n + N;) e log(p + Np). 


log (p.n,NB,NA) log (p,n,NG.NA ) 


Figura 8.9 Ilustração do método de Schockley para a determinação de E,. No ponto 
de interseção, p + Nj =n + Ná 


A abscissa do ponto de interseção é o valor de E, que satisfaz à Eq. 8.40. 

Os valores de n, p, N5 e N4 também podem ser obtidos diretamente do gráfico. 

Na prática, entretanto, é mais conveniente obter do gráfico apenas E, (como 

se trata de uma escala linear, o erro não é muito grande) e calcular n e p através 
das equações 

n= Ne” (E ErlikT, (Ea. 8.46) 


pot, (Eq. 8.47) 
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ou então através das equações 


p= Ne rt, (Eq. 8.48) 
2 
n= Es (Eq. 8.49) 


Do mesmo modo, é mais conveniente calcular Nj e N ; através das equações 


| 
N$ = Nous (Eq. 8.50) 
a 1 
Na =NTr o EE (Eg. 8.51) 


Com referência à Fig. 8.8 convém notar o seguinte: 


1) A escala das ordenadas é logaritmica. Assim para E, = E,, p+ Nj =p; 
para Ep = E.,n+N, =»n. Na interseção da curva de p com a curva de 
No, Pp + Ná = 2N5 (em escala logaritmica). 

2) A inclinação das curvas de n c p e das assíntotas oblíquas de Nj e N; de- 
pende da temperatura. 


EXERCÍCIOS 


88.1, Considere um semicondutor tipo N. em equilíbrio térmico. que contém uma con- 
centração Ny de átomos doadores. Chamando de n, e p, as concentrações de 
portadores em maioria e minoria no material, respectivamente: a) calcule n, e p,, 
supondo que todas as impurezas estão ionizadas: b) a aproximação p, = n;/N, é 
muito usada na prática, Para um cristal de germânio com N, = 7,5 x 10 em”?, 
calcule o erro cometido usando essa aproximação. 

8.2. Calcule a localização do nível de Fermi intrínseco no silício a -78ºC, 27"C e 
300ºC. É razoável supor que o nivel de Fermi está no centro da faixa proibida? 

8.3. Desenhe (em escala) o diagrama de niveis de energia para o silício com 10! átomos 
de boro por centimetro cúbico. Indique as faixas de condução e de valência, o 
nivel de Fermi intrinseco e o nivel de Fermi real. 

8.4. O ouro no silício é uma impureza peculiar, pois pode comportar-se tanto como 
doador como aceitador. Em outras palavras, o ouro introduz dois níveis na faixa 
proibida, um a 0,35eV da banda de valência e outro a 0,54eV da banda de con- 
dução. Pergunta-se: a) Um cristal de silício com uma certa concentração N, de 
ouro é tipo N ou tipo P?; b) Qual é a posição do nivel de Fermi? 

8.5. Quando discutimos o fenômeno da neutralidade macroscópica, consideramos 
apenas as impurezas mais simples. Existem impurezas, como o zinco no silicio, 
que podem se comportar como um aceitador simples ou como um aceitador duplo, 
isto é, que podem capturar um ou dois elétrons. No caso de um cristal de silício 
com uma concentração N, de átomos de zinco e N, átomos de fósforo, escreva 
a equação de neutralidade macroscópica. 

!*8.6. O zinco introduz no silício os seguintes níveis na faixa proibida; E, , 0,35 eV acima 
do topo da banda de valência; E,. 0,2eV acima de E,. Para um cristal de silício 
com 2 x 10!º átomos/cm” de As e 0,9 x 10'* átomos/cm” de Zn, calcule; a) a 
posição do nível de Fermi; b) a concentração de elétrons; cj a concentração de 
buracos. 
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8.7. a) Explique o fenômeno da compensação de impurezas, isto é a diminuição do 
número de portadores e da condutividade de um cristal quando este contém si- 
multaneamente impurezas tipo P e tipo N; b) Podemos verificar se um cristal é 
intrinseco simplesmente medindo a sua condutividade? Justifique. 

8.8. Calcule n,p e E, para um eristal de silício contendo 1,1 x 10'º átomos/em” de 
boro c 9 x 10'* átomos/cm* de fósforo a 27ºC. 

8.9. Analise a variação da concentração de aceitadores ionizados com o nivel de Fermi. 
Desenhe a curva de log,, N, em função de E,. 
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9 Transporte de cargas nos 
semicondutores 


Neste capitulo vamos estudar o movimento dos portadores de carga nos 
semicondutores e apresentar as leis que regem os dois processos básicos de 
transporte de cargas, que são a deriva” e a difusão. 


9.1. ESTUDO DA DERIVA 


9.1.1. O processo de deriva 


Consideremos um semicondutor tipo N com uma concentração uniforme 
de doadores. Na ausência de campo elétrico externo, os elétrons do cristal estão 
em movimento térmico aleatório, interrompido por colisões (Fig. 9.1a). 


——————— 
tg ts 
tg y go O & 
to 
to 
b A . HR dá b N 
total; totsta ta taty to Deslocamento 
ts t totalemt=to 
(a) (b) E (c) 


Figura 9.1 Deriva de portadores. (a) Movimento térmico dos portadores na ausência 
de campo elétrico externo. (b) Componente de deriva, produzida pelo campo externo. 
(c) Eleito combinado da deriva e da agitação térmica 


É fácil demonstrar que o deslocamento total do elétron, depois de um 
grande número de colisões, é zero, pois não há uma direção preferencial para 
o movimento do elétron depois de cada colisão. 

Se aplicarmos um campo elétrico ao cristal, o elétron será acelerado pelo 
campo elétrico entre cada duas colisões. Essa aceleração tem uma direção pre- 
ferencial: a do campo elétrico aplicado. Assim, superposta à velocidade tér- 
mica, haverá uma pequena componente da velocidade devido ao efeito do campo 
elétrico [Fig. 9.1(b)]. Essa componente é chamada de velocidade de deriva; se 
somarmos os efeitos da velocidade térmica aos da velocidade de deriva, teremos 
o movimento real do elétron sob a ação do campo elétrico [Fig. 9.1(c)]. Para 
termos uma idéia de ordens de grandeza, a velocidade térmica dos elétrons no 
silício à temperatura de 27ºC é de 10º m/s, enquanto a velocidade de deriva 
sob o efeito de um campo elétrico de 10 V/cm é de 390 m/s, muito menor, 
portanto. 

Como a velocidade de deriva está sempre superposta a uma velocidade 
térmica muito maior, não faz sentido acompanharmos o movimento de um 


'“Em inglês, drift 
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elétron em particular; o que devemos fazer é calcular o movimento médio de 
muitos elétrons em um momento específico de tempo. 

Suponhamos que a probabilidade de que um elétron sofra uma colisão em 
um intervalo de tempo dt é dt/t, onde z é uma constante de tempo. Seja n, O 
número de elétrons que sofreram uma colisão no tempo t, e suponhamos que 
desses elétrons, n não sofreram colisões subsequentes no intervalo de tempo 
t- ty; então: 


dn = ne, (Eg. 9.1) 


n=nge to, (Eq. 9.2) 
De acordo com a Eq. 9.1, o número de elétrons que sofrem colisões entre 
tet+ dt é igual a ndt/t. Para esses elétrons, o tempo entre colisões é t-t;. 
Assim o tempo médio entre colisões” para os ny elétrons é dada por 
dude dt 
<i-lg> = — | (U-to)n— = 1. (Eq. 9.3) 
no q, T 
0 
A velocidade média de deriva pode ser calculada de forma semelhante. À 
velocidade de um elétron que sofreu a última colisão em t, é dada por 
qelt—to) 
v= Do OZ 2% o", (Eq. 9.4) 
onde o primeiro termo representa a velocidade inicial depois da colisão e o 
segundo termo representa a velocidade causada pelo campo elétrico s**. O 
simbolo m* representa a massa efetiva dos elétrons no cristal. A massa efetiva 
é uma grandeza que substitui a massa de um elétron livre (veja o Cap. 6) em 
vários cálculos. Ela reflete o efeito da rede do cristal sobre o comportamento 
do elétron. 
Como a velocidade v, depois da colisão não tem uma direção preferencial, 
a média de w, para muitos elétrons é zero. Para calcularmos o valor médio de » 
no tempo t, entretanto, precisamos conhecer a distribuição de elétrons em função 
de to» 
Seja w a porcentagem de elétrons que ainda não sofreram uma colisão no 
tempo t depois de sofrerem uma colisão em t,. De acordo com a Eq. 9.2, 
w=D =e tok, (Eq. 9.5) 
no 
Assim o número de elétrons que sofreram uma colisão no intervalo de 
tempo entre ty € ty + dt, e ainda não sofreram outra colisão no tempo t é dado 
por 
dio nowdt, = dt 


=(t-to)r Po. 
E: a noe E (Eq. 9.6) 


dn=n 


WEstamos representando o valor médio de uma grandeza / pelo símbolo <f> 
'**Chamamos a atenção do leitor para a simbologia que estamos utilizando: 
E = permissividade, 
E = campo elétrico, 
E = energia, 
e = base dos logaritmos neperianos 
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Nesse caso, a velocidade média de deriva dos elétrons será: 


! 
<p> = | nowit-t)— =— St =—ue. (Eg. 9.7) 
o = 0 


A constante yu = gr/m* é definida como a mobilidade dos elétrons e será 
estudada na Sec. 9.1.3. 

Naturalmente, a velocidade média de deriva é independente do tempo. O 
sinal negativo da Eq. 9.7 indica que a velocidade de deriva dos elétrons tem 
o sentido contrário ao do campo elétrico c. 

De acordo com a Eq. 9.7, a densidade de corrente de deriva dos elétrons 
é dada por 


J,=-ng<v,> = ng pet. (Eq. 9.8) 


No caso de buracos, a equação correspondente à Eg. 9.7 será 
<v,> = H,£, (Eg. 9.9) 


onde o segundo membro é positivo, porque a velocidade de deriva dos buracos 
tem o mesmo sentido que o campo elétrico. Assim, a densidade de corrente 
de deriva dos buracos é dada por 


Jp =Pq<v,> = pqu,e. (Eq. 9.10) 

Quando a corrente de deriva consiste em elétrons e buracos, podemos 
combinar as Egs. 9.8 e 9.10, obtendo 

J=4,+J,=qnu,+puje (Eq. 9.11) 


9.1.2. Resistividade de um semicondutor 


Considere um bloco de um material qualquer de seção reta 4 e compri- 
mento |, ao qual é aplicada uma tensão V. produzindo uma corrente | (Fig. 9.2). 


Figura 9,2 Delinição de resistividade IT 


Por definição, a densidade de corrente J é dada por 

d. 

4 

(A equação J = gnv é uma equação básica que relaciona a densidade de corrente (J), em 


qualquer material, à velocidade dos portadores de carga nesse material (v), à densidade de portadores 
(n) e à carga de cada portador (q) 


Jd = 
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Por outro lado chamando de p a resistividade do material, a resistência entre 
as duas seções transversais será 


Assim, temos 


ou 


porém, 1/4 = Je V/l = e, onde & é o campo elétrico no interior do bloco. Con- 
sequentemente, 


| 
=—"E=õógEé, . 9.12 
: E = 0 (Eq ) 


onde q = 1/p é a condutividade do material. 
Comparando essa expressão com a Eg. 9.11, verificamos que em um se- 
micondutor 


o =aqlun+ up), (Eq. 9.13) 
1 


="———: Eq. 9.14 
qlun + Hp) di 


p 


De acordo com a Eg. 9.14, a resistividade dos semicondutores varia inver- 
samente com a mobilidade e com a concentração de portadores. Assim, a uma 
dada temperatura, a resistividade depende apenas da concentração de impurezas 
(doadoras ou aceitadoras) do material. Na prática, a variação da resistividade 
com a concentração de impurezas só obedece a uma equação simples para 
baixas concentrações de impurezas, caso em que a mobilidade dos elétrons e 
buracos é praticamente constante; para concentrações elevadas, tanto a mo- 
bilidade como a concentração dos portadores dependem da concentração de 
impurezas (veja 9.1.3). 

A Fig. 9.3 mostra a variação da resistividade do silício com a concentração 
de impurezas para cristais tipo P e tipo N. A mobilidade dos portadores se 
mantém constante até a concentração de impurezas atingir valores da ordem 
de 10!º átomos/em”. 


9.1.3. Mobilidade dos portadores em um semicondutor 


A mobilidade dos portadores de carga em um cristal depende das interações 
desses portadores com os átomos da rede cristalina, com os átomos de impu- 
rezas, com defeitos da rede, com outros portadores, etc. Todas essas interações 
são conhecidas pelo nome genérico de mecanismos de espalhamento (scaltering 
mechanisms), 

No caso dos semicondutores usados na prática, que estão relativamente 
livres de defeitos e cuja densidade de portadores não é muito grande, os me- 
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CONCENTRAÇÃO DE IMPUREZAS (CM"?) 
Figura 9.3 Resistividade do silício a 300 K em função da concentração de impurezas. 
[Trvin, J.C., “Resistivity of Bulk Silicon and of Diffused Layers in Silicon”, Bell Syst. 
Tech. J. 41, 387 (1962). Reproduzido com a permissão da American Telephone and Te- 
legraph Company.] 


canismos de espalhamento mais importantes são o espalhamento pela rede (lattice 
scattering) e o espalhamento por impurezas ionizadas (ionized impurity scattering) 
(Fan, 1955). 

O espalhamento pela rede resulta da interação entre os portadores e os 
átomos que constituem a rede cristalina. Como vimos no Cap. 6, os átomos 
da rede cristalina não estão imóveis, mas oscilam em relação à posição de equi- 
líbrio. A amplitude dessas oscilações é diretamente proporcional à temperatura 
do cristal. É fácil ver portanto que o tempo médio entre colisões dos portadores 
com a rede, tu, deve ser inversamente proporcional à temperatura. Por ana- 
logia com a Eq. 9.7, podemos definir uma constante 4, através da equação. 


p= ER, (Eg. 9.15) 


onde uy é a contribuição do espalhamento pela rede para a mobilidade total, 
Ta é O tempo médio entre colisões com a rede e m* é a massa efetiva do portador 
considerado. 

Outro mecanismo importante de espalhamento, principalmente em baixas 
temperaturas, ou para grandes concentrações de impurezas, é o espalhamento 
por impurezas ionizadas. Esse tipo de espalhamento se deve ao fato de que 
tanto os portadores como os átomos ionizados das impurezas possuem cargas 
elétricas. Ao passar pelas proximidades de um átomo ionizado, o portador 
sofre um desvio em sua trajetória, devido à interação entre o campo elétrico 
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do átomo ionizado e a carga do portador. Naturalmente, quanto maior a ve- 
locidade do portador, menor a deflexão sofrida; assim, o espalhamento por 
impurezas ionizadas é mais importante em baixas temperaturas, que implicam 
em uma menor velocidade térmica dos portadores. 

Como no caso do espalhamento pela rede, podemos associar ao espalha- 
mento por impurezas ionizadas uma constante de tempo, t,, o tempo médio 
entre colisões dos portadores com os átomos das impurezas ionizadas. Podemos 
definir também uma constante, 4, através da equação: 

w = El. (Eq. 9.16) 
m 
onde 4, é a contribuição do espalhamento por impurezas ionizadas para a mo- 
bilidade total. Naturalmente, 7, e 4 são inversamente proporcionais à concen- 
tração de impurezas ionizadas. 

Os outros mecanismos de espalhamento só são importantes em circuns- 
tâncias especiais. O espalhamento eletron-eletron, o espalhamento buraco-buraco 
e o espalhamento eletron-buraco, caracterizados por diferentes tempos médios 
entre colisões, 7, .., Ty-n € Te -n+ fesultam da interação direta de portadores e 
só são importantes quando a densidade de portadores é muito elevada. O espa- 
lhamento por impurezas neutras (neutral impurity scattering), caracterizado por 
um tempo entre colisões 7, resulta da interação de portadores com impurezas 
neutras e só é importante em temperaturas muito baixas, quando nem todos 
os átomos das impurezas estão ionizados. 

Se supusermos que os mecanismos de espalhamento que ocorrem no ma- 
terial são independentes entre si (o que nem sempre é verdade), a mobilidade 
total dos portadores pode ser calculada a partir das mobilidades parciais através 
de uma equação bastante simples. Com efeito, basta observar que se T,, Ty, Tas--- 
são os tempos médios entre colisões para os diversos mecanismos de espalha- 
mento, então 1/t,, 1/t,. 1/t,.... representam os números médios de colisões 
por unidade de tempo para os referidos mecanismos. Assim, supondo que os 
mecanismos são independentes, o número médio total de colisões por unidade 
de tempo é dado por 


l/t=1/t + 1/0 + 17, +... (Eg. 9.17) 


Ora, como por definição u = qx/m*, |, = q7,/m*, UH, = q7,/m*,..., temos 


e (Eq. 9.18) 
Hoy Ho Hs 
onde 4 é a mobilidade total e u,, ,, 4;,--- São as contribuições para a mo- 
bilidade total dos diversos mecanismos de espalhamento, 
Quando os mecanismos de espalhamento considerados são apenas o espa- 
lhamento pela rede e o espalhamento por impurezas ionizadas, a Eg. 9.18 se 
reduz a uma expressão muito usada na prática 


datação: (Eq. 9.19) 


A mobilidade que aparece na tabela de constantes no início deste livro é 
a mobilidade do material intrínseco à temperatura ambiente e corresponde 
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praticamente à mobilidade resultante do espalhamento pela rede. Entretanto, 
esse valor também é usado na prática para materiais tipo P e tipo N, desde que 
a concentração de impurezas não seja muito grande. 

Finalmente, completando nosso estudo da mobilidade, vamos discutir a 
influência do campo elétrico sobre a mobilidade dos portadores. 

Voltemos à equação básica que relaciona a velocidade de deriva dos por- 
tadores, v, ao campo elétrico aplicado, &: 

v = ue. (Eq. 9.20) 

Se 4 não dependesse do campo elétrico, o gráfico de v em função de & seria uma 
linha reta de coeficiente angular igual a pu. 


À Fig. 9.4 mostra o que acontece na prática à medida que o campo elétrico 
aumenta. 
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Figura 9.4 Variação da velocidade de deriva com o campo elétrico no silício (Ryder, 
1953) 
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Para campos elétricos fracos, a velocidade aumenta linearmente com o 
campo: para campos elétricos intensos, a velocidade dos portadores passa a 
aumentar mais lentamente, tendendo para um valor constante. Esse fenômeno, 
que é chamado de efeito dos elétrons quentes, se deve ao fato de que os portadores 
não se mantêm em equilibrio térmico com a rede cristalina quando sua energia 
cinética é muito grande (Ryder, 1953). 


9.2. ESTUDO DA DIFUSÃO 


9.21. O processo de difusão 


Ao calcularmos a densidade de corrente de deriva, Eq. 9.11, não levamos 
em conta a variação espacial das concentrações de portadores, ne p. Sen ep 
variarem ao longo do cristal, haverá uma outra contribuição para a corrente 
que será somada à corrente de deriva. 

Consideremos o movimento térmico de partículas no espaço unidimensional 
da Fig. 9.5, onde v, é a velocidade térmica das partículas e | é o caminho livre 
médio, definido por | = v,t. (1 = tempo médio entre colisões). 

O número de partículas que atravessam a linha A da esquerda para a direita 
por unidade de tempo e por unidade de área é igual a n,0,/2, enquanto o nú- 
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W “ 


vt vt Figura 9.5 Processo de difusão 


ma o n, 


mero correspondente da direita para a esquerda é n,v,/2, onde n, e n, são as 
concentrações de partículas à esquerda e à direita de linha A, respectivamente. 
O fluxo total de partículas por unidade de tempo e por unidade de área é dado por 


=. (Eq. 9.21) 


S=u, 


Como apenas as partículas que estão a uma distância igual ou menor que o 
caminho livre médio da linha 4 conseguem atravessá-la, a diferença n, —n, deve 
ser calculada a essa distância. Assim, n, = ny + (ôn/ôx)l e n;=no (ôn/ôx)l, 
tomando apenas os primeiros dois termos da expansão em série de Taylor. 
Substituindo na Eq. 9.21, obtemos 


S = me (Eg. 9.22) 


Se as partículas possuem carga elétrica, o movimento dessas partículas 
constitui uma corrente, cuja densidade é dada por 


J,=-g5= quis! = e E = qD = , (Eq. 9.23) 
onde, por definição, | = v,t, m*v; = kT. No caso tridimensional, a Eg. 9.23 
fornece a componente x da corrente de difusão se substituirmos v, pela com- 
ponente x da velocidade, v,.. A corrente J definida pela Eg. 9.23 é chamada 
de corrente de difusão. A constante D é chamada de coeficiente de difusão. Re- 
pare que de acordo com as Egs. 9.7 e 9.23 existe uma relação entre D e q: 
D = kTu/q. Essa relação é chamada de relação de Einstein. 


9.3. SUPERPOSIÇÃO DOS PROCESSOS DE DERIVA E DIFUSÃO 


De acordo com as Egs. 9.11 e 9.23, a densidade total da corrente de elétrons 
ao longo do eixo dos x em um semicondutor é dada por 


ên 


Ja =4D, x + qu,NE,- (Eq. 9.24) 
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onde o primeiro termo representa a corrente de difusão e o segundo a corrente 
de deriva. Do mesmo modo, a corrente de buracos é dada por 


ô 
Jo =D, is +quppê,- (Eq. 9.25) 


EXERCÍCIOS 


9.1. Sugira um método para medir a condutividade de um semicondutor. 

9.2. Descreva dois métodos para verificar experimentalmente se um semicondutor é 
tipo P ou tipo N. 

9.3. Calcule a resístividade de um cristal de silício com 1,4 x 10!º átomos de fósforo 
por emº e 1,0 x 10*º átomos de boro por cm”, à temperatura de 27ºC. 

frig,4. Uma pequena concentração de portadores em minoria é injetada em um ponto 
de um cristal semicondutor homogêneo. Um campo elétrico de 10 V/cm é apli- 
cado ao cristal; esse campo faz com que os portadores injetados percorram uma 
distância de 1 cm em 250 us. Calcule a velocidade de deriva e o coeficiente de di- 
fusão dos portadores em minoria. 

9.5. Em um cristal semicondutor, os elétrons têm uma mobilidade de 500 cm?/V -s. 
Calcule o tempo médio entre colisões. Para um campo elétrico de 100 V/cm, calcule 
também a distância média percorrida pelos elétrons entre duas colisões sucessivas. 
Nesses cálculos, use m* = m, onde m é a massa do elétron livre. 

89.6. A resistência de um cristal de germânio tipo P (figura a seguir), entre os planos 
A e B. é igual a 1000 a 300ºK. Calcule a concentração de elétrons e buracos na 
amostra. 

Dados: 
para o germânio a 27ºC, n,= 24 x 10 cm”?: 
H, = 1900em*/V-s: 4, = 3900 cm?/V-s. 


> 


1cm 


9.7. Considere um bloco retangular de GaAs de dimensões | mm x 1 mm (seção) e 
comprimento / = 1 cm, ao qual é aplicada uma tensão de 100 V no sentido do 
comprimento do bloco. Sabendo que o material possui 10!* elétrons/cm? e que 
a mobilidade dos elétrons é de 8 x 10º cm*/V-s, calcule a corrente que circula 
no material. 

(*9.8. Considere um semicondutor intrínseco. Raciocinando qualitativamente, podemos 
dizer que a condutividade do material deve aumentar com a temperatura, já que à 
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medida que a temperatura aumenta o número de ligações covalentes que são 
rompidas também aumenta. Prove matematicamente que isso é verdade, isto é, 
que da/dT é maior que zero para qualquer temperatura. 

9.9. Mostre que a velocidade média dos elétrons em um semicondutor submetido a 
um campo elétrico alternado da forma « = Re [spe'” | é dada por 


<p> = Re Ra - gl 
mt I+imor* 


REFERÊNCIAS 


FAN,H. Y., “Valence Semiconductors, Germanium and Silicon”, em F. Seitz e D. Turnbull 
(eds.), Solid State Physics, Vol. 1 (New York: Academic Press, 1955). 

RYDER, E.J., “Mobility of Holes and Electrons in High Electric Fields”, Phys. Rev. 90, 
766 (1953). 


BIBLIOGRAFIA 


Conwell, E. M., High Field Effects in Semiconductors (New York: Academic Press, 1967). 
Paige. E.G.S., The Electrical Conductivity of Germanium, em Progr. Semiconductors, 
Vol. 8 (1964). 


10 Junções PN 


Vamos agora aplicar os conhecimentos desenvolvidos nos capítulos ante- 
riores ao estudo de um dispositivo semicondutor muito usado na prática: a 
junção PN, cuja estrutura está representada de forma esquemática na Fig. 10.1. 

Como está indicado na figura, uma junção PN é um bloco de material 
semicondutor que consiste em duas regiões distintas, uma tipo P e outra do tipo N. 


Figura 10.1 Representação esque- 
mática de uma junção P-N 


10.1. JUNÇÃO PN EM EQUILÍBRIO 


Suponhamos, apenas para efeito de discussão, que acabamos de formar 
uma junção PN, simplesmente “justapondo” um bloco de material tipo P em 
um bloco de material tipo N. É claro que existem mais elétrons do lado N e 
mais buracos do lado P; assim, existe uma região em que ên/ôx e ôp/ôx são 
diferentes de zero. Surge naturalmente a pergunta: haverá corrente, mesmo na 
ausência de campo elétrico externo? A resposta é sim, mas por tempo limitado. 
Se os elétrons e buracos não tivessem carga elétrica, as correntes de difusão 
continuariam até que as distribuições de elétrons e buracos fossem uniformes 
em todo o cristal. Mas como os elétrons e buracos possuem cargas elétricas e 
os cristais contêm os átomos ionizados das impurezas, a situação é bem dife- 
rente. À medida que, em consegiência das correntes de difusão, os elétrons do 
lado N passam para o lado P e os buracos do lado P passam para o lado N, 
este vai se tornando mais positivo devido à presença de doadores ionizados, e 
o lado P vai se tornando mais negativo devido à presença de aceitadores ionizados. 
Assim, é criado um campo elétrico nas proximidades da junção, como vemos 
na Fig. 10.2. 

Esse campo elétrico tem um sentido tal que tende a manter os buracos no 
lado P da junção e os elétrons no lado N. Assim, depois que um certo número 
de elétrons e buracos atravessa a junção, é criado um campo elétrico que inter- 
rompe o movimento de portadores. 


Figura 10.22 Distribuição de cargas 
no interior de uma junção P-N em 


equilibrio 
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É assim que se atinge o equilibrio, no qual a corrente total é zero; as cor- 
rentes de difusão dos elétrons e buracos são equilibradas exatamente pelas 
correntes de deriva devido ao campo elétrico interno. 

Vamos agora traduzir essa afirmação em termos matemáticos. De acordo 
com a Eq. 9.25, a corrente de buracos é dada por 


ô 
J, = -aD, 5 + qu, pe. [Eg. 9.25] 


Substituindo nessa expressão a concentração de buracos por seu valor 
em termos do nível de Fermi intrínseco (Eg. 8.29) 


Rd td [Eq. 8.29] 
e sua derivada 
dp  p [dE dE, | 
ak | (Eq. 10,1) 
temos 
D dE, dE 
= -—g—E pai Rs , 5 
Jo 4kT DE | + qu, Pê. (Eq. 10.2) 
Sabemos da eletrostática que, em uma dimensão, 
1 0E : 
“q ôx (Eg. 10.3) 


onde : é o campo elétrico, E é a energia potencial da particula considerada c q 
é a carga dessa partícula. Como vimos nos capítulos anteriores, a energia po- 
tencial do elétron é representada pela energia na borda da faixa de condução, 
E.. e a energia de um buraco é representada pela energia na borda da faixa de 
valência, E, . Entretanto como estamos interessados apenas na derivada da ener- 
gia potencial, podemos usar qualquer nível de energia que seja paralelo às faixas 
de condução e de valência. Em particular, podemos escolher o nível de Fermi 
intrinseco e escrever 


p=— CE. (Eq. 10.4) 


Substituindo a Eq. 10.4 em 10.2 e usando a relação de Einstein q, = qD KT, 
obtemos 
D. dE dE 
o Yi me 1 ado ts SO R 
J,=q REP de = tor (Eg. 10.5) 


Assim, para que a corrente total de buracos seja zero é preciso que o nível 
de Fermi não varie ao passarmos de um lado para outro da junção. 

Se aplicarmos o mesmo tratamento matemático à corrente de elétrons, 
obteremos a expressão 


=p E = gies=! 4/4 
no — le (Eq. 10.6) 
e portanto para que a corrente total de elétrons seja zero o nivel de Fermi tam- 


bém deve ser constante. Assim, em uma junção PN em equilibrio o nivel de Fermi 
e constante ao longo de todo o cristal. 
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Para calcularmos as caracteristicas da região de carga espacial que é for- 
mada nas proximidades da junção, vamos supor que ela consiste inteiramente 
nas cargas dos aceitadores e doadores ionizados. Isso equivale a dizer que a 
região de carga espacial está totalmente desprovida de portadores, ou melhor, 
que as concentrações de portadores n e p são desprezíveis em comparação com 
a concentração de impurezas na região de carga espacial. 

Na Fig. 10.3a) vemos a distribuição de cargas em uma junção PN em 
função da distância, para o caso em que a concentração de impurezas nas regiões 
p en é uniforme. 


Eres 


Ec 


ATE E | 
: O Ma Ei 
Figura 10.3 Distribuição de cargas, campo E ad 


elétrico e níveis de energia em uma junção 
abrupta (c) 


Uma junção desse tipo é chamada de junção abrupta. As densidades de carga 
nos lados N e P da junção são dadas por qN, e —yN,. respectivamente, Às 
larguras das duas partes da região de carga espacial são x, e x, e a largura total 
(x, + x,) é chamada de W A distribuição de campo elétrico, obtida por inte- 
gração da distribuição de carga, está representada na Fig. 10.3(b). Integrando 
a distribuição de campo elétrico, obtemos a energia E, dos elétrons em função 
da distância [Fig. 10.%c)]. Naturalmente, de acordo com a Eg. 10.4, as bordas das 
faixas de condução e de valência acompanham as variações do nível de Fermi 
intrínseco E,. 

Vamos agora calcular a curvatura das faixas de energia. isto é, a variação 
da energia quando passamos do lado P para o lado N. Examinando a Fig. 10.30), 
vemos que, como o nível de Fermi é constante, a variação de energia é a soma 
de duas parcelas, E,, e E,,, que são definidas do seguinte modo: 


Epp = E,— Eo)etada é » Em = E Eliomar (Eq. 10.7) 


Para calcularmos os valores de E,, e E,,, vamos supor que na região P 
neutra, isto é, fora da região de carga espacial, p = N,, enquanto que na re- 
gião N neutra, n = N,. Nesse caso, as Egs. 8.28 e 8.29 fornecem 


Es, =—kTln Na. (Eq. 10.8) 


Es, = +Tin Do. (Eq. 10.9) 

i 
A variação total de energia E, em equilíbrio será a soma de E,, € E,,: 
E; = Epp + Em (Eg. 10.10) 
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Como sabemos da eletrostática, o potencial elétrico de um elétron é dado 
por 
E 
Fo) = “a (Eq. 10.11) 
onde E é a energia do elétron e q é a sua carga. Assim, a variação total do po- 
tencial elétrico & através de uma junção PN em equilibrio é dada por 


E KT, NaNp 
=- = - 1012 
Es e La a de pass 


O potencial p, é chamado de potencial de difusão da junção PN. 
Como o campo elétrico nas regiões neutras do semicondutor deve ser zero, 
as cargas totais dos dois lados da junção devem ser iguais e opostas, isto é, 
NoXa = NaXp (Eg. 10.13) 
O campo elétrico máximo da junção é dado por 
o 4N pXa = GN aX, 
E SR O 
onde « é a permissividade do material”). 


A variação total de potencial elétrico é dada pela área sob a curva da Fig. 
10.3(b), ou seja 


(Eq. 10.14) 


PRE, (Eg. 10.15) 


2 Cmax 


Combinando as Egs. 10.13, 10.14 e 10.15 e lembrando que x, + x, = W 
vemos que nesse caso a largura da região de carga espacial é dada por 


W = (Eae d Na É ra, (Eq. 10.16) 


onde, em equilíbrio, pp = Pp. 
Voltando à Eq. 9.25, vemos que em equilibrio térmico, ou seja, quando 
Jp, = O, temos 


qu, pe — ab, é? (Eq. 10.17) 


ou, usando a relação de Einstein Ho = D,/kT, 
E: X go 4 Ea. 10.18 


Integrando a Eqg. 10.18 de 4 a B [Fig 10.3%a)], já que no resto do cristal 
o campo elétrico é zero, temos 


B r â B 
id O edx; (Eq. 10.19) 
Me. kT P 
PB) abr. 
m AV et? (Eg. 10.20) 
P po = etómikT, (Eq. 10.21) 
Pno 


'*"'Novamente lembramos que e = permissividade, ; = campo elétrico 
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onde p,, € Pr, São as concentrações de buracos em equilibrio térmico nas regiões 
Pe N, respectivamente. Para os elétrons, a relação equivalente a Eq. 10.21 será 
n 


—Po  g-aerikT (Eq. 10.22) 


Ro 


onde n,, € Ny, São as concentrações de elétrons nas regiões P e N 


10.2. JUNÇÃO POLARIZADA 


10.2.1. Polarização direta e polarização inversa 


Quando aplicamos à junção PN uma tensão V tal que o lado N fique po- 
sitivo em relação ao lado P, a variação total de potencial elétrico através da 
junção aumenta de W. Assim, 


dr =p + 4 (Eg. 10.23) 
Em conseqiiência, a largura da região de carga espacial aumenta. Essa 


situação está ilustrada na Fig. 10.4. Nessas condições, dizemos que a junção 
está polarizada inversamente. 


(a) 


Be+VI VI 
Figura 10.4 Distribuição de carga, campo 
elétrico e níveis de energia em uma junção fem — 
PN polarizada inversamente 


= (e) 

Substituindo a Eq. 10.23 na Eg. 10.16, podemos calcular a largura da região 
de carga espacial em função da tensão V, aplicada. 

Quando a junção é polarizada diretamente, isto é, quando aplicamos uma 
tensão V, tal que o lado N fique negativo em relação ao lado P, a variação total 
de potencial elétrico diminui e portanto a largura da região de carga espacial 
tende a diminuir. Entretanto como veremos mais tarde, em uma junção po- 
larizada diretamente a corrente é muito grande — imensamente maior do que 
no caso da polarização inversa. Assim a concentração de portadores na região 
de carga espacial se torna considerável e portanto a aproximação de que a carga 
espacial se deve apenas aos doadores e aceitadores ionizados, usada para de- 
duzir a Eg. 10.16, perde sua validade. 


10.2.2. Junção polarizada inversamente 


Em condições de equilibrio, pares elétron-buraco são gerados constan- 
temente no cristal. Na ausência de campo elétrico externo, os pares elétron- 
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-buraco se recombinam e a corrente total é zero. Entretanto quando aplicamos 
uma tensão positiva à região N da junção PN, como na Fig. 10.4, os pares, uma 
vez formados, são imediatamente separados e sua probabilidade de recombi- 
nação diminui. Esse é o mecanismo básico que dá origem à corrente nas junções 
PN polarizadas inversamente. Todas as correntes em uma junção polarizada 
inversamente são devidas a pares elétron-buraco gerados no semicondutor. 
Esses pares podem ser gerados em duas regiões distintas do cristal: a região 
de carga espacial e as regiões neutras P e N longe da junção (Fig. 10.5). 


Figura 10.5 Junção PN polarizada inversamente, nO) = pO) = O 


A corrente que resulta da geração de pares clétron-buraco na região de 
carga espacial é chamada de corrente de geração. A corrente que resulta de pares 
criados nas regiões neutras é chamada de corrente de difusão. A principal dife- 
rença entre essas correntes é que enquanto os pares gerados na região de carga 
espacial estão submetidos a um campo elétrico e se movem principalmente por 
deriva, os pares gerados nas regiões neutras não sofrem a ação de campo elé- 
trico algum e portanto podem mover-se apenas por difusão. 

Demonstra-se que a densidade de corrente de geração é dada aproxima- 
damente por 


Dm 


E (Eq. 10.24) 


onde n; é a concentração intrinseca de portadores do material, 7 é o tempo 
médio de vida dos portadores (tempo que os portadores levam para se recom- 
binarem. depois que são gerados) e W é a largura da região de carga espacial. 

A corrente de geração é proporcional à intensidade da tensão aplicada: 
quanto maior a tensão, maior a largura da região de carga espacial, W, e por- 
tanto maior o número de pares produzidos nesta região. Esta corrente também 
aumenta com a temperatura, pois n; aumenta com a temperatura. 

A densidade de corrente de difusão é dada aproximadamente por 


D, [D 
Jau = NE [non AO] + q e Lp po — PAO]. (Eq. 10.25) 


onde n,, € Pa, São as densidades de portadores em minoria longe da região de 
carga espacial e n,(0) e p,(0) são as densidades de portadores em minoria na 
fronteira entre a região neutra e a região de carga espacial (Fig. 10.5). 
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Quando pares elétron-buraco são gerados na região N neutra, os buracos 
que chegam à região de carga espacial são levados imediatamente para a região 
P pelo campo elétrico que existe na região de carga espacial. Do mesmo modo, 
os elétrons da região P neutra que chegam à região de carga espacial são logo 
carregados para a região N. Assim, para uma junção polarizada inversamente, 


nÃO) = p (0) = 0. (Eg. 10.26) 
E a Eq. 10.25 se reduz a 


D, D, , 
Jar = Mp E + GP = (Eq. 10,27) 


onde n,, € Py São as densidades de portadores em minoria. 

Para avaliarmos a importância relativa dessas duas correntes, vamos 
calcular a relação entre elas. De acordo com as Fgs. 10.24 e 10.27, temos 
ai Npo + Pno Det ' 


=2 


ger i 


2 
Como nn, = ProPpo =; € tomando n,, = Na» Pro = Np» temos 


J NA NN/D 
e = (EN * 5 Dr. (Eq. 10.28) 


J Na Np 


Essa relação é portanto diretamente proporcional à concentração intrinseca 
de portadores, n,. Assim à medida que a temperatura aumenta, a corrente de 
difusão se torna mais importante. 

A corrente total em uma junção PN polarizada inversamente, que é cha- 
mada de corrente inversa, é dada por 


= AJ, + Tou): (Eq. 10.29) 


ger 


onde 4 é a área da junção. 


10.2.3. Junção polarizada diretamente 


Quando polarizamos diretamente uma junção PN, como na Fig. 10.6, a 
variação total de potencial elétrico, &,., diminui. Em equilíbrio, a condição de 
corrente zero era estabelecida por uma compensação exata entre as correntes 
de deriva e difusão de cada portador. Como em uma junção polarizada di- 
retamente a variação de potencial elétrico e portanto o campo elétrico na região 
de carga espacial é menor do que em equilíbrio, a corrente de deriva dos por- 
tadores diminui e o equilibrio entre as correntes de deriva e as correntes de 
difusão é rompido, o que torna a corrente total diferente de zero. 

É importante observar que embora a diferença de potencial entre os dois 
lados da junção seja menor do que em equilíbrio, não desaparece totalmente. 
Em consegiiência, o campo elétrico na região de carga espacial ainda se opõe 
a que os elétrons no lado N e os buracos do lado P atravessem a junção. Assim, a 
corrente será contra o sentido do campo elétrico. Isso é possivel graças aos 
enormes gradientes de concentração ôn/0x c dp/ôx que existem perto da junção. 

Outro modo de encarar o fenômeno é o seguinte: quando aplicamos um 
potencial negativo ao lado N e um potencial positivo ao lado P da junção, os 
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E ae (a |, na 
Figura 10.6 Níveis de energia em 
to: E a Eve uma junção PN polarizada direta- 


EVN mente 
buracos da região P são “empurrados” em direção ao lado N e os elétrons da 
região N são “empurrados” em direção ao lado P. Os buracos e elétrons se 
encontram nas proximidades da junção, onde se recombinam em grandes quan- 
tidades. 


À primeira vista, pareceria que em pouco tempo não haveria mais por- 
tadores no interior do cristal. Entretanto a fonte de tensão externa está cons- 
tantemente injetando elétrons na região N e buracos na região P. Assim, os 
elétrons e buracos que se recombinam são imediatamente substituídos por 
outros. Esse movimento de portadores constitui uma corrente. 

De acordo com essa segunda interpretação, a corrente direta resulta da 
recombinação de pares elétron-buraco, do mesmo modo que a corrente inversa, 
examinada no parágrafo anterior, resultava da geração de pares elétron-buraco 
no cristal, 

Como no caso da corrente inversa, podemos dividir a corrente em duas 
componentes. À corrente que resulta da recombinação de pares elétron-buraco 
na região de carga espacial é chamada de corrente de recombinação; a corrente 
que resulta de recombinação nas regiões neutras P e N é chamada de correnie 
de difusão. 

Demonstra-se que a densidade de corrente de recombinação é dada apro- 
ximadamente por 

ue qr We Dizer, (Eq. 10.30) 
onde V, é a tensão aplicada. 

A corrente de recombinação varia exponencialmente com a tensão aplicada. 

A densidade de corrente de difusão é dada pela mesma expressão que a 
usada para o caso da polarização inversa, Eg. 10.25. Entretanto, no caso da 
polarização direta, as concentrações dos portadores em minoria na fronteira 
entre a região neutra e a região de carga espacial, n,(0) e p,(0), são diferentes 
de zero (Fig. 10.7). 


Figura 10.7 Junção PN polarizada diretamente, p,(0), n,(0) + O 
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Vamos calcular esses dois valores. Em primeiro lugar, é evidente que a 
diferença entre o nível de condução do lado P e o nível de condução do lado 
N. E,,— E,» que era igual a &, em equilibrio, passa a ser igual a (), — Vo) quando 
polarizamos diretamente a junção (Fig. 10.6). Se supusermos que uma equação 
semelhante à Eq. 10.22 ainda é válida, a concentração de elétrons em minoria 
do lado P será dada por 


nO) = ne” “= Vo)lkT = Dis AMT gra FE = RT, (Eq. 10.31) 
Do mesmo modo, a concentração de buracos do lado N será dada por 
PO) = Duo 8" RT, (Eq. 10.32) 


Substituindo as Egs. 10.31 e 10.32 em 10.25, temos 


Jari = and + Aros) + de | paraar 9 à (Eq. 10.33) 


Como no caso da polarização inversa, vamos calcular a relação entre as 
duas correntes. De acordo com as Egs. 10.30 e 10.33 temos 


Jair No T Pro E Pro. o A A D,t VokT = 
=2 + le li de 


Fe h; 
D, VokT 1 
sd ( Mo ya Pas [AF A, (Eq. 10.34) 
Na No) Wº ez 


No caso da polarização direta, o fator mais importante da Eg. 10.34 é o 
fator exponencial. Em geral, para pequenos valores de tensão aplicada, a cor- 
rente de recombinação é a mais importante. À medida que V, aumenta, a influ- 
ência da corrente de difusão se torna cada vez maior. 

A corrente total em uma junção PN polarizada diretamente, que é cha- 
mada de corrente direta, é dada por 


Ip = AJ + Tam): (Eq. 10.35) 
onde 4 é a área da junção. 


10.3. EQUAÇÃO DO DIODO IDEAL 


Vamos agora calcular a relação entre a corrente direta e a corrente inversa. 
De acordo com as Eqs. 10.24, 10.27, 10.29, 10.30, 10.34 e 10.35, temos 


dpi ii We gVp/ZkT + Nou [= E C Pao [ne 2) err DKkT 1) 
Y ln, w 5 f To e 
2 T + No tas + Dao - 


Se supusermos que tanto a corrente direta como a corrente inversa são 
dominadas pelas correntes de difusão, teremos 


Ip 


I 
Ip = ("287 1), (Eq. 10.36) 


== o” 


sas esF DRT | ' 
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D, D, 
L=-h="|n5/ + JA. 


A relação entre corrente e tensão dada pela Eq. 10.36 é chamada de equação 
do diodo ideal e está representada graficamente na Fig. 10.8, Para tensões inversas 
W > kT/q. a corrente tende para o valor constante 1, = -1,. A corrente [, é 
chamada de corrente inversa de saturação. Para tensões diretas V, > KkT'/q, a 
corrente aumenta exponencialmente: 1, = 1, e”, 


onde 


Figura 10.8 Curva característica do 
diodo ideal (Eq. 10.36) 


10,4. CAPACITÂNCIA DE UMA JUNÇÃO PN 


A capacitância dinâmica por unidade de área de um material qualquer de 
seção reta uniforme é definida por 


Cate, (Eq. 10.37) 


onde dQ é a variação incremental de carga por unidade de área e dV é a va- 
riação incremental de tensão aplicada. 

No caso de uma junção PN, uma variação na tensão aplicada faz variar 
a carga espacial e portanto a junção PN possui uma capacitância diferente 
de zero. 

Para calcularmos o valor desta capacitância, suponhamos que a carga 
espacial total de um dos lados da junção é Q (do outro lado da junção, natu- 
ralmente, a carga será —Q). Nesse caso, quando aumentamos a tensão aplicada 
à junção de uma quantidade dV. a carga aumenta de dQ. De acordo com a 
Equação de Poisson, a variação de campo elétrico é dada por 

d 


de=, (Eq. 10.38) 


onde e é a permissividade do material. 
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Por outro lado, a variação dV da tensão aplicada é dada por 
dV = Wide, (Eq. 10.39) 


onde W é a largura da região de carga espacial. 

Combinando as Eqs. 10.37, 10.38 e 10.39 obtemos 

E 
ai (Eg. 10.40) 

A Eg. 10.40 é válida para qualquer distribuição de impurezas. 

Como exemplo, vamos calcular a capacitância de uma junção abrupta. 

De acordo com a Eg. 10.16, a largura da região de carga espacial de uma 
junção abrupta é dada por 


— PEN+N 
W=.| E ER br. [Eg. 10.16] 


Combinando a Eg. 10.16 com a Eq. 10.40, obtemos 


NAN qe NAN 
PA EL 22/44 PS RE, o Ea. 10.41 
VAN, END VMB FD) Na rn (Ea 1041) 


onde V é a tensão externa aplicada à junção. Repare que W é positivo quando 
a polarização é inversa e negativo quando a polarização é direta. 

A Eq. 10.41 mostra um fato muito importante: a capacitância de uma 
junção PN depende da tensão externa aplicada. Embora esse fato tenha sido 
demonstrado apenas para a junção abrupta, ocorre também em outros tipos de 
junção. Esse fenômeno é utilizado em uma classe de dispositivos semicondutores, 
os chamados varactores, que são usados como capacitores variáveis. Os varac- 
tores não passam de junções PN cuja capacitância é controlada por uma tensão 
externa. 


10.5. RUPTURA DE UMA JUNÇÃO PN 


Como vimos, a junção PN apresenta uma baixa resistência quando é po- 
larizada diretamente e uma resistência elevada quando é polarizada inversa- 
mente. Entretanto todas as junções PN também apresentam uma baixa resis- 
tência quando são polarizadas inversamente com uma tensão suficientemente 
elevada. Esse fenômeno é chamado de ruptura da junção. A ruptura pode ser 
causada por dois mecanismos, que estudaremos a seguir. 


10.5.1. Efeito Zener 


Quando um semicondutor é submetido a um campo elétrico clevado, como 
o que existe por exemplo na região de carga espacial de uma junção PN pola- 
rizada inversamente, as ligações covalentes ficam distorcidas. Se o campo clé- 
trico for suficientemente intenso, as ligações podem ser rompidas, com a formação 
de pares elétron-buraco, Essa ruptura de ligações covalentes por efeito direto 
do campo elétrico é chamada de efeito Zener. 
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10.5.2. Efeito de avalanche 


Outro mecanismo que pode provocar a ruptura da junção é o efeito de 
avalanche, o qual ocorre quando os portadores gerados termicamente no interior 
da região de carga espacial são acelerados pelo campo elétrico aplicado e adqui- 
rem suficiente energia cinética para romperem ligações covalentes. A ruptura 
de ligações covalentes produz novos portadores, que por sua vez são acelerados 
pelo campo, rompem novas ligações covalentes, etc. Devido a essa “reação em 
cadeia”, o fenômeno é chamado de efeito de avalanche. 

É evidente que a energia cinética adquirida pelos portadores entre duas 
colisões é diretamente proporcional à tensão inversa aplicada e ao caminho 
livre médio dos portadores. Para baixas concentrações de portadores, o ca- 
minho livre médio é maior e a ruptura é provocada pelo efeito de avalanche; 
para altas concentrações, a ruptura, em geral, é provocada pelo efeito Zener. 

A Fig. 10.9 mostra a curva característica de uma junção PN real, onde 
podemos ver claramente o fenômeno da ruptura. 


Vo 


Figura 10.9 Curva caracteristica de 
uma junção PN. (a) região de condução 
direta. (b) região de condução inversa; 
(c) região de ruptura 


Ir 


EXERCÍCIOS 


810.1. Em uma junção PN abrupta de germânio, a concentração N , de átomos aceita- 
dores do lado P é igual à concentração N, de átomos doadores do lado N: 
Na = Np = 10"* átomos/cm'. 
Calcule a largura da região de carga espacial: 
a) para uma polarização direta de 0,1 V; 
b) para uma polarização inversa de 0,1 V. 
Suponha que todas as impurezas estão ionizadas e que a temperatura da junção 
e 300K. 
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“EO. 


10.4. 


Dados: para o germânio a 300 K, n, = 2,4 x 10 em”?; 
e=163€e,= 14x 107!2 F/em. 


- Na figura a seguir, vemos a distribuição de cargas em uma junção PN em função 


da distância, para o caso em que a distribuição de impurezas varia linearmente 
entre as duas regiões. Uma junção desse tipo é chamada de junção gradual. A 
distribuição de cargas (e de impurezas) em uma junção gradual é caracterizada 
pelo gradiente de concentração «, definido pela equação 


plx)=-gax (-WI<x< W2) 


onde p(x) é a densidade de carga e q é à carga do elétron. 


Px) 


Para uma junção gradual. pede-se 

a) o potencial de difusão à, da junção; 

b) a distribuição de câmpo elétrico s(x) na região de carga espacial, 

c) a largura da região de carga espacial quando uma tensão inversa V, é aplicada 
externamente à junção, 

Dados 

1. O gradiente de concentração, a. 

2. A largura da região de carga espacial em equilíbrio, W,. 

3. A temperatura da junção, T. 

Suponha conhecidos os parâmetros do material de que é feita a junção. 

Duas junções PN são ligadas em oposição, como na figura a seguir. Pede-se 
a) a tensão entre os terminais de cada uma das junções; 

bj a corrente no circuito. 

Dados 

1. As duas junções obedecem à equação do diodo ideal. 

2. A corrente inversa de saturação das duas junções é IO MÁ. 

3. À temperatura das junções é 27ºC. 


+ Ah & do Vo = 


10V 


Um meio de detetar a presença de radiações — raios gama, por exemplo — é 
utilizar uma junção inversamente polarizada. Se a tensão inversa é conveniente- 
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mente escolhida, a largura da região de transição é tal que os fótons dos raios 
gama perdem toda a sua energia dentro da região de transição. 

a) Explique como funciona este detetor de radiação. 

bj Qual a necessidade de a partícula perder toda a sua energia dentro da região 
de transição” 

c) Que acontece quando a corrente de fuga é muito elevada? 

d) Que é mais vantajoso, dopar mais fortemente a região N (lado em que a'radiação 
incide) ou a região P? Justifique. 
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11 Transistores de junção 


O dispositivo semicondutor mais usado hoje em dia é sem dúvida o tran- 
sistor de junção. Sua invenção revolucionou a eletrônica, tornando possível a 
fabricação de aparelhos mais compactos, mais eficientes e mais econômicos. 
Neste capítulo vamos estudar os princípios de funcionamento dos transistores 
de junção, baseando-nos na teoria das junções PN, estudada no capítulo anterior. 


11.1. Estrutura do transistor de junção 


A estrutura do transistor de junção está representada dé forma esquemática 
na Fig. 11.1. Como podemos ver na ilustração, o transistor de junção é um bloco 
de material semicondutor que consiste de três regiões distintas. Existem dois 
tipos de transistores de junção. No primeiro, o tipo PNP, Fig. 11.1(a), a região 
central é do tipo N e as regiões externas são do tipo P. No segundo, o tipo NPN, 
a região central é do tipo P e as regiões externas são do tipo N. 


Figura 11.1 Representação es- 
quemática da estrutura de tran- EEE ERES 
sistores de junção: (a) PNP; (b) 

NPN (a) Cb) 


De acordo com o que vimos no capitulo anterior, o transistor pode ser 
considerado como a combinação de duas junções PN (Fig. 11.2). 


BENGEGRAR 
efa pe o = [5 pes 


Figura 11.2 Combinação puramente esquemática de duas junções, compondo um 
transistor 


11.2. TRANSISTOR EM EQUILÍBRIO 


Quando o transistor não está submetido a tensão externa alguma, ele se 
comporta como duas junções PN em equilíbrio. Nesse caso, como vimos no 
capitulo anterior, aparecem campos elétricos nas proximidades das junções que 
equilibram as correntes de difusão de elétrons e buracos. A Fig. 11.3 mostra 
o diagrama de níveis de energia de um transistor PNP em equilibrio. 


P N Pp 
Figura 11.3 Diagrama de niveis de Ec 
energia em um transistor PNP em A; SR es Er 
equilibrio Rae 1 te Ev 
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11.3. TRANSISTOR POLARIZADO 


Suponhamos agora que polarizamos diretamente uma das junções e pola- 
rizamos inversamente a outra junção (Fig. 11.4). 


Figura 11.4 Polarização das junções 
de um transistor P N 


Figura 11.5 Diagrama de níveis de energia em um 
transistor polarizado 


Nesse caso, o diagrama de níveis de encrgia será o representado na Fig. 11.5. 

A parte central do transistor é chamada de base; a parte externa da junção 
polarizada diretamente é chamada de emissor; a parte externa da junção pola- 
rizada inversamente é chamada de coletor. 

À primeira vista poderia parecer que em um transistor polarizado, como 
na Fig. 11.4, a corrente de emissor seria muito maior que a corrente de coletor, 
Já que a junção emissor-base está polarizada diretamente e a junção base-coletor 
está polarizada inversamente. Entretanto, se a base for suficientemente estreita, 
a corrente de coletor será quase igual à corrente de emissor. 

Para compreendermos por que isso acontece. vamos examinar a Fig. 11.6. 

Emissor Base Coletor 
TipoP TipoN TipoP 


Figura 11.6 Funcionamento de um 
transistor PNP 


Repare que quando a junção emissor-base está polarizada diretamente, um 
grande número de buracos atravessa a junção c penetra na região da base. Alguns 
desses buracos se combinam com os elétrons da base tipo N, como em uma 
junção PN comum; entretanto, se a base for muito estreita, a maioria dos buracos 
chegará a junção base-coletor. Assim que cruzam a junção base-coletor, esses 
buracos são atraídos para o terminal do coletor, 

Como vimos no capitulo anterior, os buracos que deixam o emissor são 
imediatamente substituídos por novos buracos criados nas vizinhanças do ter- 
minal de emissor, que retira elétrons da região tipo P. Da mesma forma, os 
elétrons perdidos pela base ao se recombinarem com os buracos que cruzam a 
junção emissor-base são substituídos por novos elétrons, introduzidos através 
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do terminal de base. Os buracos que chegam ao terminal de coletor se combinam 
com elétrons injetados por esse terminal na região tipo P. 

Ora o movimento dos elétrons constitui uma corrente elétrica. Usando o 
sentido convencional da corrente (sentido contrário ao do deslocamento dos 
elétrons), podemos apresentar as correntes em um transistor PNP como na 
Fig. 11,7. 


Emissor | Base | Coletor 
Tipo N | Tipo P| Tipo N 


Emissor | Base | Coletor 
Tipo P | TipoN | Tipo P 


Figura 11.7 Correntes em um tran- Figura 11.8 Correntes em um tran- 
sistor PNP sistor NPN 


No caso de um transistor NPN, a situação está representada na Fig. 11.8, 
Repare que nos dois casos devemos ter 


Ie= Lotta. (Eg. 11.1) 


11.4. GANHO DO TRANSISTOR DE JUNÇÃO 


Quando usamos um transistor em um circuito eletrônico, em geral estamos 
interessados em obter um ganho, isto é, em obter uma corrente ou tensão de 
saida maior que a corrente ou tensão de entrada. 

No caso de um transistor polarizado como na Fig. 11.7 ou 11.8, consi- 
derando o terminal de emissor como o terminal de entrada e o terminal de 
coletor como terminal de saída, temos 


Canto de porrénia sm SE E (Eq. 11.2) 
Vi le + dj 


Esse ganho de corrente é representado pelo símbolo q. 

Examinando a Eqg. 11.2, podemos ver que x é sempre menor que 1. Entre- 
tanto, se a base for muito estreita, quase toda a corrente de emissor chegará 
ao coletor e « será muito próximo de 1. Assim, por exemplo, se 98% dos por- 
tadores injetados pelo emissor chegarem ao coletor e apenas 2% se recom- 
binarem na base, teremos 


Surge imediatamente a pergunta: se a corrente de saida é sempre menor 
que a corrente de entrada, como podemos usar um transistor para amplificar 
um sinal? Existem duas respostas para essa pergunta: em primeiro lugar, não 
somos obrigados a considerar o emissor como entrada e o coletor como saida; 
existem outras maneiras de aplicar o sinal ao transistor, como veremos daqui 
a pouco, que permitem um ganho de corrente maior que 1. Em segundo lugar, 
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mesmo considerando o emissor como entrada e o coletor como saida, podemos 
obter um ganho de tensão e um ganho de potência maiores que a unidade. Vamos 
examinar essa segunda afirmação. 

Considerando apenas o transistor, sem levar em conta as resistências externas 
do circuito, a tensão de entrada é igual ao produto da corrente de emissor pela 
resistência da junção emissor-base e a tensão de saída é igual ao produto da 
corrente de coletor pela resistência da junção base-coletor. Assim, temos 

Tensão de saída 1cRgc R Roc. 
Tensão de entrada 1,Rs,g Ren 

onde « é o ganho de corrente, Rgç é à resistência da junção base-coletor e Rpg 
é a resistência da junção emissor-base. 

Ora, como a junção emissor-base está polarizada diretamente, sua resis- 
tência é pequena; como a junção base-coletor está polarizada inversamente, 
sua resistência é elevada. Em consegiência, a relação Rpc/R;g é muito maior 
que a unidade e portanto o ganho de tensão é elevado. O ganho de potência, 
que é dado pelo produto do ganho de tensão pelo ganho de corrente, também 
é elevado. Vejamos um exemplo numérico. 


11.5. EXEMPLO 


No circuito da Fig. 11.7, a resistência da junção emissor-base é 250 e a 
resistência da junção base-coletor é 50kQ. Sabemos que 98% dos buracos 
injetados na base chegam ao coletor. Calcule o ganho de corrente, o ganho de 
tensão e o ganho de potência do transistor. 

Solução: 

Ganho de corrente = a = 0,98: 


Ganho de tensão = 


Ganho de tensão = & Roc = 0,98 x 20009 
Ren 25 


Ganho de potência = Ganho de corrente x Ganho de tensão = 
= 0,98 x 1960 = 1920. 
11.6. FUNCIONAMENTO DO TRANSISTOR NA 
CONFIGURAÇÃO BASE COMUM 


Quando o sinal de entrada é aplicado ao emissor do transistor e o sinal 
de saida é colhido no coletor, como supusemos até agora, dizemos que o tran- 


= 1960; 


Ie =Is Ic=<Te=<Ts 


RL 


[6 Te= (1-=<) Is 
+ 


Figura 11.9 Ligação base comum 


SCE ei 
[É R 


Is 
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sistor está funcionando na configuração base comum. A Fig. 11.9, mostra as 
correntes em um circuito base comum. 

Quando um transistor funciona na configuração base comum, como na 
Fig. 11.9, a corrente de sinal (15) é aplicada ao emissor. A corrente de saida, 
alç, atravessa o resistor de carga, R,, produzindo uma queda de tensão que é 
a tensão de saida do amplificador. 

O ganho de corrente «, definido como 16/1,., é chamado de ganho de cor- 
rente base comum. 


11.7. FUNCIONAMENTO DE TRANSISTOR NA CONFIGURAÇÃO 
EMISSOR COMUM 


Como vimos no inicio deste capitulo, o ganho de corrente base comum, 
a, é muito próximo da unidade. Isso significa que a corrente de base é pequena 
em comparação com as correntes de emissor e coletor do transistor. Com efeito, 
combinando as Egs. 11.1 e 11.2, obtemos 


Ip=le-le= lp; =h(l-0o) (Eg. 11,3) 
Assim, se tomarmos 4 = 0,98, teremos 
1, = 0,02 1, = 0,0204 1, 


O fato de que a corrente de base de um transistor é muito menor que as 
correntes de emissor e coletor sugere o uso da base como terminal de entrada. 
Quando o sinal de entrada é aplicado à base e o sinal de saída é colhido no co- 
letor. dizemos que o transistor está funcionando na configuração emissor comum. 
A Fig. 11.10 mostra as correntes em um circuito emissor comum. 


RL 
— = IB 
Ic= = 


e 


Figura 11.10 Ligação emissor comum 


Quando um transistor funciona na configuração emissor comum como na 
Fig. 11.10. a corrente de sinal (1; = 15) é aplicada à base. A corrente de saida, 
Ie. atravessa o resistor de carga, R,. produzindo uma queda de tensão que é 
a tensão de saída do amplificador. Vamos calcular o valor dessa corrente de 
Saida em função da corrente de entrada. 

De acordo com as Egs. 11,1 e 11.2, temos 


e= [,-1,= Idi-lo. 
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Então 
le = [eu 


x 
Iç = 1-2 I; . 

Se definirmos agora o ganho de corrente emissor comum, B, como a relação 
entre a corrente de saida e a corrente de entrada nessa configuração teremos 


too Fu, 
da lg 1-a 
Tomando « = 0,98, temos 
Ê = 49. 


Podemos ver portanto que na configuração emissor comum O ganho de 
corrente pode ser bastante elevado. O ganho de tensão e o ganho de potência, 
que também são elevados, podem ser calculados do mesmo modo que na con- 
figuração base comum, conhecidas as resistências de entrada e saída do tran- 
sistor. 


11,8. TRANSISTORES DE DERIVA! 


Até agora supusemos que a distribuição de impurezas na base do transistor 
era uniforme. Nos transistores de deriva, a concentração de impurezas na base 
é muito maior em uma das extremidades da base do que na outra (Fig. 11.11), 


Emissor | Base Coletor 


No, Na Figura 11.11 Distribuição de im- 


purezas em um transistor de de- 


| 
| 
| 
| riva NPN 


Naturalmente, a variação na concentração de impurezas produz uma va- 
riação proporcional na concentração de portadores em maioria. Em conse- 
qiência da variação na concentração de portadores, existe uma corrente de 
difusão no interior da base. Entretanto nas condições de equilibrio, a corrente 
total na região da base deve ser zero, de modo que surge um campo elétrico 
no interior da base tal que a corrente de deriva seja igual e oposta à corrente 
de difusão. 

No caso de um transistor NPN, a corrente de difusão na base tende a 
acumular buracos nas vizinhanças da junção base-coletor. Assim, o sentido do 


“Em inglês, drift transistors 
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campo elétrico deve ser tal que a corrente de deriva dos buracos na base seja 
do coletor para o emissor. Ora, um campo elétrico com esse sentido acelera os 
elétrons injetados pelo emissor em direção ao coletor. Em conseqiiência, o 
tempo que os portadores injetados pelo emissor levam para atravessar a base 
é muito menor nos transistores de deriva do que nos transistores comuns, em que 
os portadores se movimentam apenas por difusão. Levando menos tempo para 
atravessarem a base, os portadores injetados têm menos tempo para se recom- 
binarem com os portadores da base. O resultado é que o ganho de corrente 


r 


dos transistores de deriva é muito elevado. 
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12 Outros fenômenos observados 
nos semicondutores 


12.1. EMISSÃO DE LUZ 


A recombinação de um par elétron-buraco envolve a liberação de uma 
certa quantidade de energia, equivalente à necessária para romper uma ligação 
covalente. Na maioria dos casos, essa energia assume a forma de vibrações 
mecânicas (emissão de fônons). Em determinadas circunstâncias, entretanto, a 
energia pode ser irradiada sob a forma de ondas luminosas (emissão de fótons). 

Assim, por exemplo, se submetermos um cristal de fosfeto de gálio a ra- 
diação ultravioleta, o cristal emitirá uma luz vermelho-alaranjada. O que acon- 
tece é que a radiação ultravioleta rompe muitas ligações covalentes no cristal, 
criando um excesso de pares elétron-buraco. Ao se recombinarem, os pares 
emitem uma radiação visível. 

O mesmo fenômeno também pode ser observado nas junções PN de certos 
semicondutores. Quando uma junção PN é polarizada diretamente, ocorre um 
grande número de recombinações nas proximidades da junção. Quando as 
condições são favoráveis, essas recombinações podem produzir luz visível. 

Normalmente, esse efeito não pode ser observado nos semicondutores mais 
comuns, O germânio e o silício, porque nesses materiais a grande maioria das 
recombinações resulta na emissão de fônons (vibrações mecânicas). Por outro 
lado, em alguns semicondutores compostos, como o fosfeto de gálio e o arse- 
nieto de gálio, o efeito é bem acentuado. 

Finalmente, é preciso observar que a luz pode ser emitida pelo cristal de 
forma incoerente — sem uma fase definida —, ou de forma coerente — com uma 
fase bem definida. No primeiro caso, o fenômeno é chamado de luminescência; 
no segundo, de efeito laser. 


12.2. EFEITO TÚNEL 


Quando a densidade de portadores nas regiões P e N de uma junção PN 
atinge o nivel de degeneração, a junção passa a apresentar novas propriedades. 
Vejamos como isso acontece. 

Na Fig. 12.1(a) vemos o diagrama de níveis de energia de uma junção PN 
degenerada'*, Nesse caso, o nível de Fermi está acima da base da faixa de con- 
dução do lado N da junção e abaixo do topo da faixa de valência do lado P. 
Os níveis ocupados estão indicados pela parte sombreada do desenho. 

Se aplicarmos à junção uma pequena polarização (negativa do lado N, 
positiva do lado P), haverá uma corrente composta de duas parcelas: a primeira 
corresponde à corrente de uma junção PN normal; a segunda é consequência 
do chamado efeito túnel. De acordo com a Mecânica quântica, uma partícula 
tem uma probabilidade finita de atravessar uma barreira de potencial, embora 
classicamentc não tenha cncrgia suficiente para vencer a barreira. Assim, com 


'*Para uma discussão do fenômeno da degeneração, veja Wannier, G. H., Elements of Solid 
State Theory (Londres: Cambridge University Press, 1960), p. 149 
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Região P Região N 


É to de Fermi 


Não hã 
Efeito Túnel 


(a) (e) 


(b) 


Figura 12.1 Efeito túnel em uma junção PN degenerada. (a) Junção em equilibrio; (b) 
Junção submetida a uma pequena polarização direta, que produz o efeito túnel; (c) Junção 
submetida a uma polarização direta mais elevada. O efeito túnel deixou de existir 


uma pequena polarização direta, os elétrons podem passar da região N para 
a região P, ocupando os estados vagos da faixa de valência da região P, como 
na Fig. 12.1(b). Entretanto se aumentarmos a polarização de maneira tal que 
a base da faixa de condução do lado N tenha energia superior ao topo da faixa 
de valência do lado P, como na Fig. 12.1(c), a corrente devido ao efeito túnel 
será reduzida a zero, pois não haverá estados disponíveis para os elétrons do 
lado P com a mesma energia que eles possuíam do lado N (no efeito túnel, a 
energia é conservada). Em consegiiência, a corrente total diminuirá. Se aumen- 
tarmos ainda mais a polarização, a corrente tornará a aumentar, devido ao 
aumento da corrente normal da junção. 

Assim, a caracteristica corrente-tensão de uma junção PN degenerada po- 
larizada diretamente tem a forma que vemos na Fig. 12.2. Repare que a curva 
apresenta uma região em que a resistência é negativa. 


Figura 12.2 Curva característica 5 
corrente-tensão de um diodo tú- a 
nel 5 io 
o Região de 
ade Resistência 
| | Negativa 


TENSÃO 
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O diodo túnel, como é chamado o dispositivo baseado nesse efeito, pode 
ser usado como elemento ativo nos circuitos, por exemplo, como oscilador de 
microondas. 

Os diodos túnel em peral são feitos de germânio, silício ou arsenieto de 
gálio. 


12.3. EFEITO HALL 


Suponhamos que um cristal seja submetido simultaneamente a um campo 
elétrico e a um campo magnético, perpendiculares entre si (Fig. 12,3). 


Figura 12,3 Geometria do efeito 
vm Hall 


V 


O campo elétrico pode ser aplicado, por exemplo, ligando-se uma fonte 
de tensão V às extremidades do cristal. O campo magnético pode ser aplicado 
colocando-se o cristal entre os pólos de um eletroíimã. 

Sob a ação do campo elétrico, os portadores do cristal adquirem uma certa 
velocidade de deriva, v. Ora, sabemos do eletromagnetismo que uma partícula 
carregada sob a ação de um campo magnético B está sujeita a uma força per- 
pendicular à direção do movimento e à direção do campo magnético, Essa força 
é dada pela equação vetorial 


F, = (o x B), (Eq. 12.1) 


onde q é a carga da partícula. 
A Fig. 12.4(a) mostra o que acontece quando o portador é um elétron. 
Quando o portador é um buraco, a situação está representada na Fig. 12.4(b). 


Feeeeeeeerrerrerereo 


(a) 


E a il o RR alo a a SR a o a a 


Figura 12.4 Efeito Hall para os dois tipos de portadores. (a) elétrons. (b) buracos. O 
sentido do campo magnético é para dentro do papel 
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Sob a ação do campo magnético, os portadores se deslocam para uma 
das faces do cristal. Surge então um campo elétrico transversal, e, produzido 
pelo desequilíbrio de cargas. Em equilíbrio, a força aplicada ao portador pelo 
campo magnético, F,, é igual à força aplicada pelo campo elétrico transversal, 
F,. Para a geometria das Figs. 12.3 e 12.4, podemos escrever 


F,, = quB; (Eq. 12.2) 
F. = QEg- (Eg. 12.3) 

Assim, em equilíbrio temos Í 
vB = Ey- (Eg. 12.4) 


Por outro lado, chamando de J a densidade de corrente e de a a concentração 
de portadores, temos 


J = qav. (Eq. 12.5) 
Assim, temos (Fig. 12.3) 


pes —— = —= (Eq. 12.6) 


Em conseqiiência, temos 


Eg = vB = (Eq. 12.7) 


qaht 


A diferença de potencial que aparece entre a face superior e a face inferior 
do cristal é dada por 


ou 
Va = Ru? (Eq. 12.8) 


onde Rg = 1/ga é chamado de constante de Hall do material considerado. 
Quando os portadores são todos elétrons (carga —q, concentração n), po- 
demos escrever 


Ry= e (Eg. 12.9) 
Quando os portadores são todos buracos, temos 
1 
=—. (Eg. 12.10) 
"q q 


As Egs. 12.9 e 12.10 mostram que a polaridade da tensão V, (que é cha- 
mada de tensão de Hall) depende do tipo de portador presente no cristal, O 
efeito Hall pode ser usado portanto para determinar se um semicondutor é 
tipo N ou tipo P. 
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EXERCÍCIOS 


'*Considere a experiência esquematizada na figura a seguir. 


Sabendo que 


Va =0,5V; 

1 = IOmA; 
V=-1Y; 

B =2x10"3T. 
Pede-se 


a) o tipo dos portadores (elétrons ou buracos), 
b) a resistividade do material, 

c) a concentração dos portadores, 

d) a mobilidade dos portadores. 
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APÊNDICE A Teorema de Bloch 


O teorema de Bloch diz que a solução mais geral da equação de Schroedinger 
para um elétron em um cristal é 


Val) = e Tur), (Eq. A.1) 
onde uz(7) é uma função com a mesma periodicidade espacial que a rede cris- 
talina. 

Prova: Consideremos o operador T;, definido por 

TEM = Fr +R), (Eq. A.2) 

onde R; é o vetor de translação da rede, isto é 
Vi) = Ve + Rj. 
Aplicando o operador T; à Eq. 5.12, temos 


h»T Go q? 0? o «sá É 
las +ãa + 22] vs; + vw) = ET) 
como 
TAVE 2] = VIP + RITA) = VOTA lo). 
temos 


E vita] + VOLTA] = EMT. 


Assim, a função Tab) satisfaz à mesma equação que a função Welr) (Eq. 
5.12), Nesse caso, T.W-z(F) só pode diferir de Wy ;(r) por um fator constante. Vamos 
chamar esse fator constante de a,. Assim, temos 


Tg) = balé + R) = ay; (Eq. 4.3) 


Vamos fazer 

a; = et b BR. (Eq. A4) 
onde k é um vetor complexo. Trata-se de uma maneira perfeitamente geral de 
escrever a. 


Consideremos agora a aplicação sucessiva de dois operadores T: 


T Ty lr) = Tal”) = asa bl”) = dê (Rot Romy (7), (Eq. A.5) 


Como a aplicação sucessiva de dois operadores T corresponde à translação 
de um vetor igual à soma dos vetores de translação individuais, a Eg. A.5 mostra 
por que escolhemos para a, a forma da Eq. A.4. Vemos também que é preciso 
que k seja um vetor real, caso contrário a função de onda cresceria indefinida- 
mente com a distância ou seria atenuada rapidamente. 
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Vamos agora definir uma função Hg(P) do seguinte modo: 
nl) = emp). (Eq. A.6) 
Usando as Egs. A.2, A3 e A.4, temos 
Tu) = ur + R)= Te pad] = 
Fe eTikir+ *Tur() == e RARO RR p (7) == 
= np?) = po(r). (Eg. A.7) 
Assim, a função g(r) tem a periodicidade da rede cristalina, isto é, 
nl” + R$) = gy(?). Multiplicando os dois membros da Eq. A.6 por exp (ik 7), 
temos finalmente 
Val) = ul) (Eq. A.8) 


o que completa a prova do teorema de Bloch. 


APÊNDICEB Experiências em semicondutores 


O desenvolvimento da física dos semicondutores deve-se sem dúvida a um 
exaustivo trabalho de pesquisa experimental, que prossegue até hoje em labo- 
ratórios de todo o mundo. 

Julgamos importante que o estudante se familiarize desde cedo com as 
técnicas experimentais, utilizadas para investigar o comportamento dos semi- 
condutores. Por esse motivo, decidimos apresentar neste apêndice uma lista de 
referências a experiências simples em semicondutores. Essas experiências po- 
derão ser usadas, a critério do professor, para ilustrar e/ou expandir os conceitos 
teóricos desenvolvidos nesta obra. 
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EXPERIÊNCIAS 


Experiência 1. Ilustração do crescimento de um cristal utilizando blocos cúbicos de ma- 
deira ou papel. Ref. 1. p. 58. 

Experiência 2, Ilustração da estrutura de cristais utilizando bolas furadas e hastes ci- 
lindricas. Ref. 1. p. 59. 

Experiência 3. Determinação da frequência critica, f., na qual um cristal começa a absorver 
encrgia de um feixe de luz monocromática (photoelectric threshold). Ref. 1. p. 92. 

Experiência 4. Estudo do fenômeno da fotocondutividade, ou seja, da variação da con- 
dutividade de um material semicondutor em função do fluxo luminoso incidente. 
Ref. 1. p. 93. 

Experiência 5. Ilustração do fenômeno de luminescência em um cristal. Ref 1. p. 94. 

Experiência 6. Verificação do tipo de semicondutor (P ou N) para materiais fortemente 
dopados. Toca-se o cristal com uma ponta aquecida e observa-se a leitura de um 
voltimetro de alta impedância ligado entre essa ponta e um outro ponto do cristal. 
Ref. 1. p. 87. 

Experiência 7. Experiência de Shockley-Haynes. Essa experiência permite medir as cor- 
rentes de deriva e difusão e verificar experimentalmente a relação de Einstein (D/u = 
= kT/g). Ref. 1. p. 120, 143. 

Experiência 8. Medida da condutividade de uma amostra de germânio intrínseco, em 
função da temperatura. Essa experiência permite medir o coeficiente de temperatura 
de q e também estimar a energia necessária para criar um par elétron-buraco. Ref. 
1. p. 143. 

Experiência 9. Efeito Hall — Experiência extremamente interessante, que além de suas 
aplicações industriais (medida de corrente), permite uma confirmação da existência 
de dois tipos de portadores de carga. Ref. 1. p. 126e Ref. 2. pp. 195, 205. 
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Experiência 10. Ilustração dos fenômenos de recombinação e tempo de vida de portadores. 
Ref. 1. p. 147. 

Experiência 11. Medida da capacitância de uma junção PN em função da tensão inversa 
aplicada. Ref. 1. p. 180. 

Experiência 12. Curva caracteristica de um diodo — Experiência que tem como fina- 
lidade levantar experimentalmente a curva característica de um diodo. Rel. 3. pp. 
38, 39e 40. 

Experiência 13. Medida da condutividade de um semicondutor em função da temperatura 
(experiência mais elaborada que a experiência 8). Ref. 2. p. 198. 

Experiência 14. Produção de pares elétron-buraco e largura da faixa proibida de energias 
(energy gap). Ref. 2. p. 213. 

Experiência 15. Medida da mobilidade dos portadores em um material semicondutor. 
Ref. 2. p. 222. 

Experiência 16. Construção de um diodo tipo liga (ulloy). Ref. 2, p. 230. 


Respostas e sugestões para 
os exercícios propostos 


CAPÍTULO 1 
1.1. ay = 5,292 x 107 !!m. 
1.2. Use a Fig 1.5 ca Eg. 1.10. 


CAPÍTULO 2 


2.5. É possível provar que a seguinte relação é válida para um polinômio de Hermite 
de ordem n: 
[ug dê =Y-niz 


CAPÍTULO 4 


41. a) 5,8 x 1072 
b) 407 À. 
4.4. 109º 28". 
4.7. Observe que [100], [010] e [001] são direções de “mesmo tipo”. 


O plano (011) contém as direções [100]. [011] e [111]. 
O plano (101) contém as direções [010], [101] e [111]. 
O plano (110) contém as direções [001], [110] e [111]. 


CAPÍTULO 7 


7.3. Considere o recipiente onde está o gás suficientemente raso e portanto despreze a 
variação da energia potencial. Nesse caso 


2 2 2 
Misa 23 0270 BP +. 
E 3 (+ o, + uv) spa 
CAPÍTULO 8 
-No + No + nº 
Bio à) fe SEA 
bj=10%. 


8.6. Use o método gráfico de Schockley. 


a) Nível de Fermi cerca de 8kT acima de E,. 
b)n=2x104cm. 
e) Desprezivel. 


CAPÍTULO 9 


94. 0, =4x 10º cm/s. 
D = 10 cm?/s. 


124 Introdução à física dos semicondutores 


9.6. Suponha N, > n, e verifique depois se a suposição é correta. 
p=33x 10'*em'. 
n=1,7 x 10!2em”3. 

9.8. Prove que 


do qEs . 


dr” “2kT 


9.9. Prove primeiro que 
QE tor aituto 
HE) é — qrestimet, e ' 
Ao esa Le =D ] 


CAPÍTULO 10 


10.1. a) 5,0 x 10 “ em. 
b) 9.9 x 107º cm. 
10.3. a) V, = 0017 V; 
V, = 9,983 V. 
b) 10 ua. 


CAPÍTULO 12 


a) Elétrons. 

b) 1Q cm, 

c) 6,2 x 10! cm 2. 
d) 10º cm?/V's. 
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